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ABSTRACT 


During  this  report  period,  work  nas  continued  on  a  survey  cf  the 
structure,  electrical  and  magnetic  properties  of  transition  metal  oxide- 
phosphate  glasses  and  glasses  in  the  As^e^-As^e^  system  which  possess 
electrical  or  magnetic  device  potential.  Results  of  magnetic  and  elec¬ 
trical  observations  in  several  transition  metal -phosphate  glasses  have 
revealed  a  high  degree  of  magnetic  and  structural  order.  The  pronounced 
influence  of  glass-glass  phase  separation  and  compositional  segregation 
has  also  been  noted. 

Examination  of  a  series  of  FeO-P^O^  glasses  has  revealed  a  glass 
forming  range  which  extends  to  80  mole  %  FeO.  Low  temperature  suscep¬ 
tibility  and  MOssbauer  spectroscopy  results  on  these  glasses  are  reported 
herein.  Low  temperature  susceptibility  of  the  manganese,  vanadium  ano 
iron  glasses  indicate  that  these  glasses  are  all  amorphous  anti  ferro¬ 
magnetic  systems  with  a  distinctive  downward  curvature  of  the  recipro¬ 
cal  susceptibility  versus  temperature  plot. 

Detailed  studies  of  the  As2Tej-As2S6j  system  has  shown  switching 

behavior  which  can  be  controlled  with  compositional  variation.  The 

compositional  variation  of  the  microstructure  Is  noted  in  micrographs 

i  nc  I  uded  in  this  report .  Defoils  of  illustrations  in 

iKis  document  may  bo  bottor 
Hu. lied  on  microfithe. 
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ol  relationships  between  glass  prapara.lon  v.rl.bl„  tM 

trlcal  and  magnetic  properties. 

Electronic  conduction  In  amorphous  solids  has  become  tha  subject  of 
Interest  to  a  number  of  th.or.tlcl.ns  end  hes  bean  reviewed  by  Mott'. 2 
Gubanov3  and  numerous  others.  Virtually  ell  0f  these  works  have  begun 
with  an  assumption  that  emorphous  solid,  ere  uniformly  rendom.  even  though 
they  recognize  glasses  are  ganarally  heterogenous.  These  theoreticians 
have  developed  analytical  descriptions  of  several  systems  which  have  been 
experimental ly  verified  In  s™  cases.  Attempts  to  extend  this  approach 
to  microscopically  heterogenous  systems  have  bed  notably  little  success. 
There  remains  a  considerable  body  of  experimental  results.  Including 

H°"  3nd  SMb8Ck  C08ff are  not  rationalized  by  present  theory. 

Pearson4  has  suggested  that  heterogenous  structure  In  these  materials 
may  explain  these  anomalies  If  the  separated  phase  Is  crystalline.  It 
appears  that  heterogenous  transport  analysis  similar  to  that  of  Volger3 
or  Bube*  Is  repaired  to  ascertain  the  transport  behavior  In  each  phase. 

Another  Important  anomaly  between  theory  and  observation  concerns  the 
theoretically  predicted  Insensitivity  of  amorphous  semiconductors  to  doping. 
Early  experimental  observations  by  Kollmets.  et  el. 7,  conformed  to  the 
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Farther  *vid#»»c*  tMt  ktructural  retbrogeMtltk  !*•  •*  ,rv#  f Wif  °* 
thti#  inomllM  can  t«  Inferred  from  «ork  by  Klnwr,  of  «l.f.  In 
V^glesse*.  This  work  Mi  ir^n  that  merked  changti  In  dlalactrlc  bahavlor 
occur  during  thermal  treatments  customarily  used  to  strait  rallava  glasses. 
Tnese  changes  have  been  shown  to  be  the  result  of  structural  changes  Involving 
precipitation  of  small  amounts  of  crystals. 

Wilson  and  Klnser10  have  observed  similar,  but  somewhat  more  complex, 
behavior  in  Fe0-P205  glasses  after  thermal  treatments  corresponding  to 
annealing.  Electron  spin  resonance  (ESR)  results  have  shown  the  onset  of 

structural  changes  during  thermal  treatment  prior  to  their  observation  by 

11  12 

other  commonly  employed  techniques  » 

it  is  thus  apparent  that  homogenous  glasses,  semiconducting  or  otherwise, 
are  the  exception  rather  than  the  rule. 

Genera i  Methodoiogy 

The  electrical  and  magnetic  property  changes  accompanying  structural 
modifications  during  glass  processing  are  of  prime  interest  in  the  present 
work.  The  above  questions  can  only  be  answered  with  detailed  structural 
characterization  of  representative  glasses  from  the  oxide  and  chai cogen ide 
groups.  The  initial  oxide  glasses  examined  were  the  Fe0-P205,  V205"P205' 
Cu0-P205,  Ti02-P205  and  MnO-P^  systems.  The  chalcogenide  glasses  are  from 
the  As2Te3_As2Se3  system  with  Ag-As-S  glasses  in  preliminary  stages  of  study. 
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t»«  c  «t»  *ftc  cn'r1  •*«  ion  <e*s#*t  •#•  •ii#*?J*iiir  I  **•#•**•*»  of 
*o*iP#rel«iro  0*0f  I  to  frocvoecy  re*?*  500  fo  1 1,  COO  fer  infrerod 

fTM^rononit  c<yt«wc!»3  #t  557  wlcroe*  indicate  further  in#!  the  tren¬ 
ail, ion  constant  lor  tr»«  80  arsenic  trltellurlee  ?0  arsenic  trlsenlenlde 
glass  1 4  Independent  of  temperature  over  the  rang*  from  10  to  300«K.  Tho 
jhF  dielectric  properties  ara  compositional ly  dependant  but  ppear  to  be 
temperature  Insensitive.  These  results  are  discussed  In  further  detail 
In  the  attached  thesis  of  Mr.  J.  0.  Pearson. 

TRANSITION-METAL  PHOSPHATE  GLASSES 
Manganese-phosphate  glosses 

Magnetic  behavior  of  two-phase  manganese  phosphate  glasses  was  studied 
by  means  of  variable  temperature  magnetic  susceptibility  and  electron  spin 
resonance  techniques.  The  high  temperature  magnetic  susceptibility  follows 
a  Cur  I  e-Weiss  law.  Figure  one  shows  that  the  low  temperature  recipro¬ 
cal  susceptibility  bends  downward  Indicating  considerable  influence  by  para¬ 
magnetic  ions.  The  electrical  conductivity  of  the  glass  is  extremely  low 
suggesting  that  almost  all  of  the  manganese  ions  are  in  the  divalent  state. 

It  has  been  determined  that  the  manganese  divalent  ions  exist  in  the  glass 
in  one  of  two  distinct  magnetic  states-antiferromagnetical ly  coupled  pairs 
in  the  manganese  rich  phase,  or  Isolated  ions  in  the  phosphate  rich  phase. 
Magnetic  behavior  was  studied  as  a  function  of  glass  phase  separation  and 
as  a  function  of  glass  composition  over  the  glass  forming  range.  The 
results  were  interpreted  in  terms  of  a  model  consisting  of  paired  and  isolated 
divalent  manganese  ions.  These  results  were  reported  at  the  May  American 
Ceramic  Society  meeting  and  the  First  International  Amorphous  Magnetism 


A  cow  of  f*#  mfrom  if# 
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A  gla*t  forming  study  of  th#  Iron  oxid#  phosphate  tyiiem  has  shown 
tnat  th#  glat»  forming  rang#  extends  to  80  mole  percent  ferric  oxide  when 
fne  glasses  are  melted  In  an  oxidizing  atmosphere  end  quenched  In  smell 
quantities.  Glasses  above  80  mole  percent  ferric  oxide  spontaneously 
crystallized  to  alpha  ferric  oxide  while  Intentionally  devltrlfled  glasses 
of  lower  Iron  content  crystallized  to  Iron  phosphate.  The  magnetic  suscep¬ 
tibility  of  glasses  examined  was  described  by  a  Curle-Welss  law  over  the 
range  100  to  525°K.  Susceptibility  studies  Indicated  that  the  majority 
of  the  Iron  Ions  are  antlferromagnetlcal ly  ordered  over  the  temperature 
and  compositional  ranges  examined.  Below  I00°K,  the  reciprocal  suscep¬ 
tibility  bends  downward  and  at  approximately  I0°K  the  susceptibility 
becomes  constant  as  shown  In  Figure  two.  These  results  have  been  reported 
In  two  recent  papers.  Copies  of  the  abstracts  of  these  papers  are  appended. 

The  coordination  and  magnetic  structure  of  the  55  mole  percent  ferric 
oxide  glass  containing  various  concentrations  of  ferrous  and  ferric  ions 
were  studied  using  Mossbauer  effect  spectroscopy  at  77  and  300°K.  Room 
temperature  Isomer  shifts  In  the  range  0.107-0.110  cm/sec  and  0.023-0.047 
cm/sec  and  quadruuols  splittings  of  0.217-0.265  cm/sec  and  0.046-0.097 
cm/sec  for  the  ferrous  and  ferric  Ion,  respectively,  Indicate  tetrahedral 
coordination  states. 

Although  no  magnetic  hyperflne  structure  was  present  at  300°K, 
splitting  at  77°K  was  interpreted  as  evidence  of  either  long  relaxation 
times  or  limited  magnetic  ordering.  Mossbauer  spectroscopy  of  samples 
hear  treated  for  intervals  up  to  six  hours  at  600°C  showed  progressive 
precipitation  of  oxides  In  several  crystalline  phases.  A  detailed  dis- 
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Tr»  eiecirldal.  thermal  ..pension  a*  mlcroatructur.l  chareeiarlailc, 
at  vanadium  phosphate  glasses  hlv*  **•*  examined  •*  *  'unction  ot  ““P0* 

„<1  oxidation  atait.  Thaw  glataa.  axllalt  liquid  -  liquid  I»ikIoM'*V 
car  a  .Id.  rang.  o.  vanadium  rich  coalition,  a*  Xlghar  ot 

T„.  presence  ot  tha  miscibility  I,  also  fatlactwl  In  the  al^trlcal  and 
softening  point  observations.  The  compositional  and  oxidation  dep.nd.nc.  ot 
the  electrical  propartl.s  and  microstructur.  Is  bast  discussed  In  terms  of  . 
ternary  diagram  consisting  ot  t.o  different  vanadium  oxides  and  phosphorous 

pendoxlde.  men  th.  system  Is  an . .  from  this  point  ot  vie.  the  electrical 

conductivity  behavior  I . to  be  clearly  mlcrostructur.l ly  dep.nd.ht  .Ith 

a  valley  ot  minimum  conductivity  nearly  corresponding  to  a  splnodal  composition 
aoross  the  ternary  diagram.  This  .ork  .as  presented  at  the  May  American  Ceramic 
Society  and  Is  presently  being  prepared  for  publication. 


CHAICOUNIK  OiASStS 


fhe  Plcroftructur*  end  OC  and  AC  electrical  properties  of  •  series  of 
arsenic  triselenlde-ersenlc  trltel lurid*  glasses  **ve  bNn  investigated. 
Electron  microscopic  observations  Indicate  th*  pr*t*nc*  of  phase  separation 
In  these  glasses.  A  metestebl*  miscibility  gap  across  the  quasl-blnary/ 
diagram  Is  Indicated.  Variations  In  DC  electrical  properties  may  be  due 
to  structural  changes,  compositional  changes  or  a  combination  of  the  two. 
Structural  heterogen  1 1 1 *s  due  to  rhase  separation  appear  to  be  responsible 
for  Maxwel l-Wagner-SI I lars  heterogeneous  losses  In  these  glasses.  These 
results  are  discussed  In  further  detail  In  the  attached  thesis  of  Mr.  0.  Hill 
and  an  abstract  of  a  paper  presented  at  the  May  American  Ceramic  Society 

meeting. 

Extensive  bulk  switching  studies  have  been  performed  on  the  same  group 
of  glasses  mentioned  above.  Both  threshold  and  memory  switching  have  been 
observed  in  these  glasses.  Breakdown  voltage  was  independant  of  sample 
thickness.  Switching  delay  times  at  iow  over-voitages  were  in  the  hundreds 
of  milliseconds  range.  Results  of  breakdown  voltage  vs  temperature  studies 
were  computer  fit  to  an  approximate  solution  of  the  heat  flow  equation 
appropriate  for  the  geometry,  yielding  experimental  values  of  activation 
energy  for  conduction  and  thermal  conductivity.  These  resuits  clearly  sub¬ 
stantiate  the  thermal  mechanism  for  electrical  switching  behavior  in  bulk 
samples  of  these  glasses.  A  recently  submitted  paper  detailing  these 
experiments  Is  attached  with  an  abstract  from  a  presentation  at  the  May 
ceramic  society  meeting. 

The  UHF  and  microwave  dielectric  properties  dielectric  constant  and 
loss  tangent  of  these  glasses  have  been  examined  as  a  function  of  tempera- 


SincsuraJ  ctarfuarifavio*  at  fhtsa  «**••**  It  t#iAg  accawfi I %,•** 
atif»g  DicfroA  nicroscooy,  *-ray,  •iKiromi  if  in  i-as  ***# 

iportroscopy,  «agna?!c  tuftctpf Ibl I Ify,  aiattro *  aicroproba,  Clalttvric 
r«i«n«t io«*,  m%t\ finer  spectroscopy  and  diffarantial  fhar*al  analysis  fee* 

niQuo*. 

In  conjunction  with  th#  structural  tools,  It  It  nacassary  that  tha 
conductivity,  switching  bahavlor  and  Saaback  coafflclant  by  sonltorad  to 
allow  dlract  structura-proparty  corralatlons. 
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I#  In  ^*54*1  eif«^xi««  is  •n»«  ii#  #*  mm*  #•  tm*  f#Mii  l« 

•  wmltief  »iwt  #*«'*§  I  I»*l  *4lCft  *r#  cle#r  »*.  i<#  If* 

?.  A»  In  our  previous  reco*w**A#t l<m,  *«  wt ir*e  »©  fdlfwmmil  is*# 
turv«y  preparation  of  new  glasses,  «*  anticipate  VMt  I*#  »re*%  Horn  **lel 
ox « de-phosphate,  silicate,  borate  end  germinate  survey  presently  in  ;r^r#n 
will  be  continued. 

It  is  also  anticipated  that  results  on  a  new  system  Ay-As-S  will  be 
most  helpful  In  developing  switching  models.  A  detailed  ternary  phase  di¬ 
gram  was  recently  published  (30)  and  our  analyses  should  be  simplified  with 
this  as  a  basis,  it  Is  further  anticipated  that  Cu  and  Au  substituted  In 
the  above  system  will  be  q:iite  Informative  from  an  atomistic  and  mlcrottruc- 
turai  model  point  of  view. 

3.  The  Mdssbauer  studies  should  be  continued  to  exmine  57  Fe  and 

127  fe  in  each  of  the  systems  presently  under  examination  using  other  tech¬ 
niques.  This  will  significantly  aid  In  atomic  structure  model  development 
in  these  systems  as  an  addition  to  the  present  toots. 

4.  We  recommend  that  the  far  Infrared  "conductivity  loss  spectra.  This 
will  allow  the  loss  behavior  to  be  explicitly  attributed  (31 )  to  each 
mechanism  thus  reinforcing  both  atomic  and  mlcrostructurai  analyses. 
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Fig.  2.  (fc)  Low  Ter-perature  Behavior  of  the  Magnetic  Susceptibility 
of  55  Mol  %  Fe0  45  Mol  %  P205  Glass 
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%  l-O-lt.  itrfki  of  tlM  F»iOrPA  Giw  System 
f  J.  H.  Dayani.*  D.  L.  Kink*  and  L.  K.  Wiuon,  Vanderbilt 
University,  Nashville,  Tenn. 

A  class  forming  study  of  the  Fe,0,-P.0»  system  has  shown 
that  the  glass  forming  range  exteads  up  to  80  mo!e%  Fe,0. 
when  the  glasses  arc  melted  in  an  oxidizing  atmosphere.  All 

*  glasses  were  melted  in  sir  using  a  platinum  ribbon  furnace  and 
rapidly  quenched.  Glasses  above  80  mo!e%  Fe.O,  spontane¬ 
ously  crystallized  to  «-Fe.O,  while  intentionally  devitriAed 
glasses  of  composition  below  80%  crystallized  to  FePO,.  The 
magnetic  susceptibility  of  all  glasses  examined  was  described 

*  by  a  Curle-Weiss  law  over  the  range  100*  to  5i  >*K.  Suscepti¬ 
bility  studies  Indicated  that  the  majority  of  the  Iron  ions  arc 
antiferromagnetically  ordered  over  the  temperature  and  com¬ 
positional  range  examined.  (Research  sponsored  by  the  U.S. 
Army  Research  Offlce-Durham). 


till 8-1  i:M  Mt 

t-SJ-72.  Thermal  and  Dependanaa  ef  Hectrical 

Switching  In  Ar.Te  -AsSe.  Glasses 

II.  R.  Sanders.*  L.  K  Wilson.  D.  J  Hill  and  D  L.  Kinscr. 
Vanderbilt  L’niverjity,  Nashville.  Tenn 
Kxtcnsive  bulx  switching  studies  have  been  performed  on 
the  xAs  Te.-f  I  -x)  As.Sei  system  as  a  f  jnction  of  composition 
Both  threshold  and  memory  switching  have  been  observed  m 
all  samples.  Breakdown  voltage  arid  suitchi.r;  subduv  in 
crease  with  increasing  Se  content.  Urrahdown  voltar-  **s 
etaentlally  Independent  of  sample  thickness  Switching. k-li 
tim*s  at  low  overvoltage  were  in  the  hundndt  o'  milltoi  .  s !« 
rang.!.  Result*  of  breakdown  voltage  vs  temprrature  stud.)  * 
were  computer  tit  to  an  approximate  wiutiOA  of  the  heat  floa 
equation,  yieldiiut  experimental  values  of  activation  energy 
tor  conduction.  These  results  o<  srly  substantiate  the  therms 
mechanisms  for  electrical  switching  behavior  in  bulk  samples 
of  these  glasses. 
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M  a  function  . if  composition  These  glasses  .shif  t 
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Of  this  structural  segregate  on  the  dc  rtwttlvgy  at  <eeu.ii,. 
lemperature  at  a  (unction  of  V* TV*  in  thr  glass  u  d«cuv«d 
r^rttam?  by  m i  UJ.  Army  Wsasarcti  otf.ee- 


fiW>t«il8  an. 

r  *-**-**•  Mirreetrnrtosal  and  Rieetrical  Pnpeniea  at 
AsHMsle,  Clssen 

O.f  Hut.*  I-  K  W  i  son  II  K  Saws  its  and  U.  I..  Kink*. 
VandcrtM't  I’mvtr  .|v,  Na«hvi.!e.  Tenn 

El*ct*na  mtctriK»*,v,  dc  and  at  c.uviactisitv  ohservatiicitt 
M'e  been  eood.  ;nl  nn  a  xeilcs  of  f  i>  Asv  guw 
ThtK*  r»  Mills  indicate  the  presence  of  extensive  pha*e  upin 
POA  m  alt  glavie*  rx  inw*  .1,  iectudet?  tu-  pun-  •  v- 
The  .4x11X21  fri'fWitnw  e»J  ihw  wgfc  n  xt.ctruei.r*. 

whkii  are  scns.i.*  t.i  thermal  h.w-,»f>  ate  CS.'iscvt*  m 
presuxuly  p,.,uiirl  rcuiUs.  TV  present  res.its  md*a.s  a 
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* 


IV 

ABSTRACT 
Submitted  to 

International  Symposium  on  Amorphous  Magnetism 

ANTI  FERROMAGNETISM  IN  THE  VANADIUM,  MANGANESE 
AND  IRON  PHOSPHATE  GLASS  SYSTEMS 

by 

L.  X.  Wilson,  E.  J.  Friebele  and  D.  L.  Kinscr 
Vanderbilt  University 
Vashville,  Tennessee  37235 

Variable  temperature  magnetic  susceptibility  and  magnetic 
resonance  measurements  have  been  made  on  a  series  of  concentrated 
vanadium,  manganese  and  iron  phosphate  glasses.  The  high  tempera*  . 
ture  magnetic  susceptibility  of  all  glasses  studied  obeyed  a 
Curie* eeiis  law  with  an  anti ferromagnetic  Curie  temperature.  At 
low  temperatures,  the  magnetic  susceptibility  data  have  been 
interpreted  in  terms  of  recent  theoretical  models  of  amorphous 
anti ferromagnetism  (1,2). 

Measurements  of  magnetic  resonance  line  intensity  and  line- 
width  at  9.0  GHz  as  a  function  of  temperature  have  shown  the 
existence  of  extensive  anti  ferromagnetic  coupling  between  the 
transition  metal  ions  (3,4).  The  resonance  measurements  have 
been  shown  to  be  sensitive  to  micros true tural  inhomogeneities 
(liquid-liquid  phase  separation  and  devitrification)  in  the 
glasses.  Weak  transitions  at  temperatures  corresponding  to  the 
Neel  temperatures  of  crystalline  transition  metal  oxides  have 
been  observed. 
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ELECTRICAL  SWITCHING  PHENOMENA  IN  THE 
A*2Te3-A*2S«3  GLASS  SYSTEM 


This  paper  reports  the  results  of  electrical  switching  studies  of  a  series  of 
As.  T«3-Aa2Se3  glasses.  Kolomeits  and  Nagsrovs  (1)  and  Rolios  (2),  and  more  recently 
Hill  (3),  have  reported  electrical  conductivity  stud.es  of  this  system. 

Rinser,  at  al.  (4)  have  also  reported  electrical  property  studies  in  connection  with 
a  study  of  the  aicrostructuta  of  this  systea. 

All  of  the  glasses  studied  exhibit  bulk  threshold  and  aeaory  switching  with  the 
critical  voltage  for  switching  (breakdown  voltage  V^)  increasing  with  At.  S#  content. 
The  switching  process  observed  in  these  glasses  has  been  found  to  occur  by  a  thermal 
aachanlaa  in  agreeaent  with  studies  previously  reported  by  Warren  (S)  for  a  single 
coaposltion  of  aaorphous  As  Se^Te  ind  Tanaka  and  his  co-workers  (6-10)  in  the  .eoronous 
Aa-Te-Ge  svstea.  The  theraul  dependency  of  the  breakdown  voltage  has  been  fitted  to  a 
thermal  model  in  agreement  with  the  theoretical  work  of  Boer  and  Ovshlnsky  (11)  and 
Sheng  and  Vcstgate  (12).  In  thesw  papers,  an  approximate  solution  of  the  heat  flew 
equation  is  used  to  derive  an  equation  for  the  temperature  dependence  of  breakdown 
voltage: 


VT> 


<i) 


where 


*  •  thermal  conductivity 
AE  ■  activation  energy  of  conduction  and 


[£E  V  1 

■n*  rj 

where  exp  expresses  the  field  dependency  of  conductivity.  The  effect  of  the 

field  dependent  term  is  negllgable  for  bulk  samples  although  it  becomes  increasingly 


1G 


2 


core  significant  as  the  thickness  of  the  sar.ple  decreases.  Although  tae  results  of 
t  iese  studies  agreed  well  in  all  cases  with  a  thermal  description,  no  attempt  will  be 
made  to  generalize  these  results  to  thin  films,  for  which  there  is  still  controversy 
about  the  nature  of  tha  breakdown  mechanism.  Stocker  (13)  and  Warren  (14)  have  argued 
tnat  the  thermal  description  still  aoplies  tor  thin  films,  while  otners  such  as  Boer  and 
Ovshinaky  (11)  maintain  that  at  least  in  some  cases,  other  affects  are  predominant. 

Experimental 

All  glares  examined  ir  this  study  were  prepared  by  fusing  tr.e  appropriate 
materials  in  evacuated  Vycor  ampoules  at  800*C  in  a  rocking  turnanc*.  After  hasting 
for  one  hour,  the  ampules  wert  rapidly  quenched  in  water.  Reagent  grade  raw  materials 
were  employed  after  an  initial  study  revealed  essentially  ne  impurity  affects  between 
reagent  materials  and  99.9999*  material*.  All  samples  were  formed  into  anal!  platelets 
by  briefly  remelting  on  e  graphite  plete  and  quenching  with  a  second  plate.  Samp.es 
for  a  1  obaarvation*  had  a  thickness  of  between  0.3  end  1.5  ea.  Measurements  were 
made  on  a  sample  by  sandwiching  it  between  brass  electrodes  in  a  holding  device  designed 
to  maintain  constant  contact  presau*’*.  The  samples  studied  covered  the  corpoeitional 
range  of  80:20  (80Ae^7*j:  iOAa^Sej)  to  40:60. 

Variable  temperature  studies  were  performed  b-  piecing  the  sample  holder  in  an 
thermostatically  controlled  oven.  The  sample  temperature  was  monitored  b-  e  thermocouple 
mounted  in  one  o'  the  brasa  electrode*.  Initial  measurements  were  made  a?  a  function  of 
heating  rate  in  order  to  determine  a  sufficiently  slow  rat*  to  Insure  thermal  equilib¬ 
rium.  For  temperature  increment*  on  the  order  of  10*C,  it  was  found  that  a  20  minute 
interval  between  measurement*  was  sufficient.  Studies  wer*  performea  o  cr  the  range 
of  25*C  to  140*C. 

The  I-V  characteristics  of  the  80:2G  and  the  20:30  glass  vare  obtained  by  using 
a  closed-loop  vol tage-crntrol led  constant  current  source  whose  DC  output  impedance  was 
ir.  excess  of  15  megohms.  Tea  high  source  impedance  was  necessarv  in  order  to  observe  tne 
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negative  resistance  region  of  the  materials  and  was  also  useful  in  elim  nating  the 
tendency  of  some  samples  to  switch  at  high  rates  between  states  of  different  current. 

The  stabilitv  of  the  closed-loop  ays ter  eliminated  thermal  hystersis  caused  by  heating 
of  the  output  transistors.  The  circuit  was  designed  to  provide  output  currents  o' 

0-100  r.a  over  a  dynamic  output  voltage  range  o'  0-250  V.  The  60.40,  50:50  and  40:60 
glasses  had  breakdown  voltages  too  high  to  be  measured  with  the  constant-current  source; 
their  1-V  characteristics  were  determined  by  connecting  them  to  a  manually  controlled 
0-1500  V  po-er  supply  through  a  large  series  resistance.  It  was  found  experimentally 
tsat  sweep  rates  of  less  than  10/sec  were  necessary  in  order  to  obtain  on  accurate  repre¬ 
sentation  of  the  1-V  characteristics.  A  sweep  cycle  time  on  the  order  of  1  sec.  was 
used  for  no:t  of  the  neasuremmts,  with  the  results  being  displaced  on  a  storage  oscil¬ 
loscope. 

Switching  delay  tire  measurements  were  performec  only  on  t.ie  6):2G  and  70:30 
glasses  because  of  the  inordinately  high  voltages  required  to  study  the  other  sample*, 
rhe  delay  time  was  measured  by  applying  a  voltage  pulse  across  the  ample  through  e 
100K  aeries  resistor.  .Tie  applied  voltage  and  the  vcltage  across  the  sample  were  simul¬ 
taneously  observed  or.  e  storage  oscilloscope. 

Results 

The  I-V  charecterlatics  shown  in  Fig.  1,  are  typical  of  all  the  samples  etudlee. 
.’irgin  samples  wars  initially  found  in  the  high  resistance  or  "off"  slate  as  shown  ir. 

Fig.  la.  Measured  values  of  room  temperature  conductivity  for  this  o  f  state  varied 
from  11  x  10-6  rho/cm  for  the  80:20  glass  to  0.05  x  10‘6  mho/cm  for  t.ie  40:60  ^aas. 

If  the  samples  were  subjected  to  breakdown,  the  threshold  switching  behavior  show  in 
Fig.  lb  results ,  nnd  the  sample  would  generall  cycle  in  this  manner  inde  initely. 
however,  a  sufficiently  alow  return  sweep  rate  would  lead  to  mmoty  ro*  behavior  as 
shown  in  Fig.  lc  and  14,  for  which  the  sample  remained  in  the  low  resistance  etete  ea 
current  was  reduced  to  «ero.  The  a«pe  remained  In  the  low  reeietence  etete  until 
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pulsed  by  a  high  current  (20-50  ma)  of  several  milliseconds  duration. 

Initial  efforts  to  determine  the  breakdown  voltages  of  the  80:20  glasses  were 
unsuccessful.  A  given  sample  would  exhibit  breakdown  voltages  over  a  range  of  5  to  00 
volts,  with  a  more  or  less  random  distribution  within  this  range.  To  aid  In  the  snalysls 
of  this  Inconsistency,  a  computer  program  was  written  which  would  plot  a  histogram  of 
the  number  of  breakdowns  versus  voltage  and  which  was  also  able  to  produce  the  hlstog  am 
with  various  degrees  of  statistical  smoothing.  The  mean  value  and  standard  deviation  of 
the  original  and  {smoothed  histograms  were  calculated  and  were  used  as  a  measure  of  the 
validity  of  the  smoothing  technique.  The  program  also  calculated  the  mean  and  the  stan¬ 
dard  deviation  of  the  initial  data  set.  The  standard  deviation  of  the  data  was  taken  as 
a  measure  of  the  data  taking  technique,  with  the  assumption  that  a  valid  set  of  t'ata  should 
have  a  small  standard  deviation.  The  final  program  function  was  to  produce  a  heavily 
smoothed  time  plot  of  the  data.  A  highly  erratic  time  plot  wss  taken  as  an  Indication  of 
poor  data.  Use  of  the  progrsm  in  conjunction  with  various  data  taking  techniques  pro¬ 
duced  a  steady  improvement  in  data  until  the  oolnt  was  reached  where  a  consistent  charac¬ 
terisation  of  the  device  could  be  obtained.  The  improvement  is  evident  in  the  fact  that 
the  initial  data  set  has  a  mean  of  48  volts  with  a  standard  deviation  of  36.4  volts  while 
a  data  set  taken  using  improved  techniques  yielded  a  mean  value  of  100  volts  with  *  stan¬ 
dard  deviation  of  5.6.  The  breakdown  voltage  was  found  to  be  virtually  independent  of 
the  sample  thickness. 

"ig.  2  shows  the  results  of  measurements  of  breskdown  voltage  versus  composition. 
Although  there  is  some  scatter  in  breakdown  voltages  of  different  samples  with  the  ssme 
composition,  the  compositional  trend  is  obvious.  It  was  alsu  found  that  increasing  Sc 
content  made  the  glass  a  more  stable  switch  at  the  expense  of  much  higher  breakdown 
voltages.  The  50-50  glass,  for  example,  exhibited  a  very  stable  breakdown  voltage  over 
many  breakdown  cycles  so  long  as  the  sample  was  not  subjected  to  extremely  high  on-state 
currents. 

fig.  3  shows  the  results  of  switching  delay  time  measurements  on  a  20:30  saepie. 
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Two  Interesting  results  are  apparent  froc  this  measurement.  First,  the  threshold 
switching  voltage  for  the  pulsed  case  is  nearly  a  factor  of  two  higher  than  it  is  for 
the  DC  or  steady  state  switching.  Second,  it  is  found  that  'or  values  of  applied  voltage 
slightly  greater  than  the  threshold  value,  the  switching  delav  time  is  extremely  slew, 
ranging  from  600  ms  for  the  70:30  simple  to  several  seconds  for  some  of  the  80:20  glasses. 
The  delay  time  wan  found  to  decrease  sharply  with  increasing  applied  voltage,  reaching  the 
low  millisecond  region  with  over-voltages  of  several  hundred  volts.  The  delay  time  was 
also  found  to  decrease  when  a  given  sample  was  subjected  to  repeated  breakdown  at  the 
same  applied  voltage.  Fig.  4  shows  the  results  of  the  delay-time  versus  applied  voltage 
measurements  on  the  80:20  glass.  Curve  1  la  for  the  virgin  sample,  which  exhibits  a  very 
long  initial  switching  delay  time  which  then  decreases  sharply  with  increasing  voltage. 

Above  around  550  volts,  the  rate  of  decrease  levels  off.  Curves  2  and  3  are  subsequent 
runs  taken  on  the  same  sample.  It  is  seen  that  the  behavior  for  lo"  voltages  is  quite 
different  for  the  three  curves,  but  that  at  higher  voltages,  all  the  curves  converge. 

**ig.  3  shows  the  variatlcn  in  I-V  characteristics  of  a  30-30  sample  as  a  function 
o'  temperature.  Hie  oscillogram  3(a),  taken  at  room  temperature,  indicates  a  ver*'  high 
off  state  resistance  and  no  breakdown  with  an  applied  voltage  in  excess  of  1000  volts. 

Curves  3(b),  3(c)  and  3(d),  taken  et  6S*C,  102*C  end  123*C,  respectively,  show  successive 
cecrcasas  in  the  off-stete  resistence  end  breekdovn  voltege  es  the  temperature  rises. 

In  Curve  3(c),  taken  at  140*C,  the  off-state  and  on-state  resistance  are  approaching  one 
another  in  magnitude  and  the  breakdown  voltage  haa  fallen  to  below  100  volts.  Final!**,  as 
the  temperature  is  increased  further,  the  sample  crystallite*  and  switch. ng  ctu.sc  ».  '.he 
resistance  of  the  crystalline  state,  shown  in  Curve  3(f)  is  extremely  lo**.  After  crystalli¬ 
zation,  the  sample  remains  in  the  low  resistance  state  as  the  temperature  decreases. 

rig.  6  shew*  a  plot  of  breakdown  voltage  versus  temperature  'or  tne  30-30  glass. 
if  the  field  dependent  term  in  Eq.  1  Is  neglected,  a  valid  approximation  'or  bulk  samples, 
the  equation  can  be  restated  as  follows: 
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A  nonlinear  least-square  fit  coaputer  program  has  been  written  to  fit  the  V__  versus 

BR 

T  results  to  Eq.  2.  AE  and  <hQ  were  taken  as  the  variable  parameters  of  the  fitting 
process.  A  very  good  fit  to  the  above  equation  was  obtained  for  all  compositions. 

Table  1  shows  tha  results  of  a  least  square  fit  of  data  from  a  70:30  sample.  V__(DATA) 
Is  the  measured  value  of  the  breakdown  voltage,  V(FUNCTION)  is  the  value  predicted  by 
the  fitted  equation,  and  the  error  Is  the  difference  between  the  two.  It  can  be  seen 
that  the  error  is  in  all  cases  within  the  range  of  experimental  error.  Table  II  shows 
the  values  of  AE  and  </oQ  obtained  over  the  compositional  range  of  80:20  to  50:50. 


table  i 

Least  Squares  Fit  of  the  Versus  T  Data  for  a  70-30  Glass 
with  </cQm  3.156  x  10“5  and  AF  -  0.40  ev 


Temperature 

V(Functlon) 

V(Dr.ta) 

Error 

295 

188.7 

190.0 

1.3 

308 

141.7 

137.0 

-4.7 

317 

117.9 

120.0 

2.1 

326 

99.2 

100.0 

0.8 

338 

80.1 

82.0 

1.9 

348 

67.8 

70.0 

2.2 

357 

58.8 

58.0 

-0.8 

364 

53.0 

50.0 

-3.0 

TABLE  II 

Activation  Energy  and  </<^ Ratio  for  a  Series  of 
As^Te^-As^Se^  Compositions 


Compos  1 ilon 


AE(ev.) 


*/o»(10~6) 


50-50 

60-40 

70-30 

80-2C 


0.46 

0.42 

0.41 

0.47 


84.5 

102.9 

27.2 

2.13 
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Discussion 

As2Se^  is  a  glassy  material  and  has  a  low  tendency  to  crystallize  when  cooled 
from  the  liquid  phase.  Samples  of  As2Se3  have  been  subjected  to  voltages  as  high  as 
1.5  kilivolts  without  breakdown.  Pure  As2Te3>  on  the  other  hand,  normally  does  not 
form  a  glass,  and  accordingly  has  a  high  conductivity.  The  system  x  As0Te^: (l-x)As2Se3 
can  be  prepared  in  a  glassy  state,  the  presence  of  the  selenium  being  sufficient  to 
stabilize  the  glass  for  x  j>  0.2  (15).  However,  all  evidence  indicates  that  the  As2Se^ 
is  not  actively  involved  in  the  switching  process.  Rather,  by  forming  a  glass  with  a 
high  As^Te^  content,  the  As2Se3  sets  up  the  conditions  for  the  switching  process  to 
occur. 

The  initial  breakdown  mechanism  is  unquestionably  thermal  in  nature.  This  fact 
is  subtantiated  here  by  the  independence  of  breakdown  voltage  on  thickness,  the  good 
agreement  of  the  temperature  data  with  2q.  2,  and,  in  a  qualitative  vzay,  by  the  switch¬ 
ing  delay  time  behavior,  which  agrees  well  with  that  predicted  by  Warren  (5)  on  the 
basis  of  his  solution  of  the  heat  flow  equation.  When  a  voltage  is  first  applied  to 
a  sample  in  the  off  state,  current  tends  to  become  concentrated  in  small  or  filamentary 
paths  through  the  material,  probably  because  of  local  inhomogenities.  As  the  voltage 
is  increased,  the  power  dissipation  causes  a  localized  heat  build-up.  Because  of  the 
low  thermal  conductivity  of  the  glass,  a  point  is  soon  reached  for  which  the  heat  gene¬ 
rated  is  greater  than  that  which  can  be  conducted  away  by  the  adjacent  cooler  material 
and  thermal  runaway  occurs.  This  instability  accounts  for  the  negative  resistance 
portion  of  the  I-V  characteristic.  Once  runavray  occurs,  the  temperature  of  the  con¬ 
ducting  path  increases  rapidly,  and  if  the  current  is  allowed  to  increase  sufficiently, 
phase  separation  and  crystallization  followed  by  liquif ication  in  the  localized  region 
will  occur.  As  the  current  is  reduced,  the  liquid  material  may  return  either  to  the 
glassy  or  crystalline  state,  depending  on  the  rate  of  cooling.  This  statement  is  demon¬ 
strated  by  the  DTA  data  obtained  for  the  As2Te3: As2Se3  system.  A  typical  thermogram 
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shows  an  exothermic  reaction  at  aom«  temperature  I'1  which  represents  crystallization 
of  the  material  and  an  endothermic  reaction  at  some  higher  temperature  ?2  which  is  the 
melting  point.  As  the  heating  ratio  increases,  ?1  approaches  T2,  and  for  very  rapid 
heating  rates,  such  as  occurs  in  switching,  the  exothermic-crystallisation  reaction 
may  not  occur  at  all.  However,  the  cooling  rate  in  switching  is  controlled  by  the 
external  sweep  circuit  and  can  be  made  as  slow  as  desired.  DXA  observations  indicate 
that  upon  cooling  from  the  liquid  state  the  material  will  become  glassy  if  cooled 
rapidly,  crystalline  if  cooled  slowly.  Thus,  a  rapid  cooling  or  fast  return  sweep 
rate  leads  to  threshold  type  switching,  while  a  slot;  return  sweep  rate  leads  to  memory 
behavior,  where  the  conducting  path  is  left  in  the  crystalline  conducting  state. 

The  lack  of  consistency  in  the  delav-time  measurements  of  the  80:20  and  70:30 
glasses  prevented  any  quantative  analysis  of  this  data.  The  inconsistency  appears  to 
result  from  areas  of  remanent  crystallization  brought  about  by  the  relatively  high  cur¬ 
rents  associated  with  the  pulsed  switching.  To  a  lesser  degree,  a  similar  problem  had 
jeen  observed  in  steady  state  switching  of  the  80:20  and  70:30  glasses  and  is  not  sur¬ 
prising  in  view  of  the  relatively  unstable  nature  of  these  glasses.  The  60:40  and 
50:50  glasses,  which  were  considerably  more  stable  under  DC  switching,  would  probably 
yield  useful  delay-time  data,  but  unfortunately  these  measurements  are  unfeasible  be¬ 
cause  of  the  extremely  high  voltages  required. 

Conclusions 

1.  Both  threshold  and  memory  switching  are  exhibited  by  the  As2Te3:As2Se3 
glass  system. 

2.  Switching  stability  improves  with  increasing  Se  content  at  the  expense  of 
higher  VBR. 

3.  The  switching  mechanism  in  bulk  As^As^  glasses  is  thermal  in  nature. 
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CHAPTER  I 


THEORETICAL  ASPECTS  OF  THE  MOSSIAUER  EFFECT 


In  spite  of  now  research  tools  and  recent  theoretical  advances 
there  is  sons  confualon  concernini  the  cheaical  and  Magnetic  structure 
of  the  iron  in  iron-phosphste  glasses.  Recently  Hossbauer  Effect 
Spectroscopy,  HIS,  has  been  successfully  used  to  augnent  data  obtained 
fro*  the  conventional  optical  and  Magnetic  Methods.  This  thesis  des¬ 
cribes  the  use  of  MES  to  obtain  new  structure  and  Magnetic  information 
on  iron-phosphate  glass  systans.  Although  info  mat  ion  obtained  froa 
MU  la  sin  liar  to  that  obtained  by  nuclear  quadrupole  resonance,  NQR, 
and  nuclear  Magnetic  resonance,  IMR,  the  inportance  difference  lies  in 
the  fact  that  quadrupole  nonente  of  escited  states  as  well  as  ground 
states  nay  be  observed,  revealing  valence  and  coordination  in forest ion 
concerning  isotopes  such  as  Pe57  which  has  no  ground  state  quadrupole 

The  fact  that  a  photon,  enltted  by  the  decay  of  an  excited  state, 
is  characterised  by  the  sane  energy  as  the  intoning  photon  was  denon- 
strated  by  Ubod  in  1904  in  the  resonant  scattering  of  light.  The  finite 
lifetime  of  the  excited  state  gives  rise  to  a  line-width,  T,  of  the 
enltted  radltatlon  and  is  expressed  by 


T  •  J*h 


1 


2 


The  excitation  probability  W(E)  expressed  as  a  function  of  the  energy 
of  radiation  is  given  by  the  Breit-Wigner  relation, 


i  r 

W(E)  -  j - -2 

(E-Eo)2  +  ^ 

which  is  more  generally  referred  to  as  the  Lorentzian  shape  function. 

It  will  be  noticed  that  a  resonance  occurs  only  if  the  probability  func 
tion  of  the  source  and  that  of  the  absorber  overlap  to  an  appreciable 

extent,  that  is,  if 


E0  -  {  <  1W  <  E0  +  I 

Th«  nuclear  .c.tt.rlng  proce..  1.  characterized  by  a  cross  section 
parameter,  o,  for  photons  of  energy  E  -  hv  by  the  relation 

0-0  W(E) 
o 


where 


o 

o 


21 


2d 


ex 


gr 


+  1 

TT) 


I  and  I  are  the  spins  of  the  excited  and  ground  states  respectively 
ex  gr 
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and  \q  is  the  wavelength  at  resonance.1 

The  difficulty  In  the  observation  of  gamma-ray  resonance  In 
nuclear  systems  arises  from  the  fact  that  t.  law  of  the  conservation 
of  momentum  requires  a  non-vanishing  energy  of  recoil  for  the  emitting 
system, 


E 


r 


2Mc 


where  M  Is  the  mass  of  the  system  and  c  Is  the  velocity  of  light.  In 
the  case  of  atomic  radiation,  the  emitted  quanta  of  radiation  Is  con¬ 
siderably  smaller  than  the  llnewldth  of  the  excited  state.  Zn  nuclear 

transitions,  however,  the  anergy  of  the  emitted  and  absorbed  quanta  Is 
4 

3  x  10  larger  with  a  proportional  Increase  In  the  free  recoil  energy. 

Since  the  llnewldths  of  both  phenomena  remain  about  the  same.  It  Is 

understood  that  the  energy  lost  In  recoil  Is  enough  to  bring  the  entire 

nuclear  system  out  of  resonance.  P.  B.  Moon,  In  1955,  showed  that  It 

was  possible  to  add  energy  to  the  system,  bringing  the  Incident  and 

2 

emitted  gamma  ray  back  Into  coincidence.  He  accomplished  this  by  add¬ 
ing  exactly  twice  the  velocity  of  the  recoil  process  to  the  source  by 
placing  It  on  a  rapidly  spinning  wheel,  thereby  compensating  for  the 
velocity  lost  In  both  the  source  and  the  absorber. 

In  1958,  Rudolf  Mossbauer  discovered  a  different  method  of  ob¬ 
taining  resonance  scattering  of  high  energy  photons.  He  was  lnvestlgat- 

191 

lng  the  resonance  absorption  of  the  129  Kev  transition  of  Ir.  This 
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line  has  a  low  recoil  energy  ind  is  considerably  broadened  by  the  random 
Doppler  shifting  of  velocities  due  to  thermal  agitations.  In  an  effort 
to  observe  a  decrease  in  the  amount  of  resonance  scattering,  Mossbauer 
decreased  the  temperature  to  affect  a  decrease  in  the  Doppler  shift 
variations.  Instead  of  a  decrease  in  the  resonant  scattering,  he  ob¬ 
served  an  Increase.  At  the  apparent  failure  of  the  methods  of  classical 
mechanics,  Mossbauer  finally  interpreted  his  results  in  terms  of  the 
Lamb  theory  (1939)  for  the  resonant  capture  of  thermal  neutrons  in  crys¬ 
tals.  Thus,  Instead  of  compensating  for  a  recoil  energy  loss,  Mossbauer 
had  found  a  method  by  which  the  recoil  energy  is  reduced  to  zero  and  the 
emitted  photons  are  resonantly  absorbed. 

If  the  recoil  energy,  Er,  as  calculated  for  the  free  atom  case, 
is  less  then  hv,  the  photon  energy,  and  at  the  same  time  is  less  than 
the  energy  required  to  remove  an  atom  from  its  lattice  site,  a  fraction 
of  the  processes  that  occur  will  be  recoiless  according  to. 


f  - 


Application  of  x-ray  or  neutron  scattering  yields 

.  r  4it<x2>i 

f  •  exp  [ - s — J 

r 

where  <x2>  is  the  mean  square  vibrational  amplitude  of  the  nuclear 
•notion  and  A  is  the  wave  length  of  the  gamma  ray.  It  is  evident  that 
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for  u  large  effect  to  exist  a  small  vibrational  amplitude,  such  as 
that  found  in  a  crystal  lattice,  is  desired.  For  a  Debye  solid, 


f  - 


!xp(  -  2k  eT  Cl 


♦  «  I  >  4/T  ^dx] 

e 


where  6  »  Debye  temperature,  T  -  absolute  temperature,  and  k  - 
Boltzmann's  constant.  Inspection  of  this  equation  will  show  how  the 
reduction  of  the  absolute  temperature  enabled  Mossbauer  to  make  the 
original  observation  of  the  recoilless  emission  of  gamna  radiation. 

Due  to  the  narrow  llnewldth  of  the  nuclear  transitions  involved, 
the  resonant  absorption  is  quite  sensitive  to  small  energy  variations 
in  the  gamma  radiation.  For  this  reason  it  is  possible  to  observe  di¬ 
rectly  the  interactions  between  the  nucleus  and  the  orbital  electron 
cloud.  The  manifestations  of  these  interactions  on  the  positions  of 
the  nuclear  levels  are  what  make  Mossbauer  Effect  Spectroscope  valua¬ 
ble  in  chemical  and  solid  state  applications. 

Isomer  Shift 

In  atomic  spectroscopy  one  notices  that  the  positions  of  the 
electronic  levels  are  n,-f;  fixed  due  to  the  interaction  which  exists 
between  the  levels.  There  do  exist,  analogous  to  these  interactions, 
othea  interactions  between  electronic  levels  and  nuclear  levels.  A 
change  in  the  s-electron  density  due  to  a  change  in  valence  would  di¬ 
rectly  affect  the  probability  of  finding  that  electron  in  the  vicinity 
of  the  nucleus.  As  a  result  there  is  a  Coulombic  interaction  which  may 
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shift  the  positions  of  the  nuclear  levels.  This  effect,  properly 
termed  an  "electronic  monopole  interaction"  but  more  generally  called 
the  "isomer  shift",  arises  from  a  difference  in  energy  between  an  un¬ 
perturbed  state  and  an  electrostatically  affected  version  of  that  same 
state  called  the  isomeric  state.  The  effect  of  this  monopole  interac¬ 
tion  is  to  shift  the  nuclear  levels  without  removing  their  spin  degen¬ 
eracy.  Thus,  in  the  isomer  shifting  of  these  levels,  the  center  of 
gravity  of  the  entire  Mossbauer  spectra  is  seen  to  be  displaced  in 
velocity,  usually  in  the  positive  direction.  An  expression  for  the 
expected  isomer  shift  may  be  obtained  by  first  expressing  the  electro¬ 
static  energy  of  the  nucleus,  6E,  in  terms  of  the  nuclear  radius,  R, 
given  a  specific  electronic  configuration,  |<Ko)|,  as, 

6E  -  Ze2|<JKo)|2  R2 

It  is  recalled,  however,  that  the  nuclear  transitions  take  place  be¬ 
tween  levels  so  that  the  energy  of  the  excited  gamma  ray  becomes 

«  •  SE.x-SE,i  ■  fZe2|*(0)|2(R.x2- V2) 

At  this  point  it  is  seen  that  the  absolute  value  of  the  energy  of  a 
nuclear  state  is  not  as  important  as  the  energy  of  that  state  with 
respect  to  some  arbitrary  standard  value.  If  this  is  the  case  then 
the  isomer  shift  is  given  by  the  difference  in  the  nuclear  of  two 
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Identical  materials  caused  by  a  variation  In  electron  parameters.  Thus 
see  Figure  1. 


X.S.  -  «K2  -  Rgd2><  l*,b.<°>l2  -  I* 


source 


(o)|2} 


or  by  separating  nuclear  and  atomic  parameters. 


I.S.  -  £  2.2  „2(  &  )  {  I*  „.(o)|2  -  l*.ourcJ2  > 


abs 


source' 


nuclear 


atomic 


where  6R  ■  Rex  -  Rg(j  • 

Of  course  the  above  expression  does  not  Include  any  correction 
for  relativistic  considerations  or  distortions  of  the  electronic  wave 
function  at  the  origin.  There  have  been  more  accurate  expressions 
developed  as  reported  In  summiry  form  by  D.  A.  Shirley.3  As  mentioned 
previously,  the  electronic  contribution  to  the  Isomer  shift  Is  attri¬ 
butable  directly  to  the  activity  of  the  s-electron  density  functions. 
For  an  element  such  as  iron,  chemical  bonding  involves  only  a  small 
percentage  of  the  total  s-electron  density,  therefore  one  should  expect 
that  for  any  appreciable  shift  to  occur  there  must  be  a  dependence  of 
the  s-electron  wave  functions  on  the  wave  functions  of  the  orbitals. 

As  will  be  discussed  later,  the  study  of  iron  reveals  that  the  <l-elec- 
tron  levels  have  the  greatest  influence  of  this  type.  As  a  result,  the 
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«E  .  f  Z«2  |»(0>|2  *2  (f-Utlv.^to^oint 

E  -  E. 

Isomer  Shift  ■  a  * 


IS  -  £•  z.2  l|».(0)l4-l%»>l8)  lR«2'\a1 

5  * 


P;g.  1  Isomer  shift  and  quadruple  splitting  In  Fe” 
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values  at  the  nucleus  of  the  Is,  2s ,  and  3s  restricted  Hartree-Fock 

A  n 

wave  functions  were  calculated  by  R.  E.  Watson  for  different  d  con¬ 
figurations.  The  possible  factors  which  Influence  the  Isomer  shift 
have  been  summarized  by  Shulman  and  Sugano  as  either  a  direct  contri¬ 
bution  from  As  bonding  or  an  Indirect  contribution  from  3d  bonding. 

The  Indirect  contribution  results  from  covalency  effects  between  d- 
electrons  and  either  filled  ligand  orbitals  or  empty  ligand  orbitals. 
Covalency  Involving  filled  ligand  orbitals  tends  to  Increase  d-electron 
density  thereby  reducing  s-electron  density  at  the  nucleus  by  means  of 
electrostatic  shielding,  whereas  covalency  Involving  empty  ligand  or¬ 
bitals  tends  to  reduce  d-electron  density  In  the  target  atom  thus  In¬ 
creasing  Interaction  between  s-levels  and  nuclear  fields.  As  far  as 
application  Is  concerned,  two  quantities  are  determined  from  Isomer 
shift  Investigations.  These  are  the  nuclear  radius  dimensions,  used 
In  verifications  of  nuclear  models,  and  the  electron  density  at  the 
nucleus,  of  Interest  In  molecular  ctructure  and  magnetic  studies  of  the 
solid  state. 

Quadrupole  Splitting 

So  far,  It  has  been  assumed  that  the  electron  distributions  are 
spherically  symmetric.  If  the  more  general  case  is  considered  in  which 
both  prolate  and  oblate  distortions  are  allowed,  it  1 r  >een  that  the 
(21  +  1)  degeneracy  is  the  appearance  In  a  Mossbauer  spectra  of  two  or 
more  lines  separated  by  a  difference  In  velocity  representative  of  some 
energy  AE.  This  splitting  is  cau^sd  by  an  interaction  between  the  nu¬ 
clear  quadrupole  moment,  Q,  *nd  the  gradient  of  the  electric  field  of 

f 
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of  the  electron  cloud  and  is  given  by 

H  -  Z  •  VE 


where  ■  /  pX^X^d  x 


or 


H  - 


41(21-1)  ^3I! 


id  +  i)  +  §  u+2  +  i 2 


)] 


where  n  1.  the  ...yetry  pereheter  end  I+  end  1.  ere  raising  and  lower 
ing  operators.  The  Hamiltonian  in  this  form  has  the  eigenvalues 


Eq  "  41(21-1) 


Oj  -  I,  1-1. 
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It  is  seen  that  the  second  power  of  the  magnetic  quantum  number  insures 
that  states  that  differ  only  in  the  sign  of  mx  remain  degenerate.  Since 
worV  in  iron  involves  transitions  only  between  the  ground  and  first  ex¬ 
cited  states,  it  is  seen  that  m,  -  3/2,  1/2,  -  1/2,  -  3/2  will  result 
in  a  single  doublet  for  each  distinct  value  of  the  electric  field  gradi¬ 
ent.  An  evaluation  of  the  electrical  field  gradient  tensor  is  a  rather 
complex  problem  of  solid  state  physics  that  has  been  dealt  with  to  some 

4  7»8 

extent  with  reference  to  iron  compounds. 
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In  general,  the  factors  affecting  the  electric  field  gradient 
are  charges  on  distant  ions  and  electrons  in  incomplete  shells  of  the 
target  atom  itself.  If  the  lattice  parameters  wre  known  to  good  accu¬ 
racy  it  is  possible  to  determine  the  gradient  of  the  electrical  field 
at  the  site  of  the  target  atom  but  this  may  be  quite  different  from 
the  value  of  the  same  field  at  the  target  nucleus.  An  explanation  of 
this  is  at  hand  when  it  is  recalled  thit  the  electric  field  of  some 
ligand  charge  configurations  can  lift  the  degeneracy  of  the  five  d- 
orbital  states.  Originally  the  states  appear  as  illustrated  (see  Fig¬ 
ure  2).  The  ligands  are  usually  negative  ions  or  neutral  molecules 
with  a  prominent  lone  pair  of  electrons  and  affect  the  d-electrons  of 
a  transiton  metal  atom  in  two  ways.  One  is  through  the  electrostatic 
field  of  the  negative  charges  and  the  other  is  covalent  bonding  with 
the  d-orbitals.  Those  d-orbitals  which  have  a  large  electron  density 
in  the  direction  of  the  ligand  are  repelled  and  ha ee  a  higher  energy. 
Those  orbitals  which  avoid  the  directions  of  the  ligands  are  not  pushed 
up  in  energy  so  high.  Therefore  in  the  octahedral  coordination  the  dg 
and  d  2  2  are  increased  in  energy  by  the  same  amount  leaving  d^,  dxZ 

at  a  lower  energy.  The  situation  is  exactly  opposite  in  the  tetrahed¬ 
ral  coordination  placing  d^,  dyz>  and  d^  in  a  higher  energy  state. 

The  energy  difference  in  the  two  orbital  groups  is  called  the  ligand 
field  splitting,  A,  generally  denoted  by  10Dq.  If  the  coordination  is 
distorted  from  either  the  octahedral  or  the  tetrahedral  then  the  degen¬ 
eracy  of  the  de  and  do  are  lifted  and  there  exists  an  energy,  W,  between 
each  de  or  do  suborbital.  The  population  of  these  suborbitals  is  now 
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Fig.  2  Electron  population  of  d  orbitals  for  various  electronic  configurations  of  iron  in  the 
presence  of  octahedral  and  tetrahedral  ligand  fields. 
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determined  by  the  Boltzmann  factor,  exp(-  W/kT).  It  should  be  expected 
therefore  that  a  considerable  temperature  dependence  should  be  found 
for  the  quadrupole  splitting.  The  details  of  this  temperature  depend¬ 
ence,  from  which  the  magnitude  of  the  crystal  field  splitting  may  be 
deduced,  have  been  found  for  a  number  of  materials. ^ 

The  above  discussion  has  ignored  the  relative  intensities  of  the 
lines  which  do  contain  additional  information.  It  will  be  noted  that 
the  radiation  resulting  from  the  two  possible  transitions  of  iron  are 
anisotropic  as  shown  below,1 

Angular 

Transition  Probability  of  Transition  Dependence 

1  3/2-*±l/2  1  3/2(1+cos20) 

il/2-*+l/2  1  l+3/2sin20 

where  6  is  the  angle  from  the  axis  of  symmetry.  Note  that  there  is  no 
angle  for  which  either  line  vanishes  and  the  maximum  difference  in  in¬ 
tensity  occurs  for  0-0  and  is  3:1.  This  effect  has  been  explored  in 

9  10 

some  organic  compounds  of  both  tin  and  iron 
Magnetic  Hypcrfine  Splitting 

The  most  noticeable  part  of  the  Mossbauer  spectra  results  from 
the  interaction  of  the  nuclear  magnetic  dipole  moment,  ii,  with  the  mag¬ 
netic  field  of  the  electron  cloud,  H.  The  Hamiltonian  for  this  inter¬ 
action  may  be  written, 


H  -  -  C  •  H  -  -  gu  i  •  ii 
m  n 
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where  yn  la  the  nuclear  magneton  and  g  la  the  gyromagnetic  ratio.  The 
elgenvaluea  may  be  written, 


Em  "  " 

n>x  ■  I,  1-1 . -  I 

indicating  complete  lifting  of  the  (21  +  1)  levels  with  a  splitting  be¬ 
tween  adjacent  levels  of  g^H.  For  this  reason,  the  spectra  of  metallic 
iron  shows  a  six  line  pure  nuclear  Zeeman  effect.  It  is  at  this  point 
that  the  difference  between  the  magnetic  spectra  from  the  Mossbauer 
experiment  and  those  obtained  in  other  nuclear  resonance  techniques  be¬ 
comes  apparent.  In  NMR  experiments,  for  example,  nuclear  transitions 
take  place  between  the  split  levels  of  a  nuclear  state  (An^  ■  ±  1) 
whereas  in  the  Mossbauer  experiment  transitions  occur  between  nuclear 
levels  themselves  with  the  condition  that  the  z-component  of  the  angular 
momentum  transferred  to  the  emitted  quanta  of  gamma  radiation,  L,  must 
be  conserved.  The  final  selection  rules  may  then  be  written, 

I1!  ~  I2I  -  L  -  IX1  +  X2I  L  *  0 

Anij  -  0,  ±  1 

The  appearance  of  hyperfine  splitting  makes  possible  the  direct 
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computation  of  cha  transition  probabilities.  Thaaa  ara  obtained  from 
tba  squares  of  tha  Cleih-Corden  coafflclents  which  ara  written 


It  la  also  banaflelal  to  eoaaldar  tha  angular  dapandanca  of  thaaa  tran¬ 
sitions.  Thaaa  ara  also  glvan  in  tarma  of  tha  angla  between  tha  mag¬ 
netic  field  and  tha  diractlon  of  propagation,  0.  It  is  obaarvad  from 
thaaa  expressions  that  tha  tranaltion  for  Am  •  0  vanishes  at  0  •  0  but 
that  tha  avaraga  over  a  aphere  for  each  corns*  aent  angular  dapandanca 
gives  a  relative  probability  ratio  of  3s2sltl:2t3.  Th*  mm  of  the  angu¬ 
lar  aspraaaiona  la  apharically  aynmatrlc.1 

It  haa  been  shown  that  in  eone  natariala  tha  magnetic  field  arla- 
lng  from  tha  lattice  itself  is  sufficient  to  causa  magnetic  hyper fine 
splitting.  Investigation  of  iron  foil  has  yielded  an  internal  field  of 
3.33  x  10*  KOe  at  room  temperature.11  Th.s  value  was  derived  from  the 
fact  that  the  Internal  field  is  proportional  to  the  total  splitting  in 

the  apectrisi  (the  velocity  difference  between  the  two  most  distant 

12 

lines).  Tha  origin  of  this  internal  field  is  primarily  due  to 

a)  The  Feral  contact  interaction  (direct  coupling  between  the 
nucleus  and  unpaired  electrons). 

Ha  -  -  t5-  B  <  £(♦  *#2(o)  -  ♦  *b2(o))> 


i:> 


! 
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b)  The  orbital  magnetic  moment  contribution 


»L 
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c)  Dipolar  interaction  with  electron  spin 


«D 


r 
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The  appearance  of  these  terms  is  based  on  the  electron  configuration. 

A  complete  treatment  is  given  by  Freeman  and  Watson. 

It  ia  necessary  to  point  out  that  the  three  primary  features  of 
the  Mossbauer  spectra  isomer  shift,  quadrupole  splitting,  and  magnetic 
hyperfine  splitting  may  occur  simultaneously  without  regard  to  the  mech- 
anlsms  of  their  origin,  see  Figure  3.  In  order  to  split  a  spectra  into 
its  separate  components,  some  a  priori  knowledge  of  the  system  para¬ 
meters  in  conjunction  with  computer  fitting  techniques  must  be  used. 

An  interactive  analysis  program  of  this  type  was  used  by  this  author 
and  ia  described  in  Appendix  II.  In  order  to  successfully  predict  a 
spectra  resulting  from  a  number  of  interrelated  parameters,  it  is  nec¬ 
essary  to  construct  the  Hamiltonian  matrices  for  both  the  ground  and  ex' 
cited  states  of  the  nucleus  and  then  calculate  the  transition  probabil¬ 
ities  from  the  resulting  eigenvalues  and  eigenvectors. 


r 
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CHAPTER  II 

mOssbauer  effect  studies  of  amorphous  materials 

Before  Mossbauer's  discovery  in  1958,  the  possibility  of  observ¬ 
ing  hyperflne  structure  of  the  nucleus  hsd  been  all  but  discounted. 
However,  the  fact  that  gamma  ray  absorption  of  a  material  with  thermal 
vibrations  of  the  order  of  10^  cm/sec  could  be  affected  by  the  addition 
to  the  system  of  a  velocity  of  10~J  cm/sec  gave  new  life  to  areas  such 
as  magnetism,  lattice  dynamics,  chemical  bonding,  and  nuclear  physics. 
Although  most  of  the  early  theoretical  work  approached  the  problem  of  a 
derivation  of  the  expression  of  the  absorption  cross  section  area  by 

assuming  the  nucleus  was  bound  In  a  crystal  lattice  such  as  the  work 
14 

done  by  Maradudln,  several  experiments  were  undertaken  to  investigate 
the  possibility  of  a  Mossbauer  effect  In  glassy  materials.  Studies  of 
this  nature  were  reported  by  Gol'danskll^  In  1963,  who  successfully 
used  MES  to  study  the  chemical  state  of  nuclei  bound  in  both  the  crys¬ 
talline  and  the  vitreous  states. 

It  Is  recalled  from  Chapter  1  that  the  recoil-free  fraction  In 
the  resonant  process  may  be  written 


r  4tt2  <x2>  o 
f  -  exp  [ - 5 -  )  J 


where  <xZ>  Is  the  mean  square  amplitude  of  the  vibration  in  the  direction 
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of  the  emission  of  the  gaasu  rey  averaged  over  the  lifetime  of  the 
nuclear  state  involved.  It  la  seen  that  if  <x  >  is  large,  that  is  to 
say,  if  the  nucleus  is  unbounded,  the  expression  for  the  fraction  would 
vanish  and  no  effect  would  be  observed.  It  would  seem  clear  then  that 
no  Mossbauer  effect  could  occur  in  a  liquid  as  was  demonstrated  by  E. 
Segre  (1933). 18  Although  no  long  range  order  occurs  in  a  glass,  the 
short  range  order  is  sufficient  to  bind  a  nucleus  to  suitably  small 
dimensions  and  has  been  shown  to  give  a  substantial  effect.  In  addi¬ 
tion  to  the  work  of  Gol'danskll,  early  demonstrations  of  this  Include 
57Fe  in  fused  quartz  and  silicate  glass,17  iron  pentacarbonyl  frozen 
in  an  organic  3olvent,18  and  119 Sn  in  polymethyl  methacrylate.19 

Although  iron  is  not  genrally  a  major  constituent  of  inorganic 
glass,  except  as  a  bothersome  Impurity,  its  structural  information  may 
be  very  beneficial  in  deducing  the  overall  structure  of  the  compound. 
Because  it  is  a  transition  metal,  other  magnetic  and  optical  spectro¬ 
scopic  techniques  in  addition  to  Mossbauer  spectroscopy  make  structural 
conclusions  rather  complete.  For  this  reason,  most  of  the  recent  stu¬ 
dies  favor  iron  as  the  doping  element. 22*23*2^*27 

As  a  result  of  investigations  using  several  experimental  methods 

27 

it  has  been  generally  concluded  that  iron  in  glass,  in  addition  to 
the  divalent  and  trlvlanet  states,  may  be  present  in  various  other  coor¬ 
dination  states: 

a)  It  may  replace  the  groups 


[  03/2  ]  OM 
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to  give  the  structures 


[  F.111  04/2  ]  ‘  M+ 


or 


[  Fe11  04/2  ]2“  2M+ 

where  the  Iron  Is  the  vitrlfler  with  the  coordination  number 
U  (M  ■  LI,  Ns,  K); 

b)  It  may  appear  aa  a  modifier  In  the  coapoaltlon  of  the 
groups 


[  P.111  06  ] 
and 

[  Fe11  06  ] 

with  a  coordination  number  6,  or; 
c)  It  may  form  oxides  or  alkaline  ferrites  In  the  collold- 
dlsperse  states. 

If  we  assume  that  In  the  preparation  of  the  glass,  there  exists 

some  equilibrium  distribution  of  the  different  forms  of  the  Ion,  then  a 

27 

chemical  expression  may  be  written  of  the  form 


bO 
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C  F.111  04/2  ]  M  re*  Fe*  [  Fe11  0,/2  ]  2M 
(group  2)  (group  3) 

[  Fe111  0  ]  [  Fe11  0  ] 

(group  4)  (group  5) 

# 

3+ 

where  groups  4  end  5  cen  be  considered  to  consist  of  simply  Fe  and 
Fe2+  Ions  respectively.  The  fact  that  the  divalent  and  trlvalent  elec¬ 
tron  configurations  result  In  substatually  different  electric  and  mag¬ 
netic  environments  for  the  target  nucleus  should  Insure  that  the  spec¬ 
tra  resulting  from  these  two  states  will  be  readily  distinguishable. 

As  mentioned  previously,  the  Isomer  shift  arises  In  the  fact 
that  the  electronic  configuration  in  the  source  and  absorber  are  not 
identical.  This  shift  is  principally  attributable  to  a  perturbation  of 
the  wave  functions  of  the  3s  orbitals  in  Iron  so  that  the  difference  In 
Isomer  shift  between  spectra  of  Fe  and  Fe  ,  which  differ  only  by  one 

3d  electron,  is  not  Immediately  understood.  Since  a  single  3d  electron 

2 

should  not  appreciably  affect  the  electronic  charge  density,  |4#(o)  |  , 
the  shift  must  originate  in  the  fact  that  the  3s  electrons  are  spending 
some  of  their  time  farther  from  the  nucleus  than  the  3d  electrons.  As 
a  result,  the  electrostatic  influence  of  the  3s  electrons  are  directly 
affected  by  the  screening  effects  of  the  lower  3d  levels.  The  addition 
of  a  single  3d  electron  reduces  the  Coulombic  attraction  between  the  3s 

:  6i 
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level  and  the  nucleus  allowing  the  3s  level  to  expands  The  result  is  a 

substantial  shift.  Belyustin  et  al.  found  from  a  compilation  of  isomer 

shift  data  in  five  systems  of  glasses  that  in  variable  concentrations 

of  F*2°3  that  the  shlft  for  each  lon  remained  more  or  less  c instant  in 

concentrations  of  more  than  4  mole  X  at  0.037  i  0.003  cm/sec.  Along 

with  other  da.'a  indicating  that  the  llnewldths  also  remained  constant, 

they  concluded  that  a  "short-range"  order  exist3  for  itor.  atoms  in  both 

20 

valence  states.  Work  done  on  alkali  sll'cate  glasses  yielded  values 

*V» 

of  0.008  to  0.009  cm/ sec  for  high  spin  Fe  and  0.083  cm/sec  for  hig;i 
spin  Fe2+.  Results  of  experiments  done  on  silicate  and  phosphate 

24 

glasses  are  shown  in  Figure  4  along  with  certain  similar  crystals. 

3+ 

While  Fe(P03)3,  Fe(P04)  •  2H20,  and  NaFeSi^  are  known  to  contain  Fe 
in  octahedral  coordination,  crystalline  Fe(P04)  and  KFeSijOg  are  known 
to  have  the  trivalent  iron  in  tetrahedral  coordination.  The  grouping 
in  Figure  4  is  evident  and  the  results  are  generalized  to  show  that 
while  silicate  glasses  tend  to  have  Fe  in  4-fold  coordinated,  the  tri¬ 
valent  ion  tends  toward  6-fold  coordination  in  phosphate  glasses.  This 
correlation  between  isomer  shift  and  coordination  number  was  shown  by 
Walker  et  al.33  by  comparing  the  behavior  of  the  ion  in  both  coordina¬ 
tion  sites.  Increased  covalency  in  3d3  and  3d**  compounds  is  general. y 

34 

the  result  of  increases  in  4s  electron  contribution.  Since  isomer 

shifts  are  directly  dependent  on  s-electron  densities,  it  is  expected 

3+ 

that  a  change  of  coordination  number  from  6  to  4  for  either  the  Fe  or 
pe2+  result  in  a  decrease  in  the  isomer  shift.  Kurkjian  and  E.  A. 

Sigaty2*  found  isomer  shifts  for  tetrahedral  sites  of  the  order  of  0.001 
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Figure  4.  Comparison  of  Isomer  Shifts  in  Typical  Silicate  and  Phosphate  Glasses  (After  Buchanan 

and  Siegety)24 
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ca/sec  while  those  for  octahedral  sites  were  substantially  larger,  of 
the  order  of  0.020  cm/sec.  Although  verification  of  this  correlation 
for  the  divalent  is  less  conclusive  due  to  the  small  amount  of  work, 
some  results  are  listed  in  Tables  1  and  2 . 

Since  it  is  widely  accepted  that  the  structure  of  glasses  is 
vastly  distorted  from  cubic  symmetry,  one  would  expect  that  there  would 
be  some  value  for  the  electric  field  gradient  at  the  nucleus  from  the 
d-shell  electrons  of  the  target  atom  as  well  as  the  electrostatic  field 
of  the  ligand  structure.  It  is,  therefore,  no  surprise  that  in  most 
glasses  which  contain  iron  quadrupole  splitting  is  found  for  both  ion 
states.  The  trivalent  ion  splitting  in  silicate  glasses  has  been  to 
have  a  range  of  0.07-0.10  cm/ sec  which  can  therefore  be  considered  in¬ 
dicative  of  tetrahedral  coordination.  The  phosphate  glasses  show  a 
somewhat  wider  range  for  the  octahedral  coordination  of  0.032-0.088  cm/ 
8ec.^  Splitting  less  than  approximately  0.03  cm/sec  or  greater  than 
0.10  cm/sec  would  seem  to  indicate  some  combination  of  the  two  states. 

As  with  the  isomer  shift  of  the  divalent  ion,  the  quadrupole  splitting 
has  been  incompletely  explored,  but  reported  values  range  from  0.200- 
0.229  cm/ sec  at  300  degrees  Kelvin.  There  have  been  attempts  to  quan¬ 
titatively  analyze  the  substantial  temperature  dependence  of  the  split¬ 
ting  values  with  respect  to  distorted  octahedral  symmetry  by  fitting  the 

data  to  the  values  of  the  energy  differences  between  the  three-fold  de- 

6  21  22 

generate  de  levels  as  published  by  Ingalls.  ’  ’ 

25 

Work  done  by  A.  Lerman,  M.  Stiller,  and  E.  Hernon,  showed  that 
the  ratios  of  the  areas  under  each  of  the  two  superimposed  doublet  spec¬ 
tra  varied  linearly  with  the  relative  proportion  of  each  respective  ion 

€4 
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TABLE  1 

FE2+  PARAMETERS  IN  SILICATE  MATERIALS 


(I*S*)Cu 

(cm/sec)  C.N. 


0.053 

0.060-0.077 

0.081-0.101 

0.110 


4 (square  planar) 

4 (tetrahedral) 

6 (octahedral) 

8 (distorted  cube) 


TABLE  2 

FE2+  ISOMER  SHIFT  DATA  FOR  GLASSES  AND  REFERENCE 


Glass  Type 

(I-S'>Cu 
(cm/ sec) 

Ref. 

Phosphate 

0.115 

46 

Phosphate 

0.107 

21 

Phosphate 

0.095-0.097 

47 

Silicate 

0.089 

48 

Silicate 

0.080 

46 

Silicate 

0.067-0.089 

27 

26 


in  Che  material.  Results  of  their  experiments  with  mixtures  of  vivia- 

nite  and  amorphous  ferric  phosphate  within  3%  of  conventional  chemical 

35 

methods  were  reported.  At  about  the  same  time,  Karyagln  showed  that 
the  intensities  of  the  two  peaks  in  the  doublet  should  not  be  equal  in 
the  general  case  due  to  the  anisotropy  of  the  recoil  free  fraction  re¬ 
gardless  of  the  isotropy  of  the  crystal.  This  "Gol'danskii-Karyagin 

Effect"  may  be  employed  to  deduce  the  sign  of  V22  in  both  ^Fe  an(j 
1 1 9m  'iA 

Sn  spectra.  This  fact  would  make  the  separation  of  a  composite 
spectra  quite  difficult.  However,  it  has  been  determined  that  a  mechan¬ 
ical  mixture  of  two  glasses,  each  of  which  contain  only  one  ion  and  does 

give  a  symmetric  doublet,  yields  an  unresolved  three  line  spectra  that 

,  27 

is  quite  similar  to  that  obtained  in  most  composite  glasses.  There¬ 
fore,  the  assumptions  that  must  be  made  for  a  glass  under  investigation 

to  allow  this  type  of  analysis  are  that  the  absorber  is  thin  enough  to 

28 

allow  pure  Lorentzian  line  shape,  both  iron  valence  states  undergo 
quadrupole  splitting,  and  that  there  is  no  anisotropy  in  the  reemission 
of  gamma  radiation. 

One  of  the  outstanding  capabilities  of  the  Mossbauer  Effect  is 

the  opportunity  for  evaluating  the  magnetic  behavior  of  the  iron  in 

glassy  systems.  Hyperfine  splitting,  hfs,  has  been  observed  in  ferro- 

and  antiferromagnetic  materials  as  a  result  of  the  internal  magnetic 
29 

field.  In  di-  and  paramagnetic  substances,  the  ground  state  and  ex- 

30 

cited  state  degeneracy  is  removed  by  an  external  magnetic  field,  how¬ 
ever  it  should  be  possible  to  observe  hyperfine  splitting  without  an 
external  field  provided  that  the  electron  spin  relaxation  time  is  large 

£6 


i 


27 


compared  to  the  reciprocal  of  the  hyperfine  interaction  energy,  that 

h  26 

is,  if  T  >  —  where  A  is  the  magnetic  hfs  energy.  In  investigations 

■A 

of  dilute  concentrations  of  iron  in  sodium  trisilicate  glass,  Kurkjian 
and  Buchanan  have  shown  that  as  the  temperature  is  lowered  a  fairly 
distinct  six  line  spectra  is  evident  with  a  large  doublet  completely 
obscurring  the  center  two  peaks.  This  could  not  be  interpreted  on  the 
basis  of  a  precipitation  of  an  iron  oxide  since  the  hfs  spectra  de¬ 
creased  with  an  increase  in  iron  concentration.  The  effect  could  only 
be  explained  in  terms  of  a  long  spin  relaxation  time  as  discussed  above 

since  the  glass  was  known  to  be  paramagnetic  at  the  low  temperature. 

3+ 

As  the  concentration  is  increased,  the  proximity  of  neighboring  Fe 

causes  spin-spin  coupling  relaxation. 

The  "inverse  spinel"  effect  was  shown  in  work  by  Belyustin  et 

al.27  with  various  concentrations  of  Fe^  in  sodium  silicate  glass. 

Two  six  line  spectra  overlapping  on  the  last  three  lines  with  an  in- 

3+ 

tense  trivalent  doublet  indicated  eight  Fe  on  tetragonal  A  sites  and 

iL  24- 

eight  Fe  and  eight  Fe  ions  on  octahedral  B  sites  at  room  tempera- 

11 

ture.  Shaw  and  Heasly  have  investigated  heat  treatment  on  Na20-Fe203- 
Si02  and  Li20-Hn02-Si02  glasses  and  have  interpreted  their  results  on 
the  basis  of  superparamagnetism. 

Finally,  it  is  understood  that  a  combination  of  magnetic  hfs  and 

26 

quadrupole  coupling  may  exist  simultaneously.  Kurkjian  and  Buchanan 
investigated  the  spectra  of  Fe(P03>3  in  both  the  vitreous  and  the  crys¬ 
talline  states.  While  the  crystal  showed  no  quadrupole  splitting,  the 
value  of  the  splitting  for  the  divalent  ion  in  the  glass  (20%  of  the 


67 


28 


total  Iron  present)  Indicated  octahedral  coordination.  Below  10  degrees 
the  glass  showed  the  same  quadrupnle  spectra  superimposed  on  a  six- 
line  hyper fine  structure.  Although  the  hfs  presumably  resulted  from  a 
large  fraction  of  the  iron  atoms  in  an  internal  magnetic  environment, 
the  fact  that  the  proportions  of  the  amplitudes  were  not  3:2:1:1:2:3  as 
found  in  metallic  iron,  seemed  to  indicate  the  possibility  of  other  spin 
states  of  the  trlvalent  ion  that  were  not  Included  in  the  hfs.  It  was 
postulated  that  although  the  glass  was  known  to  be  paramagnetic,  either 
some  short  range  magnetic  ordering  is  present  or  the  degeneracies  are 
lifted  by  a  long  relaxation  time. 

Finally,  we  notice  that,  Just  as  in  most  other  methods  of  spec¬ 
troscopy,  the  llnewidths  are  substantially  Increased  in  glassy  mater¬ 
ials.  This  broadening  is  usually  explained  as  due  to  the  variation  in 
site  parameters.  For  paramagnetic  glasses  this  broadening  is  most 

32 

likely  due  to  Incomplete  relaxation  of  the  magnetic  hyperflne  fields. 

27 

Belyustln  et  al.  concluded  that  the  line  broadening  which  occurs  with 

a  reduction  in  iron  concentration  may  be  explained  by  the  formation  of 
II  2- 

Fe  °4/2  (8rouP  3)  in  which  the  iron  acts  as  a  vitrifier.  Although 
u 

Fe  usually  plays  the  role  of  a  vitrifier  due  to  its  smaller  atomic  ra¬ 
dius,  the  weak  concentration  may  allow  the  divalent  ion  to  appear  in 
the  network.  This  suggestion  has  also  been  made  by  other  authors,  but 
the  data  is  inconclusive  since  small  concentrations  give  low  effects  and 
statistical  errors  are  great.  Due  to  a  half-life  of  the  excited  nuclear 
state  of  about  100  ns.,  the  theoretical  minimum  value  of  half-width  is 
approximately  0.010  cm/sec  in  iron.  The  minimum  observable  half-width 


i 

i  . 
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29 


Is  simply  twice  this  due  to  the  additive  combination  of  source  and 
absorber  line  wldtha.  The  Inatriaental  Unewldth  obtained  by 
Belyustln  et  al.  was  0.040  cm/sec  and  the  trlvalent  doublet  waa  broad¬ 
ened  by  another  0.030  ca/sec.  This  may  be  due  to  the  presence  of  the 

Ion  In  both  tetrahedral  and  octahedral  coordination  statea.  Since  the 

3+ 

quadrupole  splitting  for  Fe  Is  of  the  order  of  the  llnewldth,  the 
superposition  of  the  two  coordinations  would  be  very  difficult  to  re¬ 
solve. 
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CHAPTER  III 

EXPERIMENTAL  TECHNIQUES 

MdTssbauer  spectra  are  obtained  by  irradiating  the  nucleus  of 
certain  isotopes  with  a  narrow  band  of  monoenergetic  gamma  radia¬ 
tion  and  tabulating  the  intensity  of  the  emitted  radiation  with  re¬ 
spect  to  some  geometry  of  the  system.  Although  detection  of  the  gamma 
radiation  is  well  accomplished  through  existing  counting  techniques, 
the  problem  of  obtaining  a  swept  range  of  ganma  ray  frequencies  is  a 
novel  engineering  problem.  A  workable  solution  is  to  use  a  radioac¬ 
tive  source  msterlal  which  decays  lso-energetlcally  with  the  osotope 
contained  in  the  absorber.  A  band  of  energies  is  the  obtained  by 
Doppler  shifting  the  affective  frequency  of  radiation  at  the  absorber 
by  adding  a  velocity  to  either  the  source  or  the  absorber.  An  ex¬ 
pression  for  this  change  in  energy  is  given  by  E  »■  (v/c)E,  where  E  is 
the  unshifted  energy  and  v  is  the  added  velocity  defined  as  positive 
when  moving  toward  the  nucleus.  The  resultant  spectra  is  then  given 
as  absorption  versus  velocity. 

Basically  there  are  two  methods  which  have  beon  employed  in 
Doppler  shifting  the  source  radiation.  The  simplest  of  these  is  the 
constant  velocity  method  which  Mossbauer  used  in  his  original  experi¬ 
ments.  In  this  method,  a  displacement  is  added  to  either  the  source 
or  the  absorber  by  means  of  mechanical  cams,  moving  screws  and  gears, 


30 


31 


rotating  wheels  piezoelectric  crystals,  or  loudspeakers  executing  a 

39 

constant  velocity  motion.  During  a  predetermined  time  interval  the 
detected  radiation  is  counted  and  recorded.  The  velocity  is  then 
altered  by  a  small  amount  and  the  process  is  repeated.  This  is  contin¬ 
ued  until  an  entire  spectra  has  been  accumulated.  An  advantage  of  this 
method  is  that  any  part  of  a  spectra  may  be  examined  in  detail  simply 
by  reducing  the  increments  of  velocity  and  the  limits  to  which  it  is 
altered. 

The  other  method,  the  one  used  by  the  present  author,  applies 
motion  to  either  the  source  or  the  absorber  in  such  a  way  that  all  the 
velocities  of  interest  are  covered  in  one  cycle.  This  motion  is  usu¬ 
ally  applied  by  means  of  a  loudspeaker  coil  upon  which  has  been  added  a 
pick-up  coil  that  Indicates  the  actual  movement  being  experienced.  The 
coil  is  modulated  with  a  triangular  wave  form  and  the  pick-up  coil  gen¬ 
erates  an  error  signal  which  is  used  by  the  drive  unit  in  a  feed-back 
37  38 

control  loop.  ’  A  multichannel  analyzer  having  the  same  sweep  rate 

as  the  drive  unit  samples  each  successive  cycle  in  increments  of  time 

equal  to  the  reciprocal  of  the  product  of  the  sweep  frequency  and  the 

number  of  channels.  The  number  of  counts  in  each  Increment  are  the 

stored  separately.  The  non-availability  of  a  multichannel  analyzer  can 

be  overcome  by  reducing  all  detected  pulses  to  the  same  height  and  then 

modulating  these  pulses  with  the  velocity  drive  signal.  The  result  is 

then  accumulated  in  a  pulse  height  analyzer.  This  method  provides  an 

40 

automatic  cancellation  of  an  non-linearities  in  the  drive  signal. 

The  apparatus  used  in  the  current  investigations  employs  the 


addressing  signal  for  the  last  bistable  address  register  of  a  TMC 
model  404  400  channel  multichannel  analyzer  to  drive  a  Nuclear  Science 
and  Data  Corporation  model  AM-1  velocity  drive  unit  (see  Figure  5) . 

The  analyzer  signal  Is  a  square  wave  with  a  period  of  100  microseconds 
corresponding  to  the  period  of  sweep  of  the  analyzer.  The  AM-1  drive 
unit  Integrates  this  square  wave  to  produce  the  triangular  signal  used 
for  the  modulation  of  the  drive  coll.  The  feed-back  control  loop  cor¬ 
rection  coll  Is  also  utilized.  The  drive  signal  Is  passed  through  a 
D.C.  power  amplifier  stage  before  driving  the  armature  of  the  trans¬ 
ducer  and  a  velocity  control  potentiometer  sets  the  peak  value  of  this 
driving  waveform  establishing  the  range  of  velocities  to  be  scanned. 

The  linearity  specifications  cover  only  95%  of  the  waveform,  introduc¬ 
ing  a  5%  nonlinearity  at  the  apex.  Two  bands  of  velocity  are  provided 
0-10  cm/sec  and  0-60  cm/sec,  of  which  the  former  was  utilized  exclu¬ 
sively. 

Detection  of  the  gamma  radiation  was  accomplished  by  means  of  a 
Reuter-Stokes  303PC331  proportional  counter.  The  tube  uses  Krypton  gas 
at  a  pressure  of  72  cm/Hg  and  a  Berylium  window  1"  in  diameter  with  a 
thickness  of  0.010".  The  applied  voltage  was  1900  volts.  Output  pulses 
were  amplified  by  a  Tennelec  model  TC133  FET  preamplifier  and  a  Tennelec 
model  TC202  BLR  linear  amplifier  with  an  Intervening  Tennelec  model 
TC441  single  channel  analyzer.  The  single  channel  Is  required  to  Iso¬ 
late  the  desired  line  of  the  source  spectra. 

There  are  at  least  32  different  element*,  in  which  the  Mbssbauer 
effect  has  been  found  to  exist  in  at  least  one  of  several  isotopes.  The 
objective  In  any  study  is  to  find  a  single  line  decay  of  an  unstable 
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parent  isotope  which  is  iso-energetic  with  the  nuclear  transition  to 
be  observed.  It  is  then  generally  desirable  to  diffuse  the  source 
material  into  a  cubic  host  matrix  in  order  to  eliminate  quadrupole 
coupling  and  magnetic  splitting.  The  source  used  for  the  present  work 
was  Co^  diffused  into  the  cubic  environment  of  copper  foil  manufac¬ 
tured  by  New  England  Nuclear  Corporation.  After  electroplating,  the 
source  was  annealed  for  three  hours  at  1000  degrees  centigrade  in  hy¬ 
drogen  and  quenched.  The  width  of  the  14.4  Kev  resonance  line  in  a 
K^Fe(CN)6  •  3H20  single  line  reference  absorber  with  a  thickness  of 
0.5  mg/cm  was  0.036  cm/sec.  In  order  to  reduce  the  Intensity  of  the 
6.5  Kev  x-ray  present  in  the  source  spectra,  a  3/8"  thickness  of  luc- 
lte  was  inserted  between  the  source  and  the  absorber. 

The  calibration  of  the  apparatus  was  accomplished  by  comparing 

the  positions  of  the  six  spectral  lines  obtained  from  0.001"  iron  foil, 

see  Figure  6,  with  the  reference  splittings  published  in  the  lltera- 
41 

ture.  The  velocity  difference  between  channels  was  assumed  constant 
although  examination  of  Figure  7  indicates  that  some  non-linearity  did 
exist.  All  measurements  of  Isomer  shift  were  made  with  reference  to 
the  single  absorption  line  of  stainless  steel. 

The  accumulated  data  was  extracted  from  the  multichannel  analy¬ 
zer  by  means  of  an  IBM  typewriter  which  printed  the  total  number  of 
counts  in  each  channel.  The  count  totals  were  then  transferred  to 
punched  cards  at  the  expense  of  considerable  time  and  effort.  As  one 
can  imagine,  the  novelty  and  practicality  of  such  a  method  diminishes 
very  rapidly  and  an  alternative  would  have  been  desirable. 


Kir,.  6  Calibration  curve  for  .001"  iron  foil 
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Except  for  the  iron  and  stainless  standard  foils,  all  samples 
were  finely  powdered.  A  mixture  of  400  mg  of  sample  powder  added  to 
2  grams  of  powdered  urea  was  then  ulaced  between  two  solid  brass  cylin¬ 
ders  and  pressed  in  a  hand  vise.  The  result  was  an  opaque  disc  1-1/2" 
in  diameter  and  1/16  inch  thick.  This  procedure  produced  an  iron  den¬ 
sity  of  approximately  13  mg/cm2.  In  order  to  determine  the  amount  of 
broadening  due  to  this  heavy  iron  concentration,  the  density  was  re¬ 
duced  to  about  0.25  mg/cm2  with  only  a  10X  reduction  in  linewidth. 
Resolution  at  this  concentration,  however,  was  very  poor. 

After  the  discs  had  sintered,  they  were  mounted  in  an  aluminum 
holder  which  was  then  mounted  to  the  bottom  work  surface  of  a  Texas 
Instrument  CLP  1  litre  rryoflask  with  brass  screws.  In  order  to  pro¬ 
vide  good  thermal  conduction  from  the  diac,  one  face  was  coated  with 
vacuum  grease  and  pressed  into  contact  with  a  thickness  of  aluminimi 
foil  which  was  then  placed  into  contact  with  the  aluminum  mounting 
bracket.  The  only  difficulty  .neutered  with  this  arrangement  was 
that  when  the  system  was  brought  to  room  temperature  from  77  degrees  K 
there  was  a  tendency  for  the  melting  frost  to  dissolve  the  urea  disc. 
This  was  remedied  by  allowing  the  evacuation  pump  to  remain  in  opera- 
t ion  during  the  warm-up  period. 

The  glasses  were  prepared  using  a  physical  mixture  of  55  mole  X 

FeO  and  45  mole  X  ?20y  The  mixture  WiS  melted  in  air  at  1300  degrees 

C  for  one  hour  and  poured  into  but tor s  one  cm  thick  and  two  cm  in  dia- 

2+ 

meter.  To  vary  the  different  relative  concentrations  Fe  and  Fe  in 
the  glass,  dextrose  (sample  preparations  I,  II,  and  III,  respectively) 
,dd,d  to  the  Mixture  to  red*.  F.*  to  F.2+  -*,11.  In  th.  -It.  *11 
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CHAPTER  IV 

RESULTS 


General 

The  Mossbauer  If  feet  Spectra  of  the  55-45  mole  X  FeO-P^u^ 
glasses  at  roon  temperature  consist  of  three  partially  resolved  lines 
of  unaqual  magnitudes.  These  Unas  are  assuaed  to  be  the  superposi¬ 
tion  of  two  independent  doublets  which  result  from  the  presence  of  tha 
two  different  ions.  When  tha  samplea  are  cooled  to  78°K  the  observed 
effect  la  substantially  dacreased.  This  la  accompanied  by  an  outward 
migration  of  the  line  positions  and  an  Increase  in  their  linewldths. 
Heat  treatment  of  tha  50  l  Fe^/Fa  sample  yielded  additional  lines 
considerably  dlaplaead  from  sero  velocity. 

Room  Temperature  Data 

Room  temperature  spectra  for  the  three  Iron  ratios  were  ob¬ 
tained  and  ara  shown  as  Figures  8,  9,  and  10.  Since  tha  daxtrosa 

which  was  added  to  the  melt  in  preparation  of  the  samples  served  only 

3+ 

to  change  tha  oxidation  state  of  a  portion  of  the  Fe  ion  concentra- 
2+ 

tlon  to  Fa  ,  one  wruld  expect  thr.t  tha  overall  iron  concentration  and 
thus  the  amount  of  total  affect  observed  In  all  three  samples  remains 
tha  same.  That  this  Is  indeed  the  case  is  demonstrated  in  the  fee* 
that  the  Integral  under  the  entire  spectrum  remains  constant  regard¬ 
less  of  the  relative  concentration.  By  employing  a  least-squares  ap¬ 
proximation  procedure,  two  symmetric  doublets  were  fit  to  tha  raw  data, 


38 


FOLDED  SPECTRA  OF 
15Z  F(» 3*  /F*»  SAMPLE  AT  300 


VELOCITY  tCM/SEC) 


FOLDED  SPECTRA  OF 


-1.60  -1.20  -0.80  -0.40  0.00  0.40  0.80 

VELOCITY  (CM/SEC) 


FOLDED  SPECTRA  OF 
3*  /Fft  SAMPLE  AT  300 


VELOCITY  fCK/SECI 


42 


and  the  r  suiting  line  positions  and  relative  amplitudes  are  shown 
under  each  plot#  The  solid  line  passing  through  each  sec  of  experimen¬ 
tal  data  points  shows  the  total  fitted  spectra.  The  actual  numerical 
results  for  each  fit  calculation  is  given  in  tabular  form  with  the  same 
title  as  the  plot  in  Appendix  II.  In  the  case  of  each  doublet,  the 
amplitudes  of  the  two  corresponding  peaks  were  varied  simultaneously 

as  one  variable,  whereas  each  position  was  varied  independently.  In 

3+ 

the  case  of  the  15  %  Fe  / Fe  sample,  it  was  found  that  the  best  fit 
occurred  for  the  case  in  which  the  negative  peaks  of  both  doublets  were 
made  to  have  the  same  position. 

Calculations  of  the  relative  ion  concentrations  present  were  made 
using  the  ratios  of  the  relative  area  values  obtained  in  the  computer- 
analysis  output  (Appendix  II).  Quadrupole  splittings  were  obtained  as 
the  velocity  difference  between  the  two  associated  peaks  and  the  Isomer 
shift  values  were  obtained  as  the  difference  between  the  center  of 
gravity  of  the  two  peaks  and  zero  velocity.  The  final  results  are  shown 
in  Table  III  and  Figures  11  and  12.  In  the  case  of  both  the  di-valent 
and  the  trl-valent  ion,  the  quadrupole  splitting  increases  with  an  in- 

3+ 

crease  in  the  Fe  /Fe  concentration.  Although  the  magnitudes  of  the 
tri-valent  splittings  are  only  about  30X  of  the  di-valent  values,  the 
two  curves  appear  to  have  an  almost  identical  second  order  dependence. 

The  similarity  does  not  appear  in  the  values  for  the  isomer  shift. 

3+ 

While  the  tri-valent  shift  increases  with  the  Fe  /Fe  ratio,  the  di¬ 
valent  shift  seems  to  remain  constant  over  the  entire  range  of  concen¬ 
trations.  A  more  complete  picutre  of  isomer  shift  and  quadrupole  split¬ 
ting  activities  could  have  been  obtained  using  more  data  points  in  the 
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extreme  high  and  low  ends  of  the  range  of  concentrations,  however, 
these  samples  were  not  available  and  were  omitted. 

Low  Temperature  Data 

A  spectra  was  obtained  for  each  of  the  concentrations  at  6  K 
and  are  shown  in  Figures  13,  14,  and  15.  In  each  spectrum  four  inde¬ 
pendent  Lorentzian  line  shapes  were  fitted  to  the  data.  The  positions 
of  these  lines  are  plotted  with  respect  to  the  concentration  in  Figure 
16.  As  the  tri-valent  ion  concentration  is  decreased,  an  almost  linear 
migration  of  the  four  spectral  lines  away  from  a  point  slightly  positive 
of  zero  velocity  occurs.  Accompanying  this  migration  is  a  slightly  non¬ 
linear  increase  in  the  line-width  ranging  from  a  value  equal  to  that 
obtained  in  the  room  temperature  observations  for  high  Fe^/Fe  concen¬ 
trations  to  nearly  six  times  as  large  in  the  low  Fe3+/Fe  concentration 
(see  Figura  17).  As  the  line  positions  expand  and  the  line  widths  in¬ 
crease,  there  is  decrease  in  the  amount  of  effect  that  is  observed. 

This  is  expected  since  the  integral  under  the  spectrum  is  indicative  of 
the  total  iron  present,  however,  differences  in  this  area  and  that  ob¬ 
tained  at  room  temperature  would  indicate  a  temperature  dependence  in 
the  recoil- free  fraction.  The  amount  of  effect  that  is  observed  becomes 
important  when  it  is  considered  that  the  amount  of  noise  present  in  all 
the  data  remains  fairly  constant  at  approximately  0.5%  absorption  while 
the  amplitudes  of  the  fitted  peaks  falls  as  low  as  0.65%  absorption.  At 
this  level  lines  are  quite  difficult  to  resolve  and  multiple  lines  are 

often  completely  obscured. 
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Hast  Treated  Data 

Mossbauer  spectra  were  obts Iced  at  root  temperature  of  Maples 
u 

of  the  SOX  Fe  /Fe  concentration  glass  which  were  hett  treated  for  45, 
90,  180,  and  360  minutes  at  600°C.  These  spectrt  are  shown  as  Figures 
18,  19,  20,  and  21,  respectively.  The  spectrua  of  the  glass  without 
heat  treatment  ms  Included  previously  as  Figure  9.  No  atteapt  was  aade 
to  coordinate  the  variables  associated  with  any  two  lines  for  the  rea¬ 
son  that  no  clear  cut  trend  was  discovered.  However,  a  graphic  display 
of  the  line  positions  and  relative  amplitudes  Is  shown  as  Figure  22. 
Lines  are  numbered  In  order  of  their  appMrance  as  the  tlae  of  heat 
trMtaent  progresses.  Examination  of  Figure  9  shows  that  line  1  was 
Initially  considered  to  be  a  superposition  of  the  two  lines  resulting 
from  the  doublets  of  lines  2  and  3.  While  this  probably  remains  true 
for  the  heat  treated  samples,  it  Is  expected  that  lines  4  and  5  are  due 
to  a  different  mechanism  so  that  a  redundancy  in  the  fitting  procedure 
was  not  altogether  necesMry. 

As  in  previous  experiments,  the  line  Intensifies  must  decrease  as 

more  lines  appear  to  conform  to  the  fixed  amount  of  iron  present.  This 

Is  consistent  with  a  chemical  analysis  which  was  completed  on  all  heat 

42 

treated  samples.  It  was  found  that  the  heat  treating  did  not  affect 
3+ 

the  Fe  /Fe  ratio  to  within  an  experimental  error  of  ±  0.5X.  Also,  It 
Is  noticed  that  as  new  lines  appear  and  become  more  completely  resolved 
the  value  of  the  fitted  line  half-width  decreases.  This  would  tend  to 
Indicate  that  each  of  the  initial  lines  is  the  sum  of  two  or  more  small¬ 
er  lines  in  very  close  proximity.  As  the  time  of  the  heat  treatment  is 
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Increased,  the  Individual  lines  bacon*  mor*  apparent  and  display  the 
natural  line  half-width.  It  la  reasonable  to  assun*  that  since  the 
final  half-width  value  for  the  heat  treated  aanple  compares  favorably 
with  the  value  before  heat  treating,  that  moat  of  the  superpositions 
have  been  disclosed.  However,  a  more  drastic  heat  treatment  using 
larger  amounts  of  the  sample  at  a  longer  time  and  at  the  aaae  tempera¬ 
ture  may  yield  other  Interesting  Information. 

Chemical  Analysis 

34* 

As  previously  mentioned,  the  ratio,  Fe  /Fe,  may  be  obtained  as 
the  ratio  of  the  trivalent  doublet  area  to  the  area  of  the  entire  spec 
trum.  Since  the  accuracy  of  this  method  has  been  verified  by  quantita 
tive  chemical  analysis,  it  was  decided  to  use  the  concentration  value 
obtained  from  the  line  area  ratio  throughout  the  analysis  of  all  the 
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CHAPTER  V 

DISCUSSION  OF  RESULTS  AND  CONCLUSIONS 

In  an  amorphous  Maple  with  a  high  concentration  of  a  paramag- 

netlc  Iona,  aa  la  preaent  In  the  55-45  uole  X  FaO-P^Oj  glass,  nuaeroua 

fora*  of  Fe  Iona  may  be  preaent  almultaneoualy  making  Interpretation 

of  the  Mossbauer  spectra  quite  difficult.  Retaining  the  assimptlon 

made  In  Chapter  4  that  all  of  the  Ion  states  are  In  equilibrium,  either 

Ion  may  act  aa  the  predominate  vltrifler  with  the  formation  of  stable 

groups  of  either  the  divalent  or  the  trlvalent  type,  dependent  on  the 

amount  of  Fe20  present.  Considering  ligand  field  stabilization  In  3d6 
2+ 

Ions,  Fe  might  be  expected  to  display  octahedral  coordination  In 

2+ 

glasses.  However,  since  It  Is  known  that  Ni  occurs  tetrahedrally  co- 

2+ 

ordlnated  In  some  glasses,  the  possibility  exists  that  Fe  may  act 

similarly.  With  regard  to  the  high  concentration  of  Iron  In  the  samples 

subject  to  the  present  discussion,  one  could  expect  numerous  interac- 

3+  3+  2+ 

tlons  between  neighboring  Fe  Ions,  between  Fe  and  Fe  Ions,  and  be- 

2+ 

tween  neighboring  Fe  Ions.  Heat  treatment  could  result  In  the  perci- 

43  44 

pitation  of  several  crystalline  oxides  as  shown  by  Asm,  Bamford,  and 
27 

Belyustin  et  al.  Finally,  It  Is  recalled  that  these  oxides  may  al¬ 
ready  be  present  In  the  colloid-disperse  state.  This  state  is  most  pro- 

44 

bable  at  high  iron  concentrations. 

Examination  of  Figures  8,  9,  10,  11,  and  12  suggests  that  the 
overall  structure  of  the  glass  Is  not  greatly  affected  by  concentration 
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MCIO  Of  either  «.c.inM  T<b|H  ,  >nd  lt 

>t»  IM*  shift  value.  of  approximately  0.10  and  0.020-0.045  for  the 
and  .rival™,  doubl.t.  or.  consistent  with  other  phosphate 

“**“•  *”  •'  oofohodrol  coordination  for  both  Ion. 

fact  that  the  d? ahlf ,  remain.  invariant  with  ,  .angea  l„  Con- 
centration  would  tend  to  show  that  the  electronic  structure  „..ocl.,ed 
with  th...  lona.  particularly  their  4.  wav.  functlona.  are  relatively 

free  fro.  co-v.l.„Cy  change,  due  to  th.  concentration  of  either  Ion. 
Hovever.  the  d.cre.a,  iso...  .hi,t  ,he  trlv.1,nt  ^  ^  ^  ^ 

ere...  in  the  F.^/F.  ratio  augg.at.  an  increase  1„  the  4.  electron 
denalty  which  wold  eas.ntl.lly  amount  to  an  Increase  l„  covalency 
...eclated  with  the  Fe*  l„„.  On.  might  therefor,  surmise  that  any  co¬ 
valency  or  coupling  effects  that  are  being  observed  would  Involve  only 
neighboring  trlvalent  This  Is,  however,  an  Idea  not  completely 


consistent  with  the  quadrupole  splitting  data. 

The  experimental  ranges  of  quadrupole  splitting  values  of  0.046- 

0.01)7  and  0.210-0.270  cm/sec  for  the  Fe3*  and  f«2+  < 

cne  *e  and  Fe  ions,  respectively, 

are  In  good  agreement  with  the  values  of  Kurkjian  and  Slgety24  of  0.032- 
0.088  and  0.200-0.229  cm/sec  for  octahedral  coordination.  This  does 
verify  the  earlier  conclusions  concerning  the  coordination  states  of  the 
two  ions.  However,  It  is  noticed  that  both  the  curves  of  Figure  5-4 
vary  in  a  quite  similar  manner.  It  is  recalled  that  the  electric  field 


gradient  for  the  ferrous  iron  case  Is  due  to  s  5d  electron  configura¬ 
tion.  whereas  the  field  gradient  In  the  ferric  Iron  with  a  6S  configura¬ 
tion  is  due  primer/' y  to  nearby  lattice  contributions.  The  similarity 


M 


In  th*  dtpMdwict  of  tho  quadrupole  splitting  on  tho  Po^/Fo  ratio  of 
tho  two  ion*  night ,  than,  ho  a  roault  of  a  loaaor  intoratonic  diotanco 

34, 

in  highly  couplod  ion  pairs  in  a  mixture  rich  with  the  analler  Ft*  ion. 

One  could  then  conclude  that  noat  of  the  ion  pairs  at  a  antifcrronagnet- 

ically  cooplex  Fs^-Pe2*,  in  agreement  with  Frlebele,  Wilaon,  Dorier, 

and  Kin ear,* 2  who  atudied  theao  glaaaea  by  aeans  of  aagnetic  resonance 

technique.  The  epectra  at  78°  are  more  difficult  to  Interpret.  While 

work  done  on  Fe(POj)j  glaaa  with  0. 1-5.0  mole  X  FejO^fBOX  Fe  /Fe)  by 

Kurkjian  and  Buchanan  showed  little  difference  in  epectra  obtained  at 

78°  and  300°K,  the  high  iron  concentrations  yield  definite  changes  at 

the  two  temperatures.  In  addition,  it  would  appear  that  the  differences 

are  more  precisely  attributable  to  the  increase  in  the  concentration  of 

the  divalent  ion  (decrease  in  the  Fe^/Fe  ratio).  While  the  curves  of 

Figures  13,  14,  and  15  were  fitted  with  four  lines  with  large  values  of 

34* 

half-widths,  it  is  conceivable  that  the  50%  and  15%  Fe  /Fe  spectra  are 

3+  24- 

rat  her  a  combination  of  central  Fe  and  Fe  doublets  and  a  6-line  hy¬ 
per  fine  set.  Since  the  glass  is  known  to  be  predominately  paramagnetic 

26  o 

this  would  resemble  the  results  of  Kurkjian  and  Buahanan  at  4  K  which 
were  explained  in  terms  of  either  a  short  range  magnetic  alignment  or  a 
long  electronic  relaxation  time  which  would  remove  the  nuclear  spin  de¬ 
generacy  and  lead  to  Hfs.  More  experimentation  with  various  other  iron 

O  A. 

concentrations  and  other  Fe  /Fe  ratios  at  more  temperatures  would  be 
valuable  in  resolving  this  question. 

Heat  treatment  of  the  glasses  showed  the  appearance  of  new  lines 
which  may  be  interpreted  as  a  large  central  doublet  superimposed  on  at 
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1m st  on*  6-1 ln»  hyp* r fin*  structure.  Tbs  slallsrity  of  this  4st*  to 

that  obtained  by  Bslyustln  #t  al.27  for  ths  Inverted  spinal  structure 

of  F*j<>4  and  the  data  for  the  crystal  FePO^  fron  Figure  4  Indicates 

that  a  Mixture  of  these  two  crystalline  phaees  Is  being  observed.  This 

42 

wan  also  the  conclusion  of  Prlebele  et  al.  on  the  basis  of  x-ray  dlf- 

34* 

fraction  neasurenents  Heat  treatment  of  the  13X  Fe  /Fe  was  not  exam¬ 
ined  as  this  would  probably  produce  other  crystalline  phases  further 
coapl testing  the  spectrn. 
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TABULATED  COMPUTE*  ANALYSIS  DATA 


MOSSDAUFR  SPECTRUM  PARAME TERS 


TITLE*  FOLDED  SPECTRA  OE  .001  IRON  FOIL  AT  ROOM  TEMP 


♦INITIAL 
(OBTAINED  FROM 

VELOCITY (CM/SEC) 


1 1 

-.6600 

2* 

-.3850 

3* 

-.1210 

4* 

.0550 

5* 

.3190 

6: 

.5940 

ESTIMATES* 

••SEARCH"  INPJT) 

AMPLITUDE  (7.) 


1 1 .5442 
9.6374 
5.9226 
6.7412 
9.5457 
12.1306 


USER  ESTIMATE  OF  LINE  HALE  WIDTW=  .0400  CM/SEC 

...ADDITIONAL  CHANNEL  CORRECTIONS  ARE* 

10#  II#  32#  43#  45#  97#1R5# 197#250#251 # 

269 #271 #272 #31 4# 31 6# 322 #335# 


--N0  VARIABLE  CORRELATIONS  REOUESTED- 


♦AETER  ITERATION  NO.  4* 


VELOCITY  (CM/SEC)  AMPLITUDES) 


RELATIVE  AREA 


l:  -.6553 

Pi  -.3236 

3*  -.1222 

4*  .0614 

5  *  .3225 

6  *  *5928 


1 1 .8303 
9.1199 
4.9096 
5.7212 
9.1054 
12.1666 


2.4096 
1 .8576 
1 .0000 
1 . 1 653 
1  .8546 
2.4781 


FIT  VALUE  OF  HALF  WIDTH= 

MEAN  SOUARF.D  ERROR  OF  FIT  PROCEDURE  = 
DATA  NOISE  LEVF.L  = 


.0452  CM/SEC 
.4876  Z  ABSORPTION 
.3440  %  ABSORPTION 
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»•****<*«>  nrrcfwu*  parameters 


titifi  fojofo  spectra  if  sot  rc  3  /*"  sample  AT  * 


♦  INITIAL  ESTIMATES* 
(OBTAINED  FROM  FXTERNAL  INPUT) 


VELOCITV<CM/SEC> 


AMPLITUDES) 


I t  -.0464 
•>1  -.0160 
31  .053 A 
4<  .Pill 


6.A350 

P4.3370 

P4.3370 

6.R3SQ 


USFR  ESTI MATE  OF  LINE  HALF  WIDTH-  .0400  CM/SEC 

...ADDITIONAL  CHANNEL  CORRECTIONS  ARE* 

I  HR. 


♦VARIABLE  CORRELATIONS  USED  IN  FIT  PROCEDURE* 


DEPENDENT  INDEPENDENT 


RATIO 


AMPC  1)  AMP<  4) 

AMP C  8)  AMP<  3> 


I  .000 
I  .000 


♦AFTER  ITERATION  N0.  5* 


VF.LOC  I  TY<  CM/SEC)  AMPLITUDES)  RELATIVE  AREA 


It  -.0307  5.P661 
r>t  -.0134  5.1548 
It  .0531  5.1548 
At  .8063  5.8661 


1  .0816 
1 .0000 
1 .0000 
1  .0816 


FIT  VALUE  OF  HALF  WIDTH- 

MEAN  SOUARED  ERROR  OF  FIT  PROCEDURE- 

DATA  NOISE  LEVEL- 


.0387  CM/SEC 
.8614  X  ABSORPTION 
.3301  X  ABSORPTION 
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*4 


MfMffMI 


VUiri  *f  |M  M  v»!  » 


•INITIAL  FMI*Af*M 

<9*fAlwro  t*t*  |9*yfl 


VEL0CITt«C*>1l«CI 


A*»MT*jOMt* 


it  -.opoo  7*14*0 

At  •  .01*30  i.ost© 

31  *0033  I *0390 

41  .9014  M4I9 


y*FR  m  LtNf  Mlf  VIOIW*  *0)00  C*^WC 

•  •  .AOOITIOMAL  CH4WIII  CPtfCTtlOS  4lfl 
1 1 3# 1 15*1 1 7# Iff *}3I #133# 


•VARIABLE  CORRflATIRNS  UKO  IN  EIT  PROCEOURE* 


DF*FMOfRT 

INOEpENOENT 

RATIO 

A  MFC  1) 

A MFC  4) 

t  •000 

VEL<  2) 

VEL<  II 

1*000 

A MFC  2) 

AHP<  3> 

1*000 

•AFTER 

ITERATION  NO* 

4* 

VFLOC ITT!  CH/5EC  ) 

AMPLITUDE!*) 

RELATIVE  AREA 

1 1 

-.0IR4 

7.1554 

5.9345 

2t 

-.01*4 

1*205* 

1  .0000 

3« 

•02*0 

1 *2058 

1 .0000 
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*  19*7 

7.1556 

5.9345 

FIT  VALUE  0F  HALF  WIDTH*  .0405  CH/SEC 

MRAN  SOUARFD  ERROR  0F  FIT  PROCEDURE*  .3390  t  ABSORPTION 
DATA  NOISE  LEVEL-  *393P  *  ABSORPTION 


4* 


*■•**•••*•'•■•*■•********•■*** 

nf. f  i  wn**  #*  »**  *?  r1^  # 


*** 

•*•«»•■■■**••*•• 
It  **04*t 

9|  *,^114 

>i  t^i»i 

f,i  •ooo* 


».w 

19*1940 

l,W) 

f.MM 


...  WlflMI 
^14. 


9#  114#  »4|i 

r*Mn*i  ^owm*** 


,©VV4  rwrffC 


or*f«i*oT 


|U<herfV}C>f 


*4f|© 


4*9*t  »l  ** 

4#(  91  3» 


I  •090 
1*900 


•irfft  HfMTIH  «.  M 


V^i  **r* 


|l  «.^4ti 

9| 

>1  *044f 

41  •910% 


9*  .439 
4.MM4 
f. «09% 
9.sm 


1*0000 

3*4«94 

0*9904 

1*0000 


9 IT  V4U#f  #9  **4lf  MIOTW* 

9940  soyooto  wmoo  §9  fit  ouw&wtf* 

04T4  0*149  irvri* 


.0399  C9/S9C 
,4444  1  49SIW9TI** 
■  9940  t  4§S*49T!<*4 


M. 


nmu*  jpr  nw'CTifj*'  HARAKTT***r, 


TiTir*  rjii  «rn  opfot^a  pf  »v  ^V»v  "•akpi f  mfat  th^m^o-as  ^in 


♦  ini 
cortainfd 

vri.#ciTYCCM/srr> 


| t  «.0|9A 

•>  t  .HAG  A 

ot  .?ir 

•JOFR  FOTIMATF  flF  I.INF  HALF 


t  A  |  FSTIKATF5* 

FROM  ••1FAf?nH‘»  INPJT) 

AMPUT'JOF<X> 


9.A97A 
5. 59 1 1 
AtAOin 

nihtm-*  .nooo  ry/^FC 


--N0  VARIAH1.F  CPRRFI.ATt  -4NS  RFO'JFSTFO-- 


♦aftfp  itfratirn  NO.  At 


uri.no  |  TYCOM/SFC)  AMP1.I  T'JOFf  Z>  RFIATIVF  ARrA 


It  -.0*^0?  A  .'MO  9 

.0066  A. 0757 

Ot  .?AAO  0.7605 


1  .RA96 

1 .or 00 

I .0000 


FIT  VAI.'JF  OF  HA1.F  WIOTMb 

f/FAN  SO'JARFO  FRRPR  OF  FIT  PROOFOURF  = 

OATA  NOT  OF  I.FVF1  * 


.0A9f>  CM/SFC 
•  P770  7.  ARS0RPT10N 
.5609  7.  AOS0PPT10N 
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MOSSRAUER  SPECTRUM  PARAMETERS 

TiTiE*  Ft;i  nrn  rprrtRA  of  sox  fc  3+/f0  sampi.e  heat  treatfd-90  MIN 


♦INITIAL  ESTIMATES* 
CURTAINED  EROM  EXTERNAL  INPUT) 


1  : 
Pt 
ri  * 

4  i 


VFL0C I  TY  <  CM/SEC  ) 

-.0| A  A 

•  0  3 1  P 
.0133 

•  30PH 


AMPLITUDES) 

7. PI  50 
A .9430 
I . 8 RAO 
5.0000 

.0300  CM/SEC 


USE?  ESTIMATE  OF  LINE  HALE  WIDTH= 

--N0  VARIABLE  CORRELATIONS  REQUESTED- - 


♦AFTER  ITERATION  N0.  A* 


VEI.0C  ITYCCM/SEC) 


AMPLITUDE  (7,) 


RELATIVE  AREA 


It  -.0176 

4. 2372 

7.2340 

Pt  .0330 

3.1 R65 

5.4401 

31  . P045 

P.09P6 

3.5725 

At  .P966 

.  5R57 

1  .0000 

EIT  VALUE  OE  HALE 

WIDTH= 

.0463 

CM/SEC 

MEAN  SQUARED  ERROR 

0E  FIT 

PROCEDURES 

.  1B52 

X 

ABSORPTION 

DATA  NOISE  l.EVEL  = 

.4035 

% 

ABSORPTION 
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M0SSBAUER  SPECTRUM  PARAMETERS 


TITLE*  FOLDED  SPECTRA  0E  SOX  Fe  ^+/Fe  SAMPLE  HEAT  TREATED-360  MIN 


♦INITIAL  ESTIMATES^ 
(0HTAINED  FROM  "SEARCH"  INPUT) 


VEL0C I TY (CM/SEC )  AMPLI TUDE(X) 


1 1 

-.1163 

P.1064 

P : 

.onoo 

7.P136 

3 » 

•  03RR 

3*6660 

a: 

•  PPP  9 

1 .7933 

Sx 

.300  4 

P . 3393 

USER  ESTIMATE 
. . .ADDITIONAL 
P99# 

0F  LINE  HALF  WIDTH® 
CHANNEL  CORRECTIONS  ARE: 

.0300  CM/SEC 

-NO  VARIARLE  CORRELATIONS 

REOUESTED-- 

♦AFTER 

ITERATION  NO. 

44 

VELOCITY (CM/SEC) 

AMPLITUDE (X) 

RELATIVE  AREA 

1  1 

-. 1P3S 

1 .S 1  S3 

| .0000 

Pi 

-.0046 

5.P796 

3.484P 

31 

.03*8 

4.095P 

P.70P6 

4: 

•  PI  P* 

1 .7647 

1*1646 

Si 

.3074 

1 .6838 

1  *U1P 

FIT  VAI.'JF  OF  HALF  WIDTH** 

MEAN  SO'JARFD  FRROR  0F  FIT  PROCEDURE® 
DATA  NOISE  LEVEL® 


.0316  CM/SEC 
.1993  X  ABSORPTION 
.P7PP  t  ABSORPTION 


1 CH 


69 


M0SSOAUER  SPECTRUM  PARAMETERS 


T  |  ti  Ft  FOLDED  SPECTRA  0E  SOI  EC  3+/Ee  SAMPLE  HEAT  TREATE.D-IRO  MIN 


♦INITIAL  ESTIMATES* 
(fintAINED  FROM  EXTERNAL  INPUT) 


VEI.fiC  I TYC  CM/SEC) 


1 1 

- • 1 P03 

Pt 

-.0034 

3 1 

.0359 

4 1 

•  P0 1  3 

S» 

.301  7 

AMPLITUDE**) 


3.66P0 
9.P060 
R • 5930 

3. P370 

4 . P070 


USER  ESTIMATE  DE  LINE  HALE  WIDTH*  .0300  CM/SEC 

...ADDITIONAL  CHANNEL  CORRECTIONS  ARE  I 
I  07* 1  HR# 


--NO  VARIARL*  CORRELATIONS  REOU*STED-- 


♦aeter  iteration  NO.  4* 


VELOCITYCCM/SEC)  AMPUTUOFCX)  RELATIVE  AREA 


It  - • | I  PR  1  .A7S7 
Pt  -.0040  4.60RJX 
St  .0340  4.6060 
At  .POOS  I.67A1 
St  .301!  P.3SPI 


I .1177 
P. 7464 
P.7447 
I .0000 
1 . 40SP 


FIT  VAI  »JE  OE  MALE  WIDTH* 

mean  SoyARED  ERROR  OF  FIT  PROCEDURE*  • 

DAT*  NOISE  (  EVEL* 


lvi!) 


03P9  CM/SEC 
1971  *  ABSORPTION 
4707  *  ABSORPTION 
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M0SSRAUrR  SPECTRUM  PA  RAMP  TERS 


t i t»  r »  poi  nrn  spectra  of  isz  Fe  3*/Fc  sampi  f  at  77  k 


•INITIAt  FSTIKATFS4 
< ORTA  I NFD  FROM  ,,SFARCMM  INPUT) 


WF|  0CITY<CM/SFC) 


AMPI  lTjOrtX) 


1 1 

-.9649 

9t 

-.0944 

It 

.9907 

4t 

•  S9*S 

P.46R7 

9.9S11 

9.17*1 

9.S774 


USPO  ESTIMATE  OF  I. INF  HALF  Wlf)TM>  .0100  CK/SFC 

.  .  .AOOI  Tl  ONAI  CHANNFL  CORRECT I  ONS  ARF* 

1H,  46,  47.  SO,  S|,  S4,  99, 107,144,146, 
144,146. I 41, 10* ,309. 149. 19*. 


--NO  VARIABLE  CORRELATIONS  REOUESTED-- 


•AFTER  ITERATION  N0.  7* 


VEt  0CITY<CM/SFC)  ANPLrrUOECt) 


RELATIVE  ARFA 


1  t 

-.1949 

.9476 

9* 

-.01*1 

1 .409S 

It 

.4190 

•  *9*7 

4» 

.4971 

•  ASA* 

I .4SS9 
9.4799 
I .1*09 
I .0000 


FIT  VALUE  0F  HALF  W!0TM» 

MFAN  SO'JARFf)  ERROR  OF  FIT  PR0CE0URE. 
OATA  NOISE  I  FVF|.« 


•  PS6I  AM/SEC 
•944*  X  ARSORPTI0N 

•  4*11  t  ABSORPTION 


1(0 
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MCSSRA'JFR  SPECTRUM  PARAMETERS 
TITIE:  FOLDER  SPFCTRA  OF  501  Fc  ^/Fe  SAMPLE  AT  77  K 


♦INITIAL  ESTIMATES* 
contained  FROM  F.XTFRNAL  INPUT) 


VFLCriTYCCM/SFC) 


1  I 
?! 

7* 

At 

•jr.FR  ESTIMATE  OF 


•  PE.7A 

•  0P63 
.P779 
.4477 

I.INF  HALF  W I  OTH* 
CCPRFCTIONS 


AMPLITUDES) 


R.0000 

P.0000 

P.0000 

9.0000 

.0700  CM/SEC 

ARF  t 


•  ..ADDITIONAL  CHANNEL 

77*P73»'>7A#7P4*3A3#3AA# 


--NO  VAR1A0I.E  CORRELATIONS  RFOUrSTEO 


•AFTER  ITERATION  NO.  0* 


VFI.CC  I  TY C CM/SFC  >  AMPLI  TUOrCf  )  RELATIVE  ARFA 


|l  .,9RAA  | .SANA 
P|  -.0P77  O.SPMP 
71  .f\7IP  | .901 R 
At  .AA3«  I *1705 


I .333P 
P. 1599 
I  .6PAP 
I  .0000 


fit  value  0F  malf  WIDTH. 

Mr  an  SOUARFO  ERROR  OF  FIT  PROCEDURE. 
DATA  NO I SF  level* 


•  IS3A  CM/SEC 
. P.7.79  2  APSORPTION 
.3174  t  ADSORPTION 


ill 
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MOSSRAUER  SPECTRUM  PARAMETERS 


TITLE!  FOLDED  SPECTRA  OF  79t  Fe  ^/Fe  SAMPLE  AT  77  K 


♦INITIAL  ESTIMATES* 
(OBTAINED  FROM  EXTERNAL  INPUT) 


VFLCCITY<CM/SFC) 


AMPLITUDES) 


|t  - .0454 
PI  .01S4 
SI  *09  79 
41  .3140 


I .0000 
S.2R70 
4.5000 
1 .4400 


USER  ESTIMATE  BE  LINE  HALE  WIDTH*  .0300  CM/SEC 

...ADDITIONAL  CHANNEL  CORRECTIONS  ARE* 

54* 109*1  PI # ISP. 109*1 9R# 304*355*396* 


--NO  VARIABLE  CORRELATIONS  REQUESTED-- 


♦AFTER  ITERATION  NO.  A* 


VELOC I TYCCM/SEC )  AMPLITUDES)  RELATIVE  AREA 


It  -.0455  4*1  TAP 
Pt  .0 1 5P  P.9P73 
3,  .0978  3.5774 
4 1  .3140  | »P 195 


3.4P45 
P • 4004 
2.9335 
1 .0000 


FIT  VALUE  OF  HALE  WIDTH* 

MEAN  SOUARED  ERROR  OF  FIT  PROOF  >URE= 
OATA  NOISE  LEVEL* 


,0449  CM/SPC 
*2P?R  X  ABSORPTION 
•  4938  X  ABSORPTION 


1 1'i 
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APPENDIX  II 

mOssbauex  spectrum  analysis  procram 

Th.  computer  progran  U.t«d  1«  thl.  .ppwdl.  -  '*•'«'  *® 

.11  .1  th.  d.ta  obt.lo.d  in  th.  pr..«vt  1  nv.it lg»tloni.  It  U  wrltt« 
in  Fortran  TV  to  h.  —d  with  th.  XDS  Sign.  7  conpotor.  With  th. 
•—option  of  oovoral  HMD  nod  WUTE  In  which  orltl-otlc 

oporotl...  nr.  porfor-d.  th.  progr-  1.  opntlbl.  with  nn, 

.golppod  -th  .  it.nd.rd  USASI  Fortr..  IV  co-pll.r  nnd  At  1—t  *5.000 
word,  of  high  .po.d  —ry.  Typlc.1  running  tin.,  nr.  two  •!">■».  for 
compilation  .nd  1.5  nlnuto.  for  ..motion.  All  input  1.  r«d  on  th. 
FtJ  DC.  .nd  .11  output  l.  wrl.ton  on  th.  Ft*  DC».  If  •  Clct-pu  pl.t 
of  th.  probl*.  1.  d..l«d.  It  1.  ohtnlMd  hy  ...lg«I«g  th.  FiFLOT  DC. 

to  th*  drwo  plotter. 

A  Mlutlon  1.  obtained  bp  eotr-tlng  th.  lnltl.l  ..tl«»to»  for 
th.  p.r«n*t.r.  of  u»  to  10  Hparat.  w.ctr.1  lino,  to  glv.  •  >•»•*- 

boat  fit  .ppro.ln.tlon  to  th.  ..p.rln.nt.1  d.M.  lnltlol  .«l- 
„t„  g.  .«pll.d  olthor  lnt.rn.1  or  o.t.rn.1  to  th.  progrr-  ond 

th.  output  1.  in  th.  font  Of  *PP«~><*  «•  Additionally.  th.  lln.  p.« 
n.t.r.  n.y  b.  «d.  lnt.rd«»nd«t.  U-  of  th.  pr.gr-  r.dulr..  th. 

following  Input  card*: 

l«t  card  -  (Format  II)  plott*r  option 

(a)  0,  no  plot  deaired 

(b)  l,  plot  on  device  aaaigned  to  FrPLOT 
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2nd  card  -  (Format  12,  15A4)  number  of  characteris  In  the  title 
atrlng  (no  more  than  20),  title  atrlng 
3rd  card  -  (Format  215,  2F10.5,  II)  number  of  points  in  raw 
data  (no  more  than  400),  number  of  peaks  to  be  fit 
(no  more  than  10),  approximation  of  line  half-width 
(in  cm/sec,  if  negative  then  not  varied  In  fit  pro¬ 
cedure),  velocity  scale  (in  cm /see,  positive  or 
negative),  option  code: 

(a)  0,  no  options  requested 

(b)  1,  options  requested  as  listed  on  next  card 
4th  card  -  (If  option  code  is  1)  -  (Format  612,  F10.5)  options 

in  order  left  to  right  are  marked  with  1  is  requested 
0  or  blank  if  not  requested,  or  aa  specified  below: 

(a)  ntaber  of  bod  channels  to  be  corrected  (default 
to  0) 

(b)  printing  of  each  iteration  of  fit  procedure 
(default  prints  only  last  iteration) 

(c)  folding  code:  -1  for  data  in  left  half,  +  1  for 
data  in  right  half,  0  if  data  is  to  be  folded 
about  middle  channel  (default  to  0) 

(d)  search  procedure  which  generates  Initial  para¬ 
meters  for  the  number  of  spectral  lines  specified 
previously  (default  requires  estimates  to  be 
Input  separately) 

(e)  parameter  correlation  option.  Input  no.  of 


lit 
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correlations  desired.  Actual  correlations  listed 
below  (deiault  is  no  correlations) 

(f)  maximum  number  of  Iterations  allowed  (default 
Is  99) 

(t)  convergence  criterion  (default  Is  0.0001) 

The  following  cards  supply  additional  information  for  the  execution  of 
options  (a),  (d),  (e).  The  additional  data  cards  for  an  option  must 
appear  If  that  option  is  reguested  and  in  the  order  listed  below: 

(1)  channels  which  are  to  be  corrected  are  listed  on 
successive  cards  (Format  1015).  The  count  value  for  a  corrected 
channel  is  made  equal  to  the  mean  count  value  of  tha  two  adjac¬ 
ent  channels.  Adjacent  bad  channels  nay  not  be  corrected. 

(it)  this  card  must  appear  regardless  of  the  nwber  of 
options  requested.  It  defines  the  Format  in  tfhlch  the  raw  data 
is  to  appear  later.  The  Format  appears  as  a  legal  FORTRAN  IV 
statement  with  enclosing  parentheses. 

(iii)  Correlations  are  expressed  in  terms  of  an  inde¬ 
pendent  vcriable,  a  dependent  variable,  and  a  multiplicities 
ratio  which  relates  then.  A  spectral  line  is  specified  by  two 
parameters,  position  and  amplitude.  Considering  (N)  lines  or¬ 
dered  by  ascending  velocity,  the  line  parameters  are  numbered 
from  1  through  (2N)  such  that  the  odd  numbers  represent  velocity 
parameters  end  the  even  numbers  represent  amplitude  paraaeters. 
Each  correlation  la  placed  on  a  separate  card,  (Format  215, 
F10.5),  giving  the  dependent  parameter  number,  the  Independent 

1  Hi 
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prueter  n«b.r,  end  the  de.lred  r.tlo  (dependent/lndependent). 
For  «upl>.  constraining  the  awplltude  of  the  third  Hue  to  bo 
twice  the  amplitude  of  the  fifth  line  would  be  written  ee, 

6  .  .  .  10 . 2.0 

three  perenet.r.  nny  be  correlated  toother  by  correlating  eny 
two  with  the  third. 

(lv)  the  Initial  astlaataa  are  Input  using  a  separata 
cerd  for  eech  petrel  line.  Th.  poeltlon  (cw/eec)  end  enplltud. 
(0.01  n  percentage  ebeorptlon)  ere  written  u.ln,  .  rotoot 

(2F10.S). 

„  .11  ce...  the  n-ber  of  entr.  deu  «1«~  or  cerd.  -net  coincide  with 
the  value,  found  on  the  option  cerd.  l«dl»t.ly  following  the  U.t  o 
the  entr.  date  cerd.,  th.  re.  deu  1.  ed*d.  «l«t«  to  the  *>™t  that 

waa  previouely  spaclifled. 
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0 1  MENS  I  0NX  (P0P)*Y<P0P)  *  C ( 22 )  *  NAME  C  1  5)  *  I  BUF(  1000>*DATA( 
$500 ) *KINDC 1 5 

A)*  ICOCPOO)*  IF0RMTC  10)  .CMATC2P *22 > 

C0MM0NNAMF  #M» IPT* MVAR.CMAT 
C0MM0N/F  J  /X»Y*C  *  OF!.*  N.  MAX  .GAMMA 

0ATAF1N0/' I 5T  S^AIF***  F0  *  »  '  LORD  *  *  'PNO  '  *  '  HALF  *  »  *  S 
SPF  '  * 

A  'CTRA  '  *  '  0F  '  *  '"SPA  •  *  '  RCWM  ' »  '  F.XTF  '  *  'RMAL  '  *  'VEL  (  '  *  *  AMPC 

f.  •/ 

CAt.L.AiJTORX 
RFA0C5*  1  1«0>  IPI.0T 

IFC IPL0T.FO.1 >CALLPL0TSCIOUF* 1000*1 > 

NAMF.C  n>*KlNDC7> 

N  A  MF  (  *1 )  =  K 1  ,M0  (  2  ) 

NAMF.  C  5  >  =  K  I  N0(  9) 

1000  prAr>(5*  1P70*F.\0=1  1  70>NCHAR*  CiMAMFC  I  >*  1=6*  1  S> 

PFADC  S*  1  P90  )  iMPTO*  NPK  *  GAMMA  *  VFL*  OPTION 
C  CP4  NPK*  1  ) =A0SC  GAMMA  ) 

CCP*NPK*P>=0. 

M=2+NPK<*P 
0010101 1=1 ,M 
DCIOlft.IJsl  *M 

TFCT I  .NF..I,I>CMAT<I  I*JJ>=0. 

1010  f  F C  1 1  .FO.J.J)CMAT<  I  T  ,  J,|)s|  . 

NR A 0=0 
1  PT  =0 
IFPI  0=0 
MAX *99 
ofl*  .ooni 

1  SFR0H*O 

IFC  OPT I  0N  •  GT  .0  >  RF.ADC  S*  1 .100  )  NR  AD*  I  PT  *  IF01.D  *  I  SERCH*  IC0RR 
t.MAX.OFL 

IFCNRAD.NF .0>RFA0C5*  1310)0  CR  Cl  I>*  11*1  *NBA0) 

RFA0C5* 1  PRO) IF0RMT 

I FC 1 SFRCH  »NF .0 )RFADC  5*1 190)CCC2*I I -1 > *C CP* !!> *  If  =  1 *NPK 

*> 

ifcicprr.fo.o.) goto  1030 

D01 OPO 1C* 1  *  f C0RR 

RFAO C  5  *  I  POO  »  JOFP  *  I NOFP .RATIO 

CMATC  .IOFP.JOFP)  *0  • 

CMATC  .IOFP*  I  NOFP  >  -RAT  1 0 
IOPO  CMATC  INOFP*.IOFP)*RATI0 
1010  RF.AOC*V*IF0RMT)COATACI)*I  *|  #NPTS> 

I F  CNRAO  *F.O  .0  )  GOT  01 070 
1040  0010601*1 *NRAO 

1FC  icon  >  .Fo.1  .0R.1CRCI  >  .FO .NPTS> G0T0I 050 
C  IFCICRCD.FO.ICRCIO  >-l  >G0  TO  SO 
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1050 


10  70 


IORO 


DATA< ICR< I >>  =  <DATA< ICR< I >-l> ♦DATA? ICBC  1 >  +  !>>*.* 

S>>*.5 

1060  C0NTINUF! 

C0NTINUF 
NRAD  =0 
NX=NPTS-1 
00 1 ORONY=P  *  NX 

A  VGF.=  <  DATA  <  NY  - 1  7  +D  ATA  <  NY  + 1  >  >  *  •  s 

IF<ARSCDATA<NY)-AVGF.)  .LT.AVfiF.+  .OUG0T0tORO 

NRAD  =NBAD+ 1 
NRADC=NBAD 
ICB<NRAD)=NY 
C0NTINUF. 

I F (NRAD . GT«0>  G0T0 1 040 
DV=.$ 

Cl.=  l  . 

CR  =  t  • 

IF< IF0LD  »LT •0)CR*0 • 

IF<IF0LD.GT.O>CL=O. 

I F< IF0LD  «NF •0>DV=1  • 

N=(iMPTS+l  >/R 
0010901  =  1  *N 
JP  =NPTS  + 1 ~ I 

Y ( i ) =OV*<CL*DATA< I ) ♦CR*DATA <  JP ) ) 

JO=C IF0LD+1 >*S+1 
NAMFC 1 )=KIND< JO) 

NAMF.(P)  =KI  ND<  JQ  +  1  > 

of.i.v/=p*arso/fl>/n 
XC1 )=-<AnS<VFL)> 

001 1001 =2»N 
XCI )=XCI-l ) +DFLV 
ymax=y< 1 > 

001 1 101=1 »N 

IF<YCI).GT.YMAX)YMAX=Y(I) 

001  IPO  I =1 *N 

IfIUf.rchIfq^callscarchcx^y^c^n^npk) 

Um<r.“NAMr<,>......S,.K.ND<lSF.RCHMO,.K.ND<.SE 

m<”*.'jk-!>.ioo.*c<2«jk>.  JK=I  .NPK) 

^fInRADC  •  WRITER  6*  1A30MICR<IY>»!YSI»NBADC> 

IF< IC0RR  *NF .0) G0T01 1 30 


090 


100 


l  10 


IPO 
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WRITF<6#1240) 

OOTOl 1 60 

mn  writf<6,i250) 

nit mnic*i #m 

ifcckat<ic#ic>.eo.i )oct«i iso 

DOI  MOM'I  #M 

I FCCMATC IC#IJ>»F0.O») GOTO I  I  40 

IP1*0 

IP2*0 

IFCMOOC IC#2).F0.0.)1PI*| 

I F( MOOC I  J»2 ) .EQ.O.  > 1PP*| 

ICl*CIC*l)/2 
I Jl *CIJ»I )/2 

WR?TEC6#1260)KINO<I4*1PI  )  #  ICl  #KINO<  I  4«IPP)  #  1.J1  #CMATC  IC 
*.  1.1} 

1140  CONTINUE 
1150  CONTINUE 
1160  CALLFIT 

i  F(  IPLOT.FO.  1  )CALLORAPH<Y#X#C#NAME#N#NPK#  IBUF#NCWAR) 

'.0T0IOOO 

1  1  70  IFCIPf  ' T.F0.1 )CALLPL0TC0.#0.#99V) 

CAl.LEXI  F 
I  1 RO  F0RMAT<  J  1  ) 

1190  F0RMAT<2F10.5) 

1200  F0RMATC2T5.FIO.5) 

1P10  F0RMATC ✓//✓/✓✓22X#29(  '  -  * )/#22X#  'M0SSBAUER  SPFCTRUM  PAR 
JAMETERS  *  /#22 

4X#29('-')//#5X#'TITLE*  ' # 1 5A4///#28X# '•INITIAL  ESTIMAT 
SES*'/#23X»  *C 

A OBTAINED  FROM  *#2A4#'  INPUT) '//# 1 6X# 'VELOCITY (CM/SEC) ' 
S#  I 4X# 

4 'AMPLITUDE CX) ' /#  1  6X#  1  6(  '  - ' )  #  1  4X#  1 2( *•')/#  C I 2X# I2#'*'#4 
SX#FR.4#2IX# 

AFP. 4)) 

I  PRO  F0RMATC/5X. 'USV  R  ESTIMATE  0F  LINE  HALF  WIDTH* ' #3X#F7.4 
S,  •  CM/SFC') 

I  POO  F0RMATC5X# '...ADDITIONAL  CHANNEL  CORRECTIONS  ARE* '/CIO 
!X#I0CI3# '# f ) 

A)) 

|  240  FORMAT  I //l  7X,»  *  --NO  VARIABLE  CORRELATIONS  REQUESTED--') 
1250  F0RMATC///I2X# '*VARIABLE  CORRELATIONS  USED  IN  FIT  PROC 
SFDURE* '//l 3X 

4#  'DEPENDENT'  #RX.‘  '  INDEPENDENT*  #  1  OX#  'RATIO  '/#  I  3X#9C  •  -  ' )  # 
JRX#1  1 ('-')#! 

A0X.5C '-')) 

1260  F0RMATC14X.A4.12# ') '.l IX#A4#I2#'> '#|0X#F7.3) 

1270  FORMATC 12# | 5A4) 
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1 280  F0RMATOOA4) 

1290  F0RMAT<2I5*2F1O.5*I 1> 
1000  F0RMAT (6I2*F10»5) 

1310  FORMAT  (  I  0  1 5  ) 

F.ND 


C 


SUB  R  OUT  1NF.SF  ARCH  (X  *Y  *C*N*NPK) 

DI  MF.NS  I  0NXC2O2)  *  Y  C 202  >  *C  (22  )  * DUM<  202  ) 
D02OOOI *1 *N 
2000  DUM<n=Y<I) 

nvar=2*npk 

D 020201 =2*NVAR*2 


C( I)=0. 

0020  1  0  J  =  1  *N 

IF(C(I) •GT.DUM(J))G0T02OtO 

c(I)=dum(J> 

C  < I -1 >aX(J> 

2010  CONTINUE 

002020 J=1 »N 

2020  OUM< J) =DUM( J) -RNTZ (X( J> #C<I)»G<I“t 
0020  50  K  =2  *  NVAR  * 2 
IF(NPK-1 ) 20 30 *20 60  *2030 
2030  O02O5OL=4*NVAR*2 

IFCCCL.-1  )  .GE.C(L-3>>G0T02O4O 
OP =0  CL -3) 


*C(NVAR+1 


CA=C(L -2) 
CCL-3) aOCL-1 ) 
CCL-2) =CCL) 


C(L-l) =CP 


C(L)=CA 
2040  CONTINUE 
2050  CONTI NUE 
2060  RETURN 
END 
C 


) 


C 


3'JRRO'JTINFFIT 

01  MENS  I 0NA<22*22) *B(22) *CC22) 


*P (202  *22 )  *R (  202 )*X(202) 


IYC202) *NAME( 

A 1 5 ) *ARFA  CIO) *C MAT (22 *22 ) 

COMMON NAME  »M* T  PT  *MVAR  *CMAT 
COMMON/FI /X*Y*C*OEL*N*MAX *  GAMMA 


COMMON /SO /A *B  » MA 


# 


AR=0  . 
MA  =  1 
I  T=0 
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3000  FL=AR 

CAI.LC0RRECTCCMAT  #C#MVAR#M) 

AR  =0  . 

0030101=1 #N 

3010  AR=AR+(Y<  I  ) -VAL<  I  )  )**{? 

IF(  IPT.  GE  •  1  .AND.  IT  ..ME  .0  )  WR I TE  C  6  #  3  1  70  )  I  T#  AR#  (  I  ,CC  I  )  #  I  =  1 
*#M) 

IFC  IT. I  .E  .0)  O0T03ORO 
IF(AR*!(AR-F.L)/AR.LE  .DFI.)  G0T03  I  10 
IF(IT»GE» MAX ) 60T031 00 
30P0  IT=IT-M 

0030301=1 ,N 
Rf  I  )  =VAI.(  I  )  —  Y Cl) 

003030 J-l #MVAR 
P( I # J)  =PART<  I  #.J) 

3030  O0NTINUF 

0030401=1 » MVAR 

n  ci  )=o.o 

003040.1=1  ,MVAR 
3040  A (  I  #  J) =0 .0 

003060  I  =  |  »  MV/AR 
003060.1=  1  »N 
D0.3OSOK*I  #MVAR 

3030  A(I#K)=A(I#K)+P(J#I)*P(,J»K) 

3060  R<I)*B(I)-P(J,I)*R(J) 

3070  CALL30LVE(  I. IK) 

G0T0<3O3O#31 10)# IJK 
30  P0  IP0INT= I 

0030901=1 #M 

IF<CMAT< I # I) .NF. 1 >G0T03O9O 
Cm=FKIP0INT)+n<I> 

IP0I NT= IP0I NT*1 
3090  CONTINUE 

IF(GAMMA.LT.0)C(M-1 )=ARS< GAMMA) 

G0T03OOO 

3100  WRITE<6#31R0)MAX 
RETURN 

3110  MHALF=(M-P)/p 

003 IPO  I ARFA  =  I #MHALF 

3  IPO  AREA ( I  AREA ) =C<M-1 ) *C <P* I  AREA ) ♦< ATAN< ( 1 00 . -C< P*I ARFA- 1 ) 
$ ) /C ( M - 1 ) ) -AT 

4AN( ( -100 • -C<P*IAREA“1 ) )/CCM-l )) ) 

ARM  I N  =  AREA ( | ) 

003 1 30  I R  =  P  #  MHALF 

31.30  IF(AREA<  IR)  .LT *ARMIN)ARMIN=AREA(  IR) 

D031 401 R=1 # MHALF 
3140  AREA< IR)=AREA(IR)/ARMIN 
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writfc6#3i5o:>it#<  jk#ccp*jk-i>#ioo-4Ccj>*jk>#areac  jk>»  jk 

S«1#MMALF> 

WRITE(6#3I60>CCM-1  >*SORT<AR/N>*100*#C004100* 


RETURN 

3150  F0RKAT<///*4X. ••AFTER  ITERATION  NO-  •#  IP*  •»•//.  10X#  *VF 
5I.0C1TVCCM/SF  .... 

40  *#4X* 'AMPLITUDE**)  **6Xr ‘RFLATIVF  AREA  V*  10X#  I  6<  -  >» 

S4X#|5»C,-,># 

44X*13< •-*>/* AX# IP#* I '#9X*F«.4»| |X#FR.4#I0X#F0.4>> 

3140  FORMAT* /5X# ‘FIT  VALUE  OF  MALE  til DTW* •  •  1PX#FT .4*  •  CM/SE 
JCV.SXi'HrtN 

4  OQUARFO  FRR0R  OF  FIT  PROCEDURE*  *  #F7.4#  *  *  ABSORPTION 
S# /SX# 'DATA  N 

40ISF.  LEVEL*  * #|9X#F7»4#  *  %  ARSORPTION*) 

31  tO  F0RHATC//*ITERATI0NM3/lRC*-*)/*SUW  SQUARED  ERROR  «*IP 
0DI3.R/ 

4//<lOX*C< *OPIP# 1F0I3. 7>> 

Hf f  P0RMAT* //•••••NO  CONVERGENCE  INM3.*  I TERATIONS****  > 
END 


SURROUTINESOLVECIJK)  _ 

DIMF.NSI  ONA*PP#PP)  #R*PP)  #C*0P#P2)  #NAME*  15)#L*2P)  #  WM*J!2) 
*#F.*4S0>#CMAT 
4<PP#PP> 

C0MHONNAMF.#M#  IPT#MVAR*CMAT 

C0MMON/S0/A#R#M 

OOURt.EPRECIRI  0NE#D#C#X 

I. IK-1 

ICOUNT-l 

004000 1 *1 #MVAR 

004000.1*1  *MVAR 

E<IC0UNT>*A<J#I> 

4000  I C  0*JNT *  I C0UNT  ♦  1 

CALLMINV*E#MVAR#D#L#MM> 

ICOQNT*! 

0040101 *1 #MVAR 
0040 10  1-1  #MVAR 
C<J#I)*E<IC0UNT> 

4010  I COUNT  *  I COUNT ♦ 1 
IFCD.EO.O)IJK«f> 

IF*D.EO.O> WRITE* 4*4060) 

Z*0  *0 

0040301-1 #MVAR 
X*0. 

D04OPOJ*! #MVAR 
40P0  X*0*  I # J)*A* J# I )+X 
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4030  7«Z«ABSCX-I *0) 

IF<Z*GT*I  *F.-03>  WR!TE<6#40T0>Z 
IFCZ*OT.1  *F-03>IJK«P 
D04O4OI *1 #*VAR 


A< I # 1 >  *0  *0 
D04O4OJ*! #MVAR 

4040  A<i#n-A<i#n*c<i#.n*B<j> 
004050! *1 » KVAR 
4050  RC!)«A<I#I) 


RETURN 

4080  E0RMATC* INVERSE  DOFS 


4070  formatcpx#  • INVALID 


NOT  EX  I  ST M 
INVFRSF  Z  *  *  #ER  < 


SJ 


END 


c**** 

c**** 


THE  FOLLOWING  SUBROUTINE  SHOULD  CALCULATE  ™E 
OF  THE  APPROXIMATING  FUNCTION  EVALUATED  AT  THE  I -TH 
DATA  POINT* 

FUNCT I  ONVAI  C  1  ) 

DIMFNSI ONNAMEC 1 5J  #X(POP> #Y<BOB>#CCPP>  #CMAT<RB*PF> 
COMM0NNAMF#M# IPT#MVAR#CMAT 
COMMON/F I /X* Y  *C  *DEL#N#MAX  # GAMMA 


L*M-P 
VAL*0 • 

005000 J* I #L#2 
ox«cxcn-ccj>)/ccM-i > 
r ■ i •  /< i .♦nx*DX) 

5000  VAL*VAL*CC J*I >PR 
VAL»VAL*CCM) 


RETURN 

END 


C 

C 

C4*** 

c«-«-** 

c 

c**»* 


THE  FOLLOWING  SUBROUTINE.  SHOULD  CALCULATE  THE  VALUE 
OF  THE  PARTIAL  OF  THF  APPROXIMATING  FUNCTION  WITH 
RESPECT  TO  THE  J-TH  VARIABLE#  EVALUATED  AT  THE 
I -TH  DATA  POINT* 

F'JNCTIONPARTCI  #  J>  _ _ 

OIMENSIONNAMEC  15)  #X<ROB >  #Y<20R>  #C<2P>  #CMATC22#22 > 

COMMON/F I /X #Y#C# DEL #N# MAX# GAMMA 

C0MMONNAME#M#IPT#MVAR#CMAT 

PART=0. 

I  MX  *0 

008000 JREAL* I #M 
IMX=IMX*CMAT<  JREAL# JREAL) 

IF< IMX*EO* J) G0T06010 
8000  CONTINUE 
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6010  IFC JREAL.fO.M)G0T06O5O 

IFC  JREAL.EO.MM >G0T06O6O 
IFCM0DC JRFAL#P> .FQ.O>G0T0AO3O 
D06OROIR0W«1#M*2 

IF<CMAT< I  ROW* JRFAL) *F0 *0) G0T06ORO 
OX*<XCI)-C<IR0W))/C<M-n 


n=i ./< i .♦dxaox) 

PART*PART+f».*C<IR0W4l  >*DX*B*R/C<M-1> 

6020  CONTINUE 


RETURN 

6030  006040 l ROW*? *M*2 

IFCCMAT<IR0W* JREAL).EO.O)O0T06OAO 
OX«CX<n-C<IR0W'l>)/C<M-l  > 
PART.PART4CMATCIR0W.JREAD/C  I  •  ♦DXADX  ) 


6040  CONTINUE 
RFTURN 


6050  PART*1 * 

RFTURN 
6060  L*M-2 

D06O7OK*! #L#2 
DX*<X<I)“C<K>>/C<M-| > 
B=DX/t I •♦OXAOX) 

6070  PART«PART*C<K*1 >*B*R 
part=?.*part/c<h-i ) 
RFTURN 
FND 


SURROUTINEC0RRECT<CMAT*C#MVAR*M> 

OIMENSI0NCMAT<?2#22> »C<?2> 

MVARsM 
0070101=1 #M 

1F<CMAT<I»I).F0.1 • >  G0T070 1  0 
C(I)=0. 

D07OOO J= 1  *  M 

IF<CMAT<I#J> .EO .0 ) O0T07OOO 

n<I)=0<I)+CMATCI#J)*0CJ) 

7000  CONTINUE 

MVAR=MVAR-1 
7010  CONTINUE 
PETURN 
END 
C 
C 

r,  SUOROUTINE  MI NV 

C 

r.  purpose 
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C  INVENT  A  MATRIX 

C 

C  USAGE 

C  CALL  MINV<A#N#D*L#M> 

c 

C  DESCRIPTION  OF  PARAMF.TFRS 

C  A -INPUT  MATRIX*  DESTROYED  IN  C0MPUTATI0N  AND  REPLAC 

SED  BY 

C  RESULTANT  INVERSE 

C  N-0ROER  0F  MATRIX  A 

C  D-RES'JLTANT  DETERMINANT 

C  L-W0RK  VECT0R  0F  LENGTH  N 

C  M-W0RK  VECT0R  0F  LENGTH  N 

C 

R  REMARKS 

C  MATRIX  A  MUST  BE  A  GENERAL  MATRIX 

C 

C  SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  NEEDED 

C  NONF 

C 

R  METH0D 

C  THE  STANDARD  GAUSS -JORDAN  METHOD  13  USED.  THE  DETE 

SRMINANT 

C  IS  ALSO  CALCUATED.  A  DETERMINANT  OF  ZERO  INDICATES 

STHAT 

c  THE  MATRIX  IS  SINGULAR. 

C 

C  . 

51 .  .  .  . 

R 

SUBROUTINE^ I NV<A*N*D#L*M) 

D I  MENS  I  WtxAC  1  >  *L<  1  )  #M(1  > 

C 

G  . 

?!•••• 

c 

r,  IF  A  DOURLF.  PRECISION  VERSION  OF  THIS  ROUTINE  IS  DESIR 

SFD#  THE 

f;  C  IN  COLUMN  1  SHOULD  BE  REMOVED  FROM  THE  DOUBLE  PPRECI 

SSI  ON 

C  STATEMENT  THAT  FOLLOWS 

DBUBLEPRECI SI ONA*D*BIGA*H0LD 

C  THE  DOUBLE  PRECISION  VERSION  OF  THIS  ROUTINE  MUST  ALSO 

C  CONTAIN  DOUBLE  PRECISION  FORTRAN  FUNCTIONS.  ABS  IN  STA 

JTEMENT 

C  ID  MUST  BE  CHANGED  TO  DABS 

R 
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C3  ra  r>  cj  r»  r> 
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C  . 

ft  •  •  • 

c 

r.  RFAPCH  F0P  THF  LARGFST  FLF.MFNT 

C 

0  =  t  .0 

NK  =  -N 

D0«tfiOK  =  t  #N 

NK*NK-»N 

L(K)*K 

M(K)=K 

KK*NK+K 

RIGA=A(KK> 

D0RORO J=K#N 
IZ*N+( J-t > 

D0RO2OI *K»N 
I J  =  I  Z  ♦  I 

R000  IF<DABS<BIGA>-DABS<A< I J>>  >B010#8020*«020 
R010  RIGAsA(IJ) 

L<K)»I 
M<  K> =J 

R0f>0  C0NTINUE 

INTERCHANGE  R0WS 


J=L<  K) 

IF(.I-K>R050»R050*B030 
R030  K I =K-N 

D0RO4O I = I • N 
K I =K I +N 
H0LD=-A<KI  ) 

JI=KI -K+J 
A(KI)=A<JI> 

R040  A < J I ) =H0LD 

INTERCHANGE  C0LUMNS 


ROSO  I=M(K) 

IF<  I  -K ) BORO  # 80 RO * B060 
R060  JP=N*(I-D 
D0RO7OJ=1 #N 
JK=NK+J 
JI=JP+J 
H0LD=-A< JK> 

A<JK)=A(.JI  ) 

RO 70  A(JI)=H0LD 
C 


1£6 


I 
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C  DIVIDE  COLUMN  BY  MINUS  PIV0T  (VALUE  OF  PIVOT  ELEMEN 

ST  IS 

C  CONTAINED  IN  RIGA) 

C 

BOBO  I F (B I GA  >81 00.8090 .81 00 
8090  0=0.0 
RETURN 

8100  0081201*1 .N 

IE(I“K)B1 10  #8120  *8110 
PflO  IK*NKtI 

A< IK)«A< IK)/< -BI GA) 

8120  CONTINUE 
C 

C  REDUCE  MATRIX 

C 

DOB | SOI » 1 «N 
I K*NK+ I 
H0LD=A( IK) 

IJ»I-N 
008  I  50 J= 1 »N 
I J=I J+N 

IF(I-K>8130. 8150. 8130 
8130  IF(J-K>8140. 8150. 8140 
8140  KJ=IJ-I*K 

A ( I J) =H0LD*A<  K J) *A< I  J) 

8150  CONTINUE 
C 

C  DIVIDE  ROW  BY  PIVOT 

C 

KJ=K-N 

DOB  |  70.1=  1  »N 

K J=K J+N 

IF (J-K) 8 160.8 170. 8 160 
8160  A(K.J>=A(KJ>/BIGA 
8170  CONTINUE 
C 

C  PRODUCT  OF  PIVOTS 

C 

DaD+RIRA 

C 

C  REPLACE  PIVOT  BY  RECIPROCAL 

C 

A  <  KK)  =  1  .O/niGA 
8180  CONTINUE 
C 

C  FINAL  ROW  AND  COLUMN  INTERCHANGE 

C 


1Z7 


i 


n  n 
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K*N 

8190  K*CK-1> 

IFCK)8260#P260#R200 
8200  I *LCK) 

I F (  I  -K)  82.30  # 8230  # 82  1 0 
8210  JO*N*(K-l ) 

JR*N*CI»1) 

D0822OJ*! *N 
JK*JQ* J 
HOI  0*A(JK) 

JI  *JRO 
A( JK)*-ACJI ) 

8220  A< JI )=H0LI) 

8f*:;o  j*mck> 

I FCJ-K>B 190 #8190*8240 
8240  KI»K-N 

0082501*1 #N 
Kl *Kl ♦N 
W0LD*ACKI) 

JI*KI-K*J 
ACKI ) *"A( JI ) 

8250  A< JI ) *M0LD 
O0T08I9O 
8260  RETURN 
ENO 


SUBROUTINE'"  ,^APH<C  #V#CST#NAME  #NP  #NC# IBUF#NCHAR) 
DIMENSION'  <202) #C(202) # VC 202) # IBUFC 1006) #NAME< 15)#CSTC 

S22) #D1 < 10) #D 

42(10) #YSM(502) #X SMC 502) 


NC2*NC*2 
D09OOOI *1 #NC 
D1CI)*CSTC2*I-1> 

9000  D2CI)=CSTC2*I)*l00./2 
D09O1OI-1 #NP 
YCI)=100.*VALCI) 

9010  CCI )*100.*0CI  ) 

XSMC 1 ) *VC 1 ) 

DELX=2*ABSCVC 1 ))/500» 

D09O2O I  S=2*  ?.0Q 
9020  XSMC I S) *XaMC I S ”1 ) ♦DELX 
I)09O3OIS=1  #500 
YSMCIS)=0. 

D09O3O JI  =2#NC2#2  . 

9030  YSMC I S) *YSMC I S) +RNTZ ( XSMC IS)#CSTCJI)#C5TCJI-1)#CST. NC2 

*♦1 ) )*100. 


1<8 
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C ALLPL0T C  P • * 4 • » -3  ) 

CAU.SCALFCV*A.O*NP*l  ) 

CALLSCALFCC*3.0*NP#-I  ) 

YCNP+1 )*CCNP»1 ) 

?.iu.«ls<O--1>.^.l7HVFL0C!TYCC«/SEC».-l7.B..O..V<NPM 

14.A..W.  .YCNPM>.Y<NP 

*♦2)) 

CALLPUflT<0.#-l .25*3) 

CALLP»-0TCO.*O.*2) 

CALLPL0TCO.*3»*3) 

CALLPL0TCO.*4.25*2) 

CALLPL0TC8.#4.2S#2) 

PALLLABF*  Co!ilv7^7.5#3.75#NAMEC1).PO.O. 2*2*0. 1#0.#0) 
J*ALULABF.LC0*5#3  *75*7  *S#3  •  75»NAME<6) *NCHAR#0  *P* I  *0*  1  »0  • 
S#0) 

0090401*1 *NP 

CCI)*CC1)-CSTCNC2»2)*100. 

9040  YC!)«YCI)-CSTCNC2+2)*I00. 

,0,0  CALLSYMB.L^ CV< 1 7 *VCNP*I >)/V<NP*2>#<C(l7 “YCNP+I >  7/YCNP* 

#2) #0 .04# 1*0. 

AO*-l) 

XSMC50I )*VCNP+I ) 

XSMCS02)=VCNP*2) 
ysmcsoi >»ycnp+i ) 

YSMCSOP^YCNP+P) 

CALLFLINFCXSM#YSM#500# 1 #0*0) 

ALNMAX*DI Cl ) 

ALNMIN*D1  C  1  ) 

D09OAO I =1 *NC 

IF  CD  I  CD  .LT.ALNMIN)ALNMIN*DI  CD 
on  ao  IFCDICD.GT  .Al.NMAX)  ALNMAX*DI  Cl) 

CAll.PI.0TC CAI.NMIN”VCN'r'*l  ))/VCNP+2)*0**3) 

CAI.l.P».0TCCALNMAX-VCNP  +  D)/VCNP+2)*O.*2) 

~^(,,w,,,W»P.8,<«.,Y.HP.,,>mNP. 

9n70*RAU-PWT<<DlCI>-V(NP*1»»/V<NP*P>.0..P) 

CALI  .PI.0TC  12.*“4.*“3) 

RETURN 

FND 


C 


FUNCTI0NRNT7CX*AMP*P0S* GAMMA) 


A  ■'Hi 

jl. 
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OALLPI-0T<2»*4.  »"3> 
CAU-SCALPCV*««0#NP#1> 
eALLSCALE<C#3.0#NP#-l  > 


YCNP+1 >*CCNPM > 
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CHAPTER  I 

INTRODUCTION 

Amorphous  semiconductors  have  been  known  to  the  scientific 

community  for  many  years.  One  system  of  these  materials,  chalcogenide 

glasses,  which  consist  of  a  combination  of  a  metal  with  the  heavier 

elements  of  the  oxygen  group  of  the  periodic  table,  were  thought  to 

behave  similarly  to  intrinsic  crystalline  semiconductors.  Kolimiets 

et  al.1  in  1955,  investigated  several  amorphous  materials  made  from 

various  combinations  of  antimony,  arsenic,  selenium,  sulfur,  tellurium 

and  thallium.  They  reported  high  room  temperature  electrical  conduct!- 

2 

vities.  Some  theoretical  theories  were  advanced  and  some  additional 

3 

experimental  work  was  done,  but  in  the  most  part,  these  glasses 
remained  a  laboratory  curiosity  for  the  ten  years  after  Kolimiets'  in¬ 
vestigations.  With  the  discovery^  in  the  early  sixties  of  an  electrical 
switching  phenomenon  in  a  chalcogenlde  glass,  much  Interest  has  been 
shown  to  the  development  of  these  and  associated  glasses  for  use  as 
solid  state  devices.  Characterization  of  the  thermal,  electrical, 
structural,  and  optical  properties  has  been  undertaken  at  numerous  in¬ 
dustrial  and  university  laboratories. 

The  majority  of  the  work  done  in  these  systems  has  been  an 
attempt  to  create  stable  and  reproducible  devices  which  exhibit  switch¬ 
ing  and/or  memory  characteristics.  Attempts  to  explain  these  unusual 
phenomenon  observed  in  terms  of  conventional  band  and  transport 


1 


2 


B  t  .  7,8,9,10 

theories5,6  have  not  been  universally  accepted. 

M»ny  of  the  assumptions  of  the  bend  theory  presuppose  a  homogen¬ 
eous  material.  An  important  property  of  glasses  that  seems  to  be 
ignored  by  proponents  of  the  above  theories  is  the  possibility  of  forma¬ 
tion  of  a  metastable  immiscibllity  in  most  glass  systems.  Haller,  et 
aj.U  have  reported  an  explanation  of  this  phase  separation  as  it  occurs 
in  a  sodium  boro.lllc.te  system.  Consideration  of  this  phenomenon 
could  lead  to  explanations  of  the  inconsistencies  in  these  theories.  It 
uould  also  support  the  filamentary  conduction  hypothesis,10  which  is 
based  on  an  electro-thermal  process.  Pear:  ,n12  ha.  reported  phase  ^sepa¬ 
ration  in  several  switching  and  memory  glasses  and  Stocker,  et  al.  have 


done  theoretical  calculations  supporting  the  thermal  theories. 

This  author  chose  to  investigate  one  of  the  chalcogenld.  systems 
which  exhibit,  memory  and  switching  characteristics  to  determine  if 
structural  lnhomogeneltles  are  Indeed  present  and  if  they  have  any 
effect  on  other  properties.  The  system  arsenic  trl.elenide-ar.enic 
tritelluride  was  chosen  since  no  work  ha.  been  done  on  the  structural 
characterisation  of  these  glasses.  It  was  decided  to  study  the  electri 
cal  properties  of  the  system  for  evidence  of  possible  correlations  with 


microstructure. 

Thus  the  purpose  of  this  investigation  was  to  measure  the  A.C. 
and  D.C.  electrical  propertlea  and  examine  the  mlcroatructure  of  glasses 
in  the  As2Te,-A.2Se3  system  to  determine  if  any  correlation  exists  be¬ 
tween  these  properties.  This  paper  presents  a  report  of  this  investiga- 
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CHAPTER  II 

EXPERIMENTAL  PROCEDURE 

Six  glasses  were  studied  from  the  arsenic-selenium-tellurium 
system.  The  compositions  were  xAs2Te3-(100-x)As2Se3  with  x  equal  to 
80,  70,  60,  50,  40,  and  0.  All  compositions  lie  on  the  As2Te3-As2Se3 
join  of  the  ternary  phase  diagram,  Figure  1.  Throughout  this  paper, 
compositions  will  be  denoted  as  a  ratio  of  mole  per  cent  As2Te3  to 
mole  per  cent  As^e^.  Thus  80/20  Indicates  a  composition  consisting 
of  80  mole  per  cent  As2Te3  and  20  mole  per  cent  As^e^. 

The  glasses  were  made  by  sealing  reagent  grade  materials  of  the 
proper  proportions  In  evacuated  Vycor  ampules  and  heating  for  one  hour 
at  800°C.  The  melting  was  done  In  a  rocking  furnace  specially  designed 
to  assure  complete  mixing  of  the  components.  The  furnace  was  a  tube 
furnace  with  a  sample  ampule  holder  mounted  such  that  It  would  pivot 
under  the  action  of  a  motor  driven  cam.  This  caused  the  liquid  sample 
to  be  rocked  back  and  forth  In  the  ampule  as  one  end  was  alternately 
raised  and  lowered  with  respect  to  the  other  end.  The  melts,  still  in 
the  Vycor  capsules,  were  quenched  In  room  temperature  water.  All  sam¬ 
ples  for  the  studies  were  made  from  these  bulk  glasses. 

A  prellmlna  y  Investigation  using  high  purity  (at  least  99.999%) 
materials  handled  under  an  Inert  gas  atmosphere  and  sealed  under  vacuum 
In  fused  silica  ampules  was  conducted.  The  specific  materials  used  are 
listed  In  Table  1.  These  samples  were  given  the  same  melting  treatment 


3 


ternary  phase  diagram. 


5 


TABLE  1 


HIGH  PURITY  MATERIALS 


Te:  American  Smelting  and  Refining  Co.,  Special 

High-Purity,  Semiconductor  Grade  (99.999%)  Tellurium 


Se:  American  Smelting  and  Refining  Co.,  Special  High 

Purity  (99.999+%)  Selenium 

As:  United  Mineral  and  Chemical  Corp.,  Arsenic  Powder, 

99.9999%  Pure 


as  described  above  and  showed  no  variation  in  properties  from  samples 
made  with  the  reagent  grade  materials.  It  is  also  well  established  in 
the  literature1^*1^'1^  that  amorphous  semiconductors  in  general  are  not 
sensitive  to  impurities.  However,  Mackenzie16  did  note  some  marked 
effect  by  oxygen  impurities  to  heat  treated  selenium  samples. 

Studies  conducted  consisted  of  measurement  of  the  electrical 
properties,  replica  electron  microscopy  (E.M.),  Guinier-de  Wolff  x-ray 
analysis,  differential  thermal  analysis  (DTA)  and  dilatometric  analysis. 

Samples  of  the  bulk  glasses  were  powdered  and  studied  in  a 
Guinier-de  Wolff  powder  camera  to  verify  that  the  materials  were  indeed 
non-crystalline.  Crystals  of  0.1  weight  per  cent  can  be  detected  by 
this  method. 

Electrical  specimens  in  platelet  form  were  prepared  by  briefly 
softening  a  sufficient  quantity  of  the  bulk  glass  on  a  graphite  plate 
and  pressing  with  a  second  graphite  plate.  All  samples  had  an  area  of 
approximately  one  square  centimeter  and  a  thickness  of  about  2  millime¬ 
ters.  All  samples  had  similar  thermal  histories.  No  samples  were  sub¬ 
jected  to  further  thermal  treatment.  Conductivity  measurements  were 
made  through  silver  conductive  paste,  Walso  #36-1  Silver  Print,  in  a 
guard  ring  configuration.  The  applied  fields  were  approximately  five 
volts  per  centimeter.  This  is  much  less  than  that  required  to  switch 
these  glasses.  Electrical  measurements  were  divided  into  D.C.  and  A.C. 
segments,  with  the  A.C.  segment  being  further  split  into  audio  and  radio 
frequency  regions.  The  same  sample  was  used  for  both  A.C.  and  D.C. 
measurements.  The  specially  constructed  sample  holder  consisted  of  a 
shielded  configuration  (Figure  2),  with  an  Iron-Constantan  thermocouple 
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Sample  holder,  without  shield 
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mounted  adjacent  to  the  sample.  Electrical  measurements  on  the  As^Se^ 
composition  were  considered  unreliable  and  electrical  data  on  this  com¬ 
position  was  generally  not  included  in  this  paper. 

The  D.C.  conductivity  measurements  were  taken  over  a  range  of 
temperatures  from  room  temperature  to  about  115°C.  Sample  temperature 
control  was  facilitated  by  use  of  a  heating  tape  and  variable  transfor¬ 
mer.  A  Keith ley  high  speed  picoammeter  was  used  to  measure  current  in 
conjunction  with  a  Hewlett  Packard  constant  voltage  source.  Current 
readings  were  taken  at  approximately  10°C  intervals.  The  heating  rate 
was  determined  by  calibration  of  the  variable  transformer  and  heating 
tape  system.  The  rate  chosen  was  as  close  to  equilibrium  heating  as 
practical  without  causing  annealing  of  the  sample  at  the  higher  temper¬ 
atures  and  to  minimize  differences  In  sample/thermocouple  temperature. 

The  A.C.  conductivity  and  capacitance  measurements  were  taken  at 
room  temperature.  The  audio  frequency  measurements,  from  200  Hz  to  50 
kHz,  utilized  a  Hewlett  Packard  615  B  oscillator  coupled  with  a  Wayne 
Kerr  B221  audio  frequency  bridge,  which  had  an  internal  null  detector. 

A  Wayne  Kerr  B601  radio  frequency  bridge  and  General  Radio  external  null 
detector  were  used  for  the  radio  frequency  measurements.  Due  to  limita¬ 
tions  of  the  null  detector,  an  oscillator  and  a  diode  radio  frequency 
mixer  were  used  to  mix  a  beat  frequency  with  the  output  of  the  radio 
bridge  for  the  measurements  using  frequencies  greater  than  100  kHz. 

Dilatometric  studies  were  undertaken  to  determine  the  softening 
point  of  the  glasses.  This  was  to  establish  a  safe  temperature  limit 
for  the  D.C.  measurements  and  an  annealing  temperature  for  later  exper¬ 
iments. 
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Differential  thermal  analyeea  (DTA)  were  conducted  to  determine 
the  mm  of  devitrification  and  the  presence  of  glass  traneltione.  The 
DTA  data  was  taken  using  powdered  .ample*,  with  a  modified  Fisher  Dif¬ 
ferential  Thermal, ter.  Model  260P.  Heating  rates  of  5,  10.  and  25°C  per 

minute  were  used  on  all  compositions. 

Electron  microscopic  studies  were  made  of  fresh  fracture  surfaces 

and  etched  fracture  surface,  using  two  stage,  platinum-palladium  shad¬ 
ed.  carbon  replicas.  The  etchant  was  a  dilute  solution  of  potassium 
hydroxide  and  water.  The  microscope  was  a  Philip.  E.M.  300. 
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CHAPTER  III 

EXPERIMENTAL  RESULTS 

Guinier-de  Wolff  x-ray  analysis  showed  no  evidence  of  crystal¬ 
lization  In  any  of  the  Initially  prepared  glasses.  Results  of  the  DTA 
studies  are  summarized  In  Table  2.  Two  clearly  defined  exothermic 
reactions  were  noted  for  the  80/20  and  70/30  glasses.  The  60/40  and 
50/50  glasses  each  exhibited  only  one  reaction.  The  glasses  with 
higher  Se/Te  ratios  showed  very  broad  and  diffuse  peaks  which  made  the 
accurate  determination  of  reaction  temperatures  difficult.  No  reaction 
could  be  determined  for  the  40/60  glass  nor  for  the  ASgSe^  glass.  High 
selenium  glasses  are  generally  difficult  to  crystallize,  and  any  reac¬ 
tion  probably  would  proceed  too  slowly  to  be  noticeable  In  our  experi¬ 
ments.  The  results  Indicate  that  the  crystallization  temperature  of 
the  glasses  Increases  as  the  Se/Te  ratio  Increases.  Our  equipment  was 
too  Insensitive  to  show  any  reactions  which  Involve  small  heats  of 
reaction.  Thus  reactions  such  as  glass  transitions  were  not  detectable 
In  our  studies.  The  experiments  were  not  continued  to  melting.  Devi- 
tlficatlon  In  all  samples  except  arsenic  triselinide  was  confirmed  by 
x-ray  analysis.  The  author  was  unable  to  Identify  the  crystalline 
phases . 

The  softening  points  determined  by  the  dilatometric  studies  are 
listed  In  Table  3.  The  softening  point,  or  softening  temperature,  as 
determined  by  the  dilatometric  studies  Is  that  maximum  temperature  at 
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DTA  DATA 
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TABLE  3 


DILATOMETER  DATA: 

SOFTENING  POINTS 

Softening  Temperature 

As2Te3/As2Se3 

°C  ±  2° 

80/20 

100°C 

70/30 

110°C 

60/40 

125°C 

50/50 

135°C 

40/60 

145°C 

0/100 

185°C 

153 
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which  the  glass  can  hold  Its  shape  and  resist  deformation  due  to  small 
applied  forces.  It  can  be  viewed  as  that  temperature,  or  range  of  temp¬ 
eratures,  at  which  the  viscosity  of  the  glass  is  lowered  to  the  point 
that  the  material  is  no  longer  rigid,  but  assumes  more  of  the  character 
of  a  liquid.  The  term  actually  has  varied  definitions  and  values,  de¬ 
pending  on  the  method  of  measurement,  heating  rates,  and  magnitude  of 
force.  The  Orton  dilatometer  measures  the  linear  expansion  or  con¬ 
traction  of  the  sample  while  it  is  being  heated.  A  slight  compressive 
pressure  Is  maintained  on  the  sample  by  a  small  spring.  As  can  be  seen 
from  the  dilatometer  curve  of  Figure  3,  there  Is  a  point,  actually  a  3 
to  5°  range,  where  the  thermal  expansion  levels  off  and  drops  suddenly 
to  zero.  The  minimum  temperature  of  this  range  we  defined  as  the  soft¬ 
ening  point.  The  data  Indicates  that  the  softening  temperature  rises 
with  Increased  selenium  content. 

Electrical 

Data  from  the  D.C.  experiments  are  plotted  in  Figures  4  through 
9  as  log  conductivity  versus  inverse  absolute  temperature.  These  plots 
are  combined  in  Figure  10.  The  samples  behave  as  expected  of  electron¬ 
ically  conducting  semiconductors,  with  the  conductivity  increasing 
according  to 


a  -  0Q  exp  (-  AE/kT)  . 


A.C.  electrical  measurements  were  interpreted  in  terms  of  the  loss 
tangent,  tan  6  -  G/ojC.  This  function  is  equivalent  to  the  tangent  of 
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the  ratio  of  the  lo.e  factor,  t\  to  the  permittivity,  e". 
pan  del  can  he  considered  to  be  composed  of  two  components,  an  i  «• 
frequency  dependent  component,  and  a  D.C.  component.  Thus 


tan  6 


Total 


-  tan  6ac>  +  tan  6D 


The  ioss  at  low  frequency  Is  higher  and  the  D.C.  component  overwhelms 
the  effect  of  the  A.C.  component.  Thus,  to  examine  the  frequency  de¬ 
pendence,  we  must  subtract  the  D.C.  component,  based  on  the  room  temp¬ 
erature  D.C.  measuranents,  f»  the  tot.!  value  of  the  tan  del  function 
computed  from  the  A.C.  measurements.  Thus 


tan  6A.C.  *  tan  6Total  '  tan  6D.C. 
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A  full  discussion  on  the  tan  del  function  can  be  found  in  van  Beck. 

This  tan  del  function,  with  the  D.C.  component  subtracted,  is 

,  -  frequency  for  each  composition  in  Figures  11  through 

plotted  versus  log  frequency 

15.  These  same  plots  are  cabined  in  Figure  lb.  *11  samples  exhibited 
at  least  one  loss  peak  and  corresponding  dielectric  dispersion.  This 
indicates  that  conditions  for  haxwell-V.gner-Sill.rs  heterogeneous  msecs 

are  present  in  these  glasses. 

The  activation  ener,ie.  and  the  corresponding  pre-exponential 
term,  °o,  were  calculated  from  the  equation 
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Figure  11.  Tan  del  versus  log  frequency,  80/20  glass. 
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Figure  16.  Tan  del  versus  log  frequency. 
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using  a  non-linear  least  squares  fit.  This  data  is  shown  in  Table  4, 
along  with  the  data  of  Roilos.19  Roilos'  data  was  adjusted  by  a  factor 
of  one  half  due  to  a  difference  in  the  equations  used  by  the  two  exper¬ 
imenters. 

It  was  noted  that  surface  conductivity  is  significant  in  these 
glasses.  Measurements  with  unguarded  samples  resulted  in  conductivities 
approximately  double  those  of  the  guarded  samples.  A  guard  ring  config¬ 
uration  as  was  used  was  necessary. 

The  computations  and  plotting  of  the  electrical  results  were  by 

computer.  The  programs  are  included  as  an  appendix. 

E.M. 

EM  studies  clearly  indicate  that  liquid-liquid  phase  separation 
occurs  in  all  the  glasses  studied,  including  arsenic  triselinide.  In 
all  the  compositions  the  droplike  structure  characteristic  of  nuclca- 
tion  and  growth  type  of  phase  separation  is  apparent,  while  some  also 
show  a  structure  characteristic  of  a  splnodal  decomposition  process. 

These  processes  are  explained  in  the  following  chapter.  Some  of  the 
droplets  in  the  glasses  have  "tails,”  as  shown  by  the  arrow  in  Figure  17. 

These  tails  are  formed  as  the  fracture  front  passes  around  a  drop  and 
changes  its  plane  of  fracture.  They  are  proof  that  the  droplets  arc  a 
part  of  the  structure  of  the  fracture  surface,  and  not  simply  artifacts 
from  the  replicating  procedure.  The  conchoidal  mode  of  fracture  and  the 
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i2Te3/As2Se3 

TABLE  4 

ACTIVATION  ENERGY 

Activation  Energy 
(AE)  (e.v.) 

a  , 

(x°10J)  (fi  •  cm) 

80/20 

0.43 

0.31 

*75/25 

0.49 

*71/29 

0.51 

70/30 

0.52 

3.14 

*66/33 

0.515 

60/40 

0.55 

7.3 

*60/40 

0.53 

50/50 

0.60 

3.1 

* 50/5 r> 

0.575 

40/60 

0.62 

1.25 

*33/66 

0.64 

I 


Figure  17.  Electron  micrograph  showing  droplet  "tall." 


Figure  18.  ASgSej  fracture  surface. 
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x-ray  results  Indicate  that  both  phases  are  non-crystalline. 

The  electron  photomicrographs  of  arsenic  trlsellnide.  Figures 
18,  19,  and  20  show  the  droplike  structures.  The  drops  are  about  one 
micron  in  diameter  and  are  well  formed.  They  occupy  approximately  a  ten 
per  cent  volimie  fraction.  A  splnodal  structure  is  not  apparent  in  this 
glass.  Note  also,  in  Figure  20,  that  the  fracture  front  has  passed 
through  the  drops.  This  is  indicated  by  the  fracture  pattern  in  the 
drops  and  the  lack  of  relief  in  the  surface.  This  implies  that  the  drop 
and  the  matrix  have  similar  mechanical  properties  and  are  probably  of 
similar  composition. 

Micrographs  of  the  40/60  glass.  Figures  21,  22  and  23  show  large 
sized  drops,  primarily  about  three  microns  in  diameter,  but  with  some 
smaller  drops  apparent.  Also  apparent  is  an  interesting  structure  with¬ 
in  the  drops  not  seen  m  the  other  glasses  which  may  be  indicative  of  a 
secondary  segregation  occurring  in  the  primary  separated  phase.  20  The 
volume  fraction  of  the  large  drops  is  less  than  hen  per  cent.  No 
attempt  was  made  to  determine  the  amount  of  secondary  segregation.  The 
fracture  "tails"  are  evident  in  this  composition. 

Ihe  50/50  glass  shows  the  droplet  structure  and  also  the  aplnodal 
structure,  see  Figure  24.  The  drop  volume  fraction  is  somewhat  less 
than  ten  per  cent,  but  the  second  phase,  due  to  the  splnodal  decomposi¬ 
tion,  is  much  larger. 

The  60/40  composition  exhibits  1  to  1%  micron  sized  drops  In  a 
splnodal  matrix.  See  Figure  25.  Voltae  fraction  of  the  drops  is  about 
2  OX. 

icrograpk  of  the  70/30  glass  show  a  large  quantity,  approximately 
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Figure  19.  As2Se3  fracture  surface. 


Figure  21.  40/60  glass  fracture  surface. 


Figure  22.  40/60  glass  fracture  surface. 


Figure  23.  40/60  glass  fracture  surface. 


Figure  24.  50/50  glass,  etched  fracture  surface. 
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40%  volume  fraction.  -all  droplets.  The  Mjority  are  less  then  one 
micron  In  dimeter.  Huny  of  the  droplet,  are  combining  or  are  connnec- 
red  in  Figure.  26  through  28.  the  droplet  morphology  « 
apparent,  and  1.  distinguishable  from  a  spinodal  .tructure. 

The  80/20  ,1...  .hove  a  structure  indicative  of  a  spinodal  decom¬ 
position  process.  The  droplet  morphology  1.  discernahle.  hut  the  struc¬ 
ture  1.  primarily  that  of  tuo  interconnected  continuous  phases.  The 
volume  fraction  of  the  second  phase  1.  about  30%.  Figure.  29  snd  30  are 


of  the  80/20  glass. 
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Figure  27.  70/30  glass,  etched  fracture  surface. 


Figure  28.  70/30  glass,  etched  fracture  surface. 
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CHAPTER  IV 

DISCUSSION 

The  results  of  the  x-ray  atudlea  indicate  that  there  are  no 
crystal  phases  of  greater  than  0.1  weight  per  cent  present.  The  DTA 
studies  verify  that  the  material  is  amorphous,  by  the  occurrence  of  a 
crystallisation  reaction.  These  studies  also  indicate  that  the  glas¬ 
ses  with  higher  selenium  content  are  better  glass  formers,  as  expected. 

A  set  of  sample  DTA  curves  showing  the  compositional  dependence  is 
given  in  Figure  31.  The  reaction  temperatura  observed  for  each  sample 
were  dependent  on  the  heating  rate,  indicating  a  kinetically  limited 
process.  Figure  32  shows  the  heating  rate  dependence  of  the  reaction 
temperatures  in  the  80/20  glass.  The  occurrence  of  double  crystalliza¬ 
tion  peaks  in  the  S0/20  and  70/30  glasses  indicates  that  there  are  two 
crystalline  phases  forming  in  these  compositions.  These  phases  may  both 
crystallize  from  a  single  phase  glass,  but  the  structural  studies  indi¬ 
cate  that  there  are  two  glass  phases  present.  Crystallization  of  each 
phase  would  result  in  a  DTA  peak.  No  double  peaks  are  observed  in  the 
other  glasses  even  though  our  E.M.  studies  show  evidence  of  phase  sepa¬ 
ration.  The  lack  of  the  double  peak  could  be  due  to  the  second  ptiase  be 
ing  present  in  only  small  quantities.  This  is  supported  by  the  trend  to 
ward  larger  amount  of  phase  separation  in  the  higher  tellurium  glasses 

which  was  noted  in  the  E.M.  studies. 

Bagley  and  Bair^1  conducted  differential  scanning  calorimetry  on 
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Figure  31.  Composition  dependence  in  DTA  curves. 


182 


43 


arsenic  trl.ellnlde-.rs.nlc  trltellurlde  glasses  with  Aa^/Aa^e, 
ratio,  of  3/1.  2/1.  end  1/2.  It.,  cryatalllaatlon  ta.per.tura.  they 
deteralned  (...  Teble  2)  agree  well  with  our  data.  Houeyer.  they  noted 
.  .mall  exothemlc  reaction  preceding  crystallisation  uhlch  they  nttrl- 
buted  to  .  glee.  tr.n.ltlon.  They  alao  noted  a  doubly  peaked  endothcr. 
.t  uniting.  They  offers  no  expl.natlon  for  the  double  pe.V..  In 
light  of  the  ph...  aep.r.tion  our  structural  atudle.  ahou  1.  preaont  In 
the.,  glaaaea.  It  1.  reaaonable  to  .ttrlbute  the  double  peak,  to  a  tvo 
phase  melting  occurring.  Lack  of  a  double  peak  at  cry.tallU.tlon  coulu 


be  due  to  simultaneous  cryatalllaatlon  of  both  phaaea. 

.agio,  and  horthov.r22  did  not  detect  phase  operation  by  elec¬ 
tron  microscopy  In  the  glaa...  »««•  >>1  “d1**  B*lr'  “°W'V*r' 
primarily  u.ed  transmission  mlcro.copy  on  thin  film.,  with  some  replica 

work  on  unetched.  fracture  surfaces.  Phaae  separation  1.  difficult  to 
detect  in  ap.clm.ns  that  have  not  been  etched.  If  the  electron  densi¬ 
ties  of  the  separated  phases  are  similar,  then  It  is  difficult 
Impossible  to  detect  ph...  separation  In  transmission.  Also,  a  thin 
film  of  an  amorphous  material  does  not  necessarily  have  the  same  struc¬ 


ture  88  the  bulk  materiel. 

The  electron  microscopic  studies  Indicate  the  presence  of  separa¬ 
ted  phase,  in  all  the  compositions  studied.  The  DTA  and  x-ray  data  shows 
both  (or  all)  of  the  phases  to  be  amorphous.  Since  the  phase  separated 
glasses  include  both  binary  side,  of  the  ternary  phase  diagram,  It  Is 
likely  that  the  phase  separation,  and  therefore  an  Immlsclblllty  gap. 
extend,  across  the  entire  arsenic  trltellurlde-arsenlc  trlsellnlde  Join. 

The  morphology  of  the  phases  observed  In  the  glasses  Indicate 
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that  two  phase  separation  processes  are  occurring:  a  nucleation  and 
growth  type  of  separation  and  a  splnodal  decomposition  process.  Nuclc- 
atlon  and  growth  occurs  by  the  classically  iccepted  process  of  growth 
of  a  heterogeneously  nucleated  particle.  Formation  of  well  rounded 
spheriods  or  droplets  is  indicative  of  thin  type  of  diffusion  controlled 
growth  occurring  in  an  amorphous  material.  The  surface  tension  of  the 
glassy-liquid  controls  the  shape  of  the  drop.  This  structure  is  sia.^ar 

tc  a  suspension  of  water  and  oil. 

Splnodal  decomposition  occurs  spontaneously  or  homogeneously 

throughout  the  parent  phase.  It  has  more  the  appearance  of  a  natural 

sponge  parent  phase,  with  the  interconnected  voids  of  the  sponge  filled 

with  the  second  phase.  Both  new  phases,  of  course,  have  a  different 

composition  than  the  original  phase.  A  sismnary  of  the  thermodynamic 

13 

description  of  these  processes  by  Cahn  and  Charles  is  given  here. 

A  necessary  condition  for  stability  of  a  phase  with  respect  to 
infinitesimal  composition  or  density  fluctuations  in  that  the 
chemical  potential  of  each  component  increase  with  increasing 
density  of  that  component.  This  condition,  as  will  be  seen,  is 
not  in  conflict  with  the  existence  of  a  more  stable  phase,  and 
thus  metastability  can  and  does  occur  commonly.  In  some  systems 
there  exists  a  temperature,  which  is  a  function  of  composition, 
where  this  necessary  cri  erion  for  stability  breaks  down.  This 
limit  to  the  metastability  has  been  called  the  splnodal.  Beyond 
the  splnodal  a  single  phase  is  unstable  with  respec*  to  infini¬ 
tesimal  composition  fluctuation  and  begins  to  separate  into  two 
related  phases  differing  only  in  composition.  This  occurs  spon¬ 
taneously,  without  the  need  for  nucleation,  since  the  phase  is  no 
longer  metastable  but  unstable.  The  related  phases  that  form  may 
be  the  equilibriun  phases,  or  they  may  in  turn  be  only  metastable 
with  respect  to  still  more  stable  phases  and  subject  to  further 
application  of  stability  criteria. 

In  a  binary  system,  Gibb’s  condition  for  metastability  reduces  to 
32F/3c2>0  where  F  is  the  molar  free  energy  and  c  is  the  mole  frac¬ 
tion.  The  splnodal  is  to  be  expected  in  every  system  that  exhibits 
stable  or  metastable  equilibrium  between  two  related  phases.  It  is 


185 


45 


t  at  coHonly  to  be  expected  in  systems  bavin*  s  coneolute  teaper- 
ature.  Consider  a  binary  system  with  an  upper  t 

Hmmm  called  critical)  teaperature,  T  (Figure  33).  At  a  constant 
tLperature  above  the  coneolute  teaperature,  the  molar  ff**en*^y 
F  as  a  function  of  coaposition  aust  curve  upward  ey*ry^J*re  •  .V\ J 
1_  33b  (otherwise  separation  into  two  phases  would  reduce  free 
A*  thete.per.turc  i.  dropped  the  cur.,  n...  by  « 
proportional  to  the  molar  entropy. 

-  II  -S 
3T 

Because  the  entropy  is  expected  to  be  anal lest  for  thetwopure 

Iiauida  (an‘  for  coapounds)  and  greatest  near  the  aiddle  of  the 
Hqutd.  t.nd-for^coupo  ^  ^  curvt  th.  t«p.r.t»rc 

SSSrss  lAVZSZJ*  ws  irasz 

these  limits  a  single  phase  is  still  stable. 

In  Figure  33b  the  inflexions  in  F  occur  at  c  and  d  ***  "J1"' 

odal  compositions  for  the  temperature  Ti,  iyingsomewherebe tween 
nn .  h  Tf  the  free  energy  function  exists  and  has  no  discontinue i 

insLe  tJ.  inflexions  must  exist  at  all  temperatures  below  the 
in  siepe,  tn  Therefore,  the  unstable  region,  be- 

srssaSTSTSEi  » -  — « che  ““°iute 

temperature. 

The  metastable  region  exists  between  the  spinodal  and the  1 phase 
k  ,,j.rv  If  the  single  phase  is  cooled  from  above  the  phase  boun 

n  £u  be  mct.Bt.bl.  «d  mould  only  d.com- 

poLTtoth.  ti  Ph....  If  th.  ..cond  ph...  — 

attempt  by  the  system  to  separate  into  regions  differing  only  slig 
lS  inPcom£osition  will  raise  the  system  free  energy.  If. howeve  , 
the  single  phase  is  brought  within  the  unstable  region  <b*twee"  c 
th.  .u  ovetem  can  continuously  lower  its  free  energy  by  contin- 
SL'LMJSl “Sir-  until  it  cache,  th.  - 

ST  %£%£ ’tnilSS.S:!  compo.ltlon 
the  P  tJ  Becond  liquid  phase  must  form  by  nucleation. 

ablc'an^dlatlnctly^aeparated^partlcleB  nf'thT n.»  ph...  ublch  In 

time  grow  in  size  and  may  coalesce. 

Uh»n  A  sinale  phase  liquid  is  cooled  through  the  metastable  region 
rapidly  enough  to  prevent  appreciable  phase  separation,  t  e  mo  e  o 
^hase  separation  changes  radically  after  it  crosses  the  spinodal  into 
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.  Dori-irlea  which  had  nucleated  and  were  grow- 
the  unstable  region.  ****£”  rpntres  from  which  sinusoidal 
lng,  no.  .top  growing  and  .X^wpletely  lnter- 

eonposltlon  variation,  opr  an"  initially  homogeneous  single 

penetrating  pha.ee  .re  fern*  *  connectivity.  The  geometry 

phase.  Eech  phaM  ha.  a  high  degree  of  connectivity^^  ^  ^ 

of  thl.  tw-ph*.e.truct.re  initially  ^  phM#>  changes  ulth  time. 
MlTTuaU  S£T‘f  rt.  factors  which  differentiate  between  the 
mechanisms. 

the  variation  in  the  wage,  of  pha.e  separation  In  the  .ample. 
...died  make,  it  difficult  to  accurately  correlate  the  *~t  of  phase 
separation  with  compo.ltion.  Generally  the  glaeses  with  a  higher  sel- 
.ni*  content  tend  to  be  more  .table  and  exhibit  les.  pha.«  .eparation. 
Thl.  1.  significant  In  light  of  the  concurrent  .witching  studies  being 
don.  by  Sander.24  In  this  laboratory.  Sander.'  result.  indicate  that 
rh.  higher  ..lenl-  gl—  have  higher  breakdown  voltue.  and  are  more 
■table  switches  than  the  high  tellurl®  «U-~  This  indicate,  a  cor¬ 
relation  between  witching  and  the  structure  due  to  pha.«  «ps-atlon. 
Careful  electron  microscopic  studies  of  he.t  treated  gl...e.  are  needed 
to  determine  the  extent  of  the  pha.e  separation  and  the  compo.ltion.  of 
the  separated  phases.  Further  witching  studle.  of  the.,  heat  treat,  :, 

glasses  would  allow  any  correlation  to  be  verified. 

The  results  of  the  electrical  studle.  are  cond«.lve.  The  con¬ 
ductivity-temperature  behavior  is  typical  of  an  intrinsic  semiconductor. 
The  D.c.  conductivity  data  generally  agree  well  with  the  results  of 
Hollos.25  slight  variation,  are  probably  due  to  thermal  history 
differences.26  The  D.C.  data  as  plotted  In  Figure  10  Is  repeated  in 
Figure  3*  With  the  data  of  Sollos.  The  activation  energies  are  compared 
in  Table  4.  The  room  temperature  conductivity  also  agree,  well  with 
that  of  Mott  and  Dwla.26  Thl.  data  1.  graphed  In  Figure  35.  As  noted 
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Figure  34.  Log  conductivity  versus  temperature  with  Roilos 
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I*  ?*VU  4*  |M«  I*  «  In  |ft« 

)| 

«  m*im  at  U*  40/40  cottmIiIj*.  Il  apt**'*4  t»  Hall  aad  &•*»* 
tut  tM  wwnui  t#f»  -ruiaiiwd  «i«iunr  c«i*um."  tn* 

vartatU*  In  IN*  O.C.  *l*ctrlcal  prcr*rt !•*  with  exposition  !•  prob¬ 
ably  d« *  to  IN*  lAt«ratltl«l  •ubatitutloii  of  a*l*oti*  by  th*  heavier 
tallwriiai  atom.  Structural  rearrangement*  may  *l*o  play  a  tale,  but 
It  I*  fall  that  tc  muU  be  a  minor  rol*  In  thla  can*. 

Iha  dleparaion*  In  tha  tan  dal  varaua  frequency  curvaa  Indicate 
that  Haxwell-Vagner-Slllara  leterogeneoua  loaa  condition*  are  present. 
Ihaee  loaaaa  result  froa  tha  presence  of  structural  lnhonogeneitlcs, 
such  as  crystallites  or  other  large  particles  with  distinct  boundaries. 

There  can  be  son  doubt  that  the  dispersions  in  the  curves  arc 
valid.  By  varying  the  value  of  the  D.C.  component  subtracted  in  the 
tan  del  calculations  and  observing  the  changes  in  the  curves,  one  can 
determine  if  the  loss  peaks  are  real.  If,  as  the  D.C.  value  is  in¬ 
creased,  the  dispersions  straighten  out  and  the  curves  go  smoothly 
through  zero  to  negative  tan  del  values,  then  the  peak  is  not  a  real 
representation  of  the  dielectric  loss.  However,  if  the  peaks  remain 
in  the  curves  as  they  go  negative,  and  the  curves  only  smooth  out  with 
very  large  or  very  small  D.C.  values,  then  the  dispersions  are  real 
and  are  due  to  actual  losses  within  the  system.  All  the  samples  be¬ 
haved  in  the  latter  manner.  An  example  is  given  in  Figure  36.  Thus, 
there  is  a  definite  loss  effect  in  these  glasses.  In  light  of  the  mi¬ 
cros  tructural  studies,  these  losses  are  logically  due  to  the  structural 
inhomogeneities  which  are  a  result  of  the  phase  separation  processes 
occurring  in  the  glasses.  The  variation  in  the  frequency  at  which  the 
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Figure  36.  Effect  of  D.C.  contribution  on  tan  del. 
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differ*!***  mtm  f#r  iW  differ**!  I ***  *******  frt* 

In  i)w  il(r*»if«(i»(t.  There  iff  de*  l#  U#  eerlaile*  l»  1** 

mmm\  of  phew  itWfiiM  1a  tte  |Iih#*. 

Tli*  rfUtl^iihlf  Ntwm  IS?  •irvtture  **1  electrical 
of  these  flatte*  (MMt  yet  So  folly  defined.  It  !•  pe*«lble  iMt  they 
are  related  eolely  through  a  coeyoeltlnnal  dependence.  Thla  Author 
fed,  that  tho  structural  dependence  t«  predominant*  but  Additional 
etudlen  with  heet  treated  glasses  are  needed  to  correlate  the  atructural 
and  corpoattlonal  variables.  The-e  atudlea  are  needed  to  fully  charac¬ 
terize  the  ayatefc.. 
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CO#Ct«IO*i 

All  of  ib#  arsenic  trlMll«U*in«<iic  trltsllurld#  gUMM 
lAijTtjO-iAijItj)  Muflid  fro*  x  •  0  to  x  •  80  art  phase  aaparattd. 
This  indicates  that  thoro  exiats  a  miscibility  gap  acroaa  tho  tut  Ire 
ternary  phaa*  diagram  of  the  arsenlc-selenium-tellurlm  ayata*.  The 
atructure  appeara  to  b«  related  to  and  vary  with  composition,  but 
without  further  atudy  the  relatlonahip  cannot  be  defined.  Analogoua 
phaae  aeparatlon  haa  been  noted  by  other  authora  In  related  chalco- 
genlde  ay a teas. 

The  electrical  properties  of  these  glasses  vary  with  composi¬ 
tion.  D.C.  conductivity  Increases  with  tellurium  content.  Activation 
energy  decreases  with  tellurium  content  and  Oq  has  a  maximum  near  the 
middle  of  the  composition  range.  The  variation  in  the  D.C.  properties 
may  be  due  to  structural  changes,  compositional  changes,  or  a  combina¬ 
tion  of  both.  Structural  lnhomogeneltles  resulting  from  glass-glass 
phase  separation  appear  to  be  responsible  for  the  occurrence  of 
Maxwell-Wagner-Sillars  heterogeneous  losses  In  these  glasses. 
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RK(rfQ?OATin#t  rm  nmnax  study 

Mfctrltil  and  atructural  studies  equivalent  to  those  described 
In  this  thesis  nssd  to  be  conducted  using  best  treated  samples.  The 
use  of  snnealed  electrical  specimens  would  allow  any  compositional 
influence  on  the  electrical  properties  to  be  defined.  Correlations 
between  the  dielectric  losses  and  the  phase  separated  structure  would 
result  from  careful  structural  analysis  of  these  sane  heat  treated 
glasses.  Extreme  care  should  be  exercised  In  the  annealing  process 
and  subsequent  handling  of  these  glasses  as  their  structure  Is  extra¬ 
ordinarily  sensitive  to  variations  in  heat  treatment.  Ideally,  E.M. 
samples  should  be  taken  from  the  actual  electrical  specimens.  Finally 
close  contact  with  the  switching  studies  should  be  maintained,  as  It 
appears  there  are  Important  correlations  In  this  area. 
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u-h-Fa  0F 

JOHN  DAVID  PEARSON 

Recent  intent  in  =  h»  " 

n  the  electrical  twitching  behavior  of  these  materials. 

— — — ■  — T. 

,.QO  r,t  radio  frequencies,  especial 
their  dielectric  properties  at  ra 

u  h.F.  and  microwave  range. 

U'H  „r  census  the  co^itio.1  variation  in  ^ 

dielectric  Parties  of  the  glass  ^  ...  at 

x.0.4,0.5,0.6,0.7,  and  0.8  over  the  frequ  -  - 

room  te,,perature  l25°C) '  .u  mixtures  of  ASj^  and 

The  materials  were  prepar 

After  hp  '  h.  ^  several  hours,  the 

As2Se3  in  a  rocking  furnace.  After 

nwlten  material  was  quenched  in  water.  technlques. 

,.Hvltv  measurements  were  made  by  three 
Complex  permittivity  mea 

0  i  nn-ROO  Mhz,  the  materials  were  b.^ 

F„r  the  frequency  r  coaxial  line.  * 

circular  disks  and  plat  of  thl.  sample,  which 

bridge  was  used  to  determine  the  complex  l(K 

,  ,  dielectric  constant  and  loss  tangent, 

was  used  to  calculate  treasure  null 

500  2  000  M».  a  slotted  coaxial  line  was  used 

‘  wttur  calculation  of  tl»  -H* 

.  ^  VSWR,  which  led  to  the  cuxu 

shirt  and  VO  ,  tbe  ran£0  12,000-18,000  Kte. 

and  the  dielectric  properties.  In  the  ranp. 
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rectangular  shaped  samples  were  used  in  conjunction  with  a  Ku-barri 
slotted  waveguide.  The  dielectric  properties  were  obtained  by  measuring 
VdWR,  wavelength,  and  position  of  the  voltage  minimum  and  using  the 
Von  Hippel  Standing  Wave  Method  to  calculate  the  proper  order  solution 
to  a  complex  transcendental  equation. 

The  average  dielectric  properties  measured  over  ^he  entire 

frequency  range  are  as  follows: 

Composition  e 1  tan  6 

80HA3  're  .-20JAs.jSe,  10.0  0.055 

70&ASjTe  -30jA3,Se,  9.3  0.0H9 

60jASjTe,-i<0jA3jSe,  7.9  0.0^2 

50jAs,Te,-50JAs,Se,  7.2  0.030 

'JOjAs^Te ’-60jAs^  7.0  0.020 

'Iho  dielectric  properties  of  the  couples  were  fourd  to  be 
essentially  frequency  independent  over  the  frequency  range  studied,  with 
si ight  variations  attributed  to  differences  In  experimental  techniques. 
By  inci-easing  the  Tellurites  composition,  it  was  found  that  the 
dielectric  constant  can  be  varied  frar  approximately  7.0  to  10.0.  Hie 
loss  tangent  v.s  also  found  to  increase  with  increasing  Te  content  rvrJ 
varied  between  approximately  0.02  and  0.06.  'fhe  dielectric  properties 
of  this  glass  system  expare  favorably  with  those  of  dielectric 
materials  used  for  microwave  substrates,  such  as  Y1G,  Silicon,  or 
AlJIMag. 
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CHAPTER  I 

IWITODUCTION 

Amorphous  semiconductors  are  disordered  materials  In  which  atoms 

have  local  organization  or  short-range  order,  but  lack  the  periodic 

structure  of  crystals,  for  distances  greater  than  a  few  atomic  radii. 

'Ilie  atomic  positions  in  the  nearest  and  next  nearest  neighbors  are  Uius 

essentially  tlie  same  in  tlx*  amorphous  and  crystalline  states  flj. 

Elements  which  can  he  foicr/xl  Into  amorphous  compounds  cxhib’tlng  tlrc.e 

chiricterlstlcs  Include  germanium  and  silicon  (Croup  IV),  antimony, 

arsenic,  onl  birxwth  (Group  V),  and  selenium,  sulfur,  and  tellurium 

(Group  VI)  (2).  Ihe  present  study  centers  around  examining  some 

properties  of  a  class  of  these  amorphous  materials,  ASgT^-As 

Amorphous  semi e orxluctors  arc  probably  most  well  known  for  their 

switching  and  memory  effects.  In  19^3  S.  R.  Ovshlnsky  began  research  on 

U«M  devices  arri  In  19ol  filed  a  patent  on  an  amorphous  device,  hByiV?- 

30£A3-12SSi-10*Ge,  which  was  Insensitive  to  polarity  In  switching 

operations  (33.  Liter  In  1961,  A.  D.  Pearson,  J.  F.  Dcwald,  W.  R. 

NorUover,  furl  W.  F.  Peck  discovered  switching  and  memory  phcnuHena  In 

As-Tt—1  glass.  Itoweve:  its  reset  properties  were  subject  to  polarity 

IU,5).  D.  L.  Eaton  reported  the  same  phenomena  in  this  glass  In 

(6).  In  1963  Gwitchlrv:  behavior  was  founl  In  glasses  of  the  class 

TlAs(3e,1*)0(  by  Kolorrdets  and  Lebedev  l?J.  Later  W.  R.  Eubanx  founl 
«• 

switching  behavior  In  Cb— 3-1  glass  (83*  Since  .his  time  considerable 
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efforts  have  been  spent  In  studying  the  switching  properties  of  other 
amorphous  materials  [9]. 

A  threshold  switch  is  a  two  terminal  device  which  can  be  in  an 
off  (almost  nonconducting)  state  or  in  ar.  on  (conducting)  state.  When 
Lite  applied  voltage  is  greater  than  the  threshold  value,  the  system 
switches  to  tie  conducting  state  and  produces  a  "negative  resistance" 
effect.  However  ir  the  current  is  reduced  below  its  critical  "holding" 
value,  the  glass  will  switch  back  to  its  original  high  resistance  state. 
Threshold  switching  has  been  seen  to  occur  in  250  picosecond,  with  a  one 
microsecond  delay  before  switching  ClO). 

The  memory  switch  is  similar  to  the  threshold  switch,  but  has  a 
3table  "on"  state,  anJ  tins  Is  a  bistable  device.  When  a  high  voltage 
pulse  Is  applied  to  the  switch  when  "off,"  heating  and  structural 
changes  occur  in  the  device  with  a  partial  ordering  of  its  structure, 
and  the  device  is  permanently  switched  to  tlx?  "on"  sUte.  It  will 
remain  there  even  If  the  current  is  completely  removed.  It  Is  believed 
that  Uic  voltage  pulse  results  in  tlie  growth  of  a  crystal  filament  or 
fvathway  between  anode  and  cathode,  resulting  In  a  low  resistance  path 
for  the  current  ( lj.  The  switch  can  be  restored  to  the  off  state  by 
applying  a  threshold  pulse  of  either  polarity  L5|.  Switching  tlnrs  for 
rr/ro:7  switches  iiave  beer,  observed  to  be  on  the  order  of  a  few 
milliseconds  [lCj. 

In  addition  to  examining  the  switching  properties  of  amorpliou; 
son! conductors ,  nany  researchers  have  studied  various  electrical  and 
physical  properties  of  these  rater  la  Is.  Properties  considered  In  Uie 
last  several  years  include  the  optical  absorption  spectrum  |  11],  the 
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variation  of  conductivity  in  the  audio  and  high  frequency  range,  the 
temperature  dependence  of  the  d.c.  conductivity  [12 ,13 ,1^ ]  and  the 
thermal  switching  properties  [151*  A  general  overview  of  research  in 
the  field  is  given  in  a  recent  National  Materials  Advisory  Board  report 
til]. 

At  Vanderbilt  University  che  switching  properties  of  the  glass 
system  (x)As2Te^(l-x)As2Se^  have  recently  been  studied  [J.6,32].  1i»e 
electrical,  thermal,  and  structural  properties  of  this  system  have  also 
been  extensively  investigated  [173.  The  dielectric  properties  of  the 
glass  in  the  frequency  range  from  d.c.  to  10  Mhz  have  been  studied  as  a 
function  of  microstj’ucture  by  Hill  Ll8J .  The  only  known  research 
ctaracterizlng  the  high  frequejicy  dielectric  properties  of  this  glass 
system  was  published  recently  by  this  writer  and  others  L 193 . 

The  purpose  of  the  research  presented  in  the  present  paper  is  to 
examine  the  U.H.P.  and  microwave  dielectric  properties  of  the  glass 
system  (x)As2*lt ^(l-x)As2Se^  for  x»0. *4, 0.5, 0.6, 0.7,  and  0.8,  and  in 
particular  to  note  their  range  or  variation  as  a  function  of  tellurium 
content . 
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CHAPTER  II 

EXPERIMENTAL 


Introduction 

Dielectric  constant  and  loss  tangent  measurements  were  'made  on 
the  glass  samples  by  three  techniques.  In  the  frequency  range  100-500 
Mhz  a  Thurston  Bridge  was  used  to  determine  the  complex  admittance  of  a 
circular  sample  placed  at  the  end  of  a  shorted  air-filled  coaxial  line. 
Tills  data  was  used  with  the  field  equations  of  a  parallel  plate 
capacitor  to  calculate  the  dielectric  constant  and  loss  tangent  of  its 
dielectric  sanple. 

For  the  range  500-2,000  Mhz  a  slotted  coaxial  line,  terminated  in 
the  sample  and  a  variable  short,  was  used  to  measure  null  shift  and 
VSWR.  This  data  was  substituted  into  the  transmission  line  equations 
for  input  admittance  to  yield  the  sample  admittance,  which  was  used  with 
the  capacitance  equations,  as  above,  to  calculate  the  dielectric 
constant  and  loss  tangent . 

In  the  range  12,000-18,000  Mhz  a  1^-band  slotted  waveguide  was 
used,  terminated  with  the  sample  and  a  variable  short  circuit. 
Measurements  taken  were  wavelength,  position  of  the  voltage  minimum,  and 
VDWR.  A  graphical  technique  given  by  Von  Hippel  was  used  to  solve  the 
equation  for  the  wave  impedance  in  the  sample,  which  led  to  the 
dielectric  constant  and  loss  tangent  of  the  sample. 

The  materials  were  prepared  by  fusing  the  appropriate  mixture  of 
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and  Ao,Se,  in  evacuated  Vycor  anpoults  41  toftPc  lit  a  Krtfclw 
furnace.  After  heating  for  one  leur  U*  eclUft  M>l«rl*l  was  **’ 

water.  The  samples  were  ttien  nlttped  and  polish"!  with  w*U«U  f  '**>*  •/> 
(g*it  carborundum  paper,  into  circular  disks  mi  rwcUivJei  rtr  Ur 
frequency  ranges  100-2,000  f*«  and  12,^00-18,000  Ku  t**p«cllveiy. 

100-500  Mhz 

In  this  frequency  range  the  admittance  of  a  cample  positioned  at 

the  end  of  an  air-filled  coaxial  line  was  measured  by  a  Thurston  l»*ld^% 

Using  this  value  in  the  equations  for  a  parallel  plate  capacitor,  IN? 

dielectric  constant  and  loss  tangent  of  its  dielectric  ample  was 

calculated.  The  parallel  plate  capacitor  arrangement  of  the  sample  1:; 

shown  in  Figure  1  and  its  electrical  equivalent  is  cJwwn  in  Pl^arc  2. 

A  is  the  area  of  the  sample  face  and  T  is  the  sample  thicknec;;. 

T  is  exaggerated  in  Figure  1  for  clarity.  However  a  very  thin  ample 

was  used  to  avoid  fringing  effects.  I  and  Ic  are  the  dissipative  arrJ 

reactive  currents  in  the  sample,  respectively.  Ihe  derivation  of  th" 

permittivity  below  follows  the  treatment  given  by  Harrington  [?0J. 

The  voltage-current  relationship  for  a  capacitor  is 

I  *  I  +1  =  YV  “  (G  +  Ju>C)V  (1) 

B  L 

and  the  field  strength  and  current  density  are,  respectively,  given  by 

E  -  J  <■"> 

and  J  ■  ^  U) 

The  constitutive  relation  for  the  field  in  the  dielectric  is 

J  =  (o  +  ue"  +  Jwe')E  ('•) 
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where  oE  represents  the  conduction  current  component  anu  Coe"  +  jwe')E 
represents  the  displacement  current  components.  Substituting  Equations 
(2)  and  (3)  into  (4), 


I  = 

[(o  +  use"  +  ju>e' )(^)]V 

(5) 

Comparing  (6)  with 

(1)  we 

find 

Y  =  G 

+  jaiC  =  (o  +  we"  +  Jwe' 

)— 

;T 

(6) 

or 

G  =  (a  +  we")^ 

(Y) 

and 

C  =  (e')| 

(8) 

Substituting  C  =  |  in  Equation  (8),  and  dropping  o  from  (7) 
because  it  is  small  for  the  glasses  considered,  and  dividing  both 
equations  by  e^  to  give  relative  values,  we  can  solve  for  the  complex 


permittivity: 


-  *  = 


e"  « 


BT 


e  wA 
o 

GT 

e  oA 
o 


Dividing  Equation  (10)  by  Equation  (9)  we  obtain  the  loss  tangent 


tan  6  = 


e"  _  G 


(9) 

(10) 

(11) 


e'  B 

The  computer  program  COMP  given  in  Appendix  I  uses  the  above  equation. 
with  the  admittance  data  from  the  Thurston  Bridge  to  calculate  the 
relative  dielectric  constant  and  loss  tangent  of  the  sample. 


The  experimental  setup  used  a  G.R.  1602-B  Thurston  Bridge  or 
admittance  meter,  to  obtain  trie  complex  admittance  of  the  sample  which 
was  placed  at  the  end  of  a  siorted  air-filled  coaxial  line.  The  samples 
were  circular  disks,  approximately  0.246"  in  diameter  and  ranged  in 
thickness  between  0.029"  and  0.042".  Figure  3  shows  a  block  diagram  of 
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the  setup  used. 

The  rf  source  was  set  to  the  operating  frequency  and  the  locuJ  rT 
source  was  set  to  a  frequency  30  Mhz  below  this.  A  30  Mhz  amplifier  and 
meter  were  used  to  adjust  the  airline  and  stub  tuner  by  detecting  the 
difference  component  between  the  two  sources,  so  that  the  sample 
admittance  appeared  at  the  bridge's  input.  By  adjusting  the  conductance 
and  susceptance  lever  arms  of  the  bridge  for  a  null  meter  reading,  thie 
sanple  admittance  was  read  directly  in  mlllimhos.  The  frequency  was 
read  by  detecting  a  harmonic  with  a  Transfer  Oscillator  and  a  10  Mhz 
Electronic  Counter.  The  technique  used  in  this  frequency  range  was  not 
repeated  for  higher  frequencies  because  of  the  unavailability  of  two 
stable  rf  sources  covering  frequencies  greater  than  500  Mhz. 

A  90°  coaxial  bend  was  placed  at  the  end  of  the  air-line  with  the 
sample  holder  and  variable  short  placed  vertically  into  this.  Hi  is  is 
detailed  in  Figure  4.  The  sample  was  placed  atop  the  brass  rod  which 
extended  up  from  the  bottom  of  the  sanple  holder.  The  sample  and  rod 
had  approximately  the  same  diameter.  The  short,  a  solid  brass  cylinder 
glued  to  the  micrometer  drive  rod,  was  carefully  lowered  until  it  made 
contact  with  the  sanple.  Admittance  measurements  were  then  taken  with 
data  treated  according  to  the  program  COMP  given  in  Appendix  I. 

500  Mhz-2,000  Mhz 

In  this  frequency  rar ire  a  slotted  coaxial  line,  terminated  in  a 
circular  sample  and  a  variable  short,  was  used  in  measuring  the  null 
shift  and  VSWR.  This  data  was  used  to  calculate  the  sanple  admittance, 
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FIGURE  3*  s^p  for  curing  dielectric  constant  and  lor,; 

tangent  In  the  range  1 00-^00  Itiz. 
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«ach  was  used  with  tl*  parallel  plate  capacitance  relatlona  t 
calculate  u»  dielectric  extant  »l  leas  taaeent  of  the  aarple. 

derivation  of  the  imM«  1U»  stance  equation  la 
flow.  It  la  a  lndlar  to  t,*t  given  for  Irpolance  hy  Collin  Ml 

ntyrv  ,  0*»a  a  transmission  11*.  »»  *  load  V  *"1" 

incident  anl  roflected  voltages  V+  anl  V  respectively.  The  Incident 
aM  reflected  waves  Interfere  to  fom  a  sinusoidal  atandlnf  wave 
pattern.  Thus  the  voltagp  at  aw  point  on  the  line  can  he  written  as 

^  v  .  vV>»*  *  Vc-J»  tl2) 

„  .  21  is  the  propagation  constant  at  so*  firequency,  and  .  >  0. 
Similarly  Uie  current  at  any  point  Is 

i  -  iVc:  -  re"j6t 

TYie  nomadized  Input  admittance  Is 

yVBt  -  V~e“^Bt 

a  -  <$-)<$>  “  -T—m 

‘c 


(13) 


Y  ~T  (IA„) 

*ln  •  YT  '  tyfy  *  “» - »  ♦  V-e'^ 


(14) 


y  is  the  characteristic  advance  of  the  line.  By  definition 
the  reflection  coefficient  at  the  load  Is  rL  -  -f-  I*"8  (l4)  can  te 
rewritten  as 


0Jflt  .  r^e-JSt 
fln  “  JJbi  t 


(15) 


l’  can  also  be  written  as 
Li 


ZL  ~  ZC 
ZL  +  Zc 


•  a  well  known  identity  derived 

f  ** 


by  Collin  [21].  This  Identity  can  be  rewritten  as 


r 

B  '  Yc  +  \  1  +  7, 


(16) 


Now  replacing  by  (cos  «t  +  Jsln  M)  and  substituting  Equation  (16) 
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(17) 


(Ifi) 


into  (15)  we  obtain 

Y  eos  Bt  4-  Jsin  Bi 

Y  a  ■■■■■■  . .  — 

ln  cos  Bl  4-  JY^sin  Bt 

'llils  can  be  solved  for  Y^*  ttie  sample  admittance, 

Y.  cos  Bt  -  Jsin  Bt 

y  =  Jn - __ - 

L  cos  Bt  -  JYlnsin  Bt 

At  tSie  voltage  minimum  the  value  Yln  is  real  and  in  magnitude  equals  the 
VSWR,  as  proved  in  Appendix  II.  The  conputer  program  LL£/I  shown  in 
Appendix  II  uses  Equation  (18)  for  admittance  with  Equations  (9)  and 
(11)  to  calculate  values  of  dielectric  constant  and  loss  tangent  for  Un¬ 
samples  . 

The  equipment  used  for  this  frequency  range  included  a  slotted 
coaxial  line  (G.R.  Type  900-LB)  which  was  used  for  measuring  null  shift 


and  VSWR.  Hie  same  disk-shaped  samples  used  in  the  frequency  range 
100-500  Mhz  were  used  in  this  range  and  the  sample  holder  of  Figure 
was  used  at  the  end  of  the  slotted  line.  The  setup  used  is  shown  iri 
Fit ure  6.  The  frequency  was  measured  by  using  a  transfer  oscillator  and 
a  10  Mhz  Electronic  Counter  to  observe  a  harmonic  of  the  fundamental 
frequency . 

The  slotted  line  scale  was  calibrated  in  millimeters  so  that 
readings  were  accurate  to  0.001  meter.  The  VSWR  was  read  from  this 
scale  by  the  two  point  method,  since  high  VSWR  values  were  encountered 
[22].  This  method  involves  finding  the  position  of  a  standing  wave 
voltage  minimum  along  the  slotted  line  using  a  square  law  detector  and  a 
1000  Hz  meter.  Then  the  positions  of  points  on  each  side  of  the  minimum 
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FIGURE  6.  Equipment  setup  for  measuring  dielectric  constant  and  loss 
tangent  in  the  range  500-2,000  Mhz. 
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were  read,  which  tod  exactly  twice  the  power  detected  as  that  at  tto 
minimum.  Tto  VSVIR  ms  then  ccrputed  by  tto  two-point  formula 

V3WH  -  (19) 

where  AX  Is  tto  distance  betvf^cn  tto  twice  power  points  and  A  13  tto 
wavelength  in  tto  slotted  line.  In  this  frequency  range,  the  two-point 
method  is  valid  for  VfWH’s  greater  than  ten,  but  a  correction  can  be 
applied  for  smaller  values  [.22).  Tto  null  shirt  or  difference  In  mil 
losltlon,  for  tto*  sample  in  and  out  of  tto  circuit,  was  measured  near 
tto  center  of  the  slotted  line.  'Itose  three  pieces  of  data,  null  shift, 
wavelength  (from  frequency),  am  VSWR,  were  read  for  each  frequency  of 
interest  in  tto  range  '300-2,000  Mhz.  Tto  data  was  treated  according  to 
the  program  SHOT  given  In  Appendix  II.  Ihis  technique  was  not  used 
above  3,000  ffoz  be  ause  the  sample  arrangement  could  ther.  no  longer  be 
coiisiderxxl  to  be  a  lumped  capacitance  [23].  This  capacitance  technique 
was  used  successfully  by  Daly,  et  al  [24]  up  throutfi  che  frequency  2, ‘,00 
Hhx,  but  they  suggest  3,000  Mbs  as  tto  upper  limit  for  obtaining,  useful 
data  by  this  method. 

12,000  Khz-18,000  Mb- 

For  this  frequency  range,  a  K^-band  slotted  waveguide,  terminated 
with  the  sample  and  a  variable  short  circuit  was  used.  Tto  rectangular 
sample,  slightly  smaller  than  waveguide  dimensions  (0.622"  x  0.311"), 
was  placed  just  inside  the  variable  waveguide  short.  A  tightly  fitting 
styrofoam  spacer’  held  the  sample  erect  in  the  guide.  The  sanple 
thickness  was  approximately  0.22".  The  setup  used  is  shown  in  Figure  7« 
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FIGURE  7. 


Equipment  setup  for  measuring  dielectric  constant  and  lo 
tangent  in  the  range  12,000-18,000  Mhz. 
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A  variable  attenuator  was  used  to  decrease  the  reflected  nifTtu 
fron  Uw  short,  a  K-baixi  isolator  not  being  available  at  the  time  of 
testing.  A  vernier  scale  allowed  readings  to  0.1  nm,  but  table  and 
waveguide  vibration  caused  by  the  oscillator  fan,  caused  tho  probe 
position  to  move  slightly,  resulting  in  an  accuracy  of  about  0.«,  run. 

Measurements  taken  for  this  method  were  wavelength,  position  ol 
the  voltage  minimum,  art  VIJWR.  'lire  two  point  method  was  again  used  to 
obtain  the  VRWR,  readings  being  taken  near  the  center  of  the  waveguide 
to  avoid  end  effects.  A  broadband  probe  was  used  to  detect  the  signal. 
The  voltage  minimum  of  the  VSWH  pattern  was  made  to  appear  on  the 
lowest  scale  of  the  1000  Hz  meter  by  shortening  the  probe  depthi.  Boforo 
any  measurements  were  taken  for  glass  sanples,  a  number  of  teflon  strips 
were  used  as  a  sample .  These  strips  were  slightly  smaller  than  the 
waveguide  dimensions  and  were  pressed  together  by  the  strofoam  spacer. 
Measured  values  of  the  teflon's  dielectric  constant  averaged  within  0% 
of  known  values. 

■me  method  used  for  tills  frequency  range  is  derived  in  complete 
form  by  Roberts  [25]  and  by  Von  Hippel  [22].  It  involves  equating  the 
wave  impedance  of  the  air-filled  waveguide  with  that  of  the  sample 
filled  section  of  waveguide,  at  the  air-sample  interface.  By  using 
charts  [22]  which  gave  solutions  to  the  resulting  complex  transcendental 
equations,  tiie  propagation  constant  in  the  sample  was  calculated.  Ibis 
then  led  to  the  complex  permittivity  of  the  sanple, 
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e  = 


<r>2  -  (— )^ 

X-  '2ir  ^ 


,1_n2 

_ 

(i-)2  +  (\f 

XC  A 


(20) 


where  Ac  is  the  cut-off  wavelength  and  equals  twice  the  width  of  the 
waveguide.  The  dielectric  constant  and  loss  tangent  are  then  given  by 


e*  =  Re[e  ] 
tan  6  = 


Im[e  ] 
“  1  T~ 

Re[e  ] 


(21) 

(22) 


The  computer  program  HIPELKU  given  in  Appendix  III  uses  Von 
Hippel's  equations  to  detennine  the  wave  Impedance  at  the  air-sample 
interface.  This  value  was  used  with  the  available  charts  [22]  and 
allowed  the  determination  of  r2,  the  propagation  constant  In  the  sa^le, 
and  the  conplex  permittivity,  In  the  computer  program  FINALKU  of 


Appendix  III. 
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CHAPTER  III 

RESULTS 

As  described  in  Chapter  II,  dielectric  constant  and  loss  tangent 
data  were  taken  in  three  frequency  bands,  100-500  Mhz,  500-2,000  Mhz, 
and  12,000-18,000  Mhz.  The  data  thus  obtained  for  dielectric  constant 
is  combined  together  on  one  continuous  graph  as  a  function  cf  frequency 
in  Figure  8.  Figure  9  shows  the  combined  data  for  loss  tangent  versus 
frequency.  The  average  values  of  data  for  the  entire  frequency  range, 
100  Mhz-18,000  Mhz,  are  indicated  in  these  figures  by  a  dashed  line  for 
each  corrposition. 

The  variation  of  dielectric  constant  with  composition  of  the 
material  is  more  clearly  indicated  in  Figure  10.  Ihis  figure  shows  the 
average  value  of  dielectric  constant  taken  over  the  entire  frequency 
range,  plotted  against  the  five  compositions  of  Elass.  by 

increasing  the  tellurium  content  of  the  glass,  it  wa3  found  that  the 
dielectric  constant  can  be  varied  from  approximately  7.0  to  10.0,  for 
the  compositions  studied.  The  values  for  dielectric  constant  found  Ik? re 
conpare  favorably  with  data  recently  given  for  ASgSe^  by  Taylor,  et  al 
[26],  by  Crevecoeur,  et  al  [27],  and  by  Lakatos,  et  al  [28],  who  give  a 
microwave  dielectric  constant  ranging  between  9.0  and  9.7  for  this 
material . 

Figure  11  shows  the  average  value  of  loss  tangent  over  the  entire 
frequency  range,  plotted  against  the  five  compositions  of  ASgTe^-As^JJc ^ 
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glass.  Hie  loss  tangent  was  found  to  also  Increase  with  increasing 
tellurium  content,  and  varied  between  approximately  0.02  and  0.06. 
Crevecoeur  and  de  Wit  [27]  have  presented  the  loss  tangent  of  As^Se^ 

glass  over  the  frequency  range  10?  to  lo'1®  Hz.  In  the  U.H.F.  and 

-4 

microwave  regions,  they  give  an  average  loss  tangent  of  10  .  However 

other  investigators  [29,j0j  give  loss  tangents  for  this  material 
averaging  between  10"  1  and  10-3.  These  values  may  be  considered  as  a 
low  limit  to  the  data  for  loss  tangent  given  in  the  present  study,  since 
they  represent  glass  naving  no  tellurium  content. 

The  dielectric  constants  and  loss  tangents  of  the  samples  were 
found  to  be  essentially  frequency  independent  over  the  frequency  range 
studied.  Slight  variations  can  be  attributed  to  certain  measurement 
errors  and  to  differences  of  method  and  equipment  used  in  the  three 
frequency  ranges.  Measurement  errors  included  errors  in  slotted  line 
null  readings  caused  by  oscillator  and  waveguide  vibration  produced  by 
the  oscillator  fan,  errors  in  VSWR  measurements,  and  slight  errors  in 
sample  height  and  position.  These  led  to  a  statistical  error  of  about 
+  ten  percent . 
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CHAPTER  IV 

CONCLUSION 

This  paper  has  presented  the  U.H.F.  and  microwave  dielectric 
properties  of  the  glass  system  (x)As2Te2(l-x)As2Se2  for  x=0.i4,0.5,0.6, 
0.7,  and  0.8.  Future  research  performed  on  this  glass  should  include 
the  determination  of  the  behavior  of  the  dielectric  properties  as  a 
function  of  temperature,  and  under  switching  conditions.  This  could 
help  to  clarify  our  understanding  of  the  basis  of  the  switching 
mechanism  in  amorphous  thin  films. 

A  possible  use  for  the  glass  system  is  an  a  dielectric  substrate 
material  for  microstrip  transmission  lines.  Within  a  given  range,  one 
may  choose  the  composition  of  the  glass  for  specified  dielectric 
properties,  since  the  dielectric  constant  can  be  varied  from  7.0  to 
10.0,  and  the  loss  tangent  and  from  0.02  to  0.06.  'The  loss  can  be 
reduced  further  by  reducing  the  tellurium  content  while  still  retaining, 
the  amorphous  properties.  Conmonly  used  microwave  substrates  with 
dielectric  properties  in  this  region  are  Silicon  (e  *  11.5,  tan  - 
0.1b),  Yttrium  Iron  Garnet  (YIG) ,  and  AlSiMag  [30].  The  glass  system 
presently  studied  has  several  properties  which  are  desirable  in  choosing 
a  substrate  material  [31].  These  include  a  fairly  high  dielectric 
constant,  high  resistivity  [11],  essentially  constant  permittivity  over 
the  U.II.F.  and  microwave  regions,  and  a  good  surface  smoothness  with  no 
s'orface  ripples  compaiable  to  the  size  of  the  thin-film  components.  The 


25 


26 


glass  also  exhibits  a  high  dielectric  strength,  dependent  on  the 
magnitude  of  the  threshold  voltage  required  for  switching. 
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APPENDIX  I 


COMP 


Line 


5 


WRITE(6,200) 

READ  ( 5 ,  10 0 )  G ,  B ,  FREQ ,  D I A ,  T 

3  G=0.001#G 

4  B=0.001*B 

5  C=B/(6.2832*FREQ) 

6  A-0.7854*((DIA)«*2) 

7  CAIR=(0.0254#A)/(I13.0973*T) 

8  EP=C/CAIR 

9  TNDLT=G/B 

10  WRITE  (6 ,300) FREQ 

11  WRI'l’E  ( 6 , 4  00 )  EP ,  TNDLT 

12  100  FORMAT (§G) 

13  200  FORMAT(  'G(MMHO)  ,B(MMHO)  ,FREQ(GHZ)  ,DIA(IN)  ,TH(IN)  ’ ) 

14  300  FORMAT ('FREQ=’F7. 3,1* 'GHZ* ) 

15  400  FORMAT(  * EP= ' F10 . r'  . "" 

16  END 


o , 5X ' LOSS  TANGENT® ' F10 . 5/ ) 


line  program  COMP  uses  the  admittance  data  from  the  Thurston 
Bridge  to  calculate  the  relative  dielectric  constant  and  loss  tangent  of 
the  sample  In  the  frequency  range  100-500  Mhz.  Values  of  G  and  B  are 
entered  in  millimhos,  the  frequency  in  Gigahertz,  and  the  sample 
dimensions  in  Inches. 

In  deriving  the  equations  for  permittivity  shown  in  the  pr-og;-?un, 
Equation  (9)  io  first  rewritten  as 

B 

I  u> 

c'  -  7T  (23) 

°T 

■Ihic  is  scon  to  be  ttie  ratio  of  the  capacitance  measured  with  the  s;uuple 
in,  to  the  capacitance  calculated  for  the  sanple  out  (free  space 
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dielectric),  when  the  thickness  T  is  equal  In  each  case,  'me  numerator 
of  Equation  (23)  Is  shown,  on  line  5  of  the  pro^am.  denominator  is 

multiplied  by  0.0234  to  change  inches  to  meters.  c0  is  normalized  and 
is  written  as  £  -  njW  factor  1C9  is  dropped  n-om  and  is 

also  dropped  from  the  frequency,  which  is  expressed  in  Oigahertz.  Hds 
eliminates  very  small  numbers  from  computer  program  calculations.  Note 
that  the  109  factors  would  have  cancelled  in  Equation  (23)  even  if  they 
were  not  dropped.  The  resulting  equations  used  to  calculate  dielectric 
constant  and  loss  tangent  are  shown  on  lines  7  and  8  of  the  program. 
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APPENDIX  II 


SLOT 


Line 

1  COMPLEX  CJ,YL 

2  WRITE(6,200) 

3  read(5,ioo)dia,t 

4  WRITE(6,300) 

'j  READ(5,100)FREQ,X1,X2,B0,B1 

6  EmG=30./FREQ 

7  VSWR=ELMG/(3.l4l6*(X2-Xl)) 

8  YIN=VSWR 

9  BETA=6 . 2832/ELMG 

10  ELL=(B0-B1) 

11  A=BETA#ELL 

12  C-SIN(A) 

13  D=COS(A) 

14  CJ=CMPLX(0. ,C) 

15  YL=0 . 02* (YIN#D-CJ )/ (D-YIN*CJ ) 

16  G=REAL(YL) 

17  B=AIMAG(YL) 

18  C=B/(6.2832*FREQ) 

19  A=0.7854*#((DTA)##2) 

20  CAIR=(0.0254»A)/(113.0973HT) 

21  EP=C/CAIR 

22  TNDLT=G/B 

23  WR1TE(6,400)EP,TNDLT 

24  100  FORMAT (9G) 

25  200  FORMAT( 'DIA(IN) ,T(IN) ' ) 

26  300  FORMAT( 'FREQ(GHZ) ,X1,X2,B0,B1(ALL  IN  CM)') 

27  400  TORMAT( ' EP= ' F10 . 5 , 5X ' LOSS  TANGENT® * F10 . 5/ ) 

28  END 


The  program  SLOT  uses  the  data  from  slotted  line  measurement:; 
and  sample  dimensions  to  calculate  the  relative  dielectric  constant  and 
loss  tangent  of  the  sample  in  the  frequency  range  500-2,000  Mhz.  The 
VSWR  was  obtained  by  tne  two-point  method  (line  7  of  program)  wherein 
the  position  of  the  voltage  minimum,  Bl,  and  its  twice  power  points  on 
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either  side,  XI  and  X2,  were  observed.  BO  is  the  position  of  the 
voltage  minimum  with  the  sample  removed  and  the  line  short  circuited, 
and  (B0-B1)  is  the  "null  shift."  Values  of  slotted  line  data  were 
entered  in  centimeters,  with  the  frequency  entered  in  gagahertz  and 
sample  dimensions  in  inches. 

The  value  Y^n  of  line  15  is  the  input  admittance  at  the  position 
of  a  voltage  minimum.  This  value  is  therefore  real  and  maximum  and  can 
be  written  -  Of^.  For  the  position  of  a  null,  Fiquation  (16)  for 
reflection  coefficient  can  be  rewritten 


1  -  G 


max 


(24) 


1  +  G. 


max 


By  definition  the  voltage  standing  wave  ratio  is 

3  -  lv+1  +  IV~l  -  1  +  I (V~/V+) I  -  1  + 

|V+|  -  |V"|  1  -  | (V”/V+) | 

Substituting  (24)  into  (25), 


1  - 


(25) 


1  +  G  +  1  -  G 
c  _  max  1  max' 


1  +  G 


max 


-  1  -  G 


(26) 


max' 


Since  G  _  is  greater  than  1,  then  1  -  G„  I 

max  ’  1  max 1 


be  rewritten 


=  G. 


max 


Gmax  -  1-  (26)  may 


(27) 


which  states  that  the  VSWR  equals  the  magnitude  of  the  normalized  Input 
conductance  (or  admittance)  at  a  null.  Ti.i^  result  is  shown  on  line  b 
of  the  program. 

'Hie  load  admittance  is  calculated  on  line  15  and  is  denormallzed 
by  multiplying  by  0.02  mho.  Lines  16-22  of  the  program,  which  calculate 
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the  relative  dielectric  constant  and  loss  tangent  of  the  sample  in  a 
capacitor  arrangement,  are  identical  in  form  and  method  to  those  given 
in  Appendix  I  in  the  program  COMP. 
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APPENDIX  III 
HIPELKU  AND  FINALKU 


Hipelku 


Line 


L4 


COMPLEX  B1,CEJZ,ELMG1 
P£AD(5,100)D1 
D-2.5^*D1 

BEAD( 5 , 100 ) XI , X2 , XI 1 , X22 
AL2=(Xll+X22)/2. 

ALl=*(Xl+X2)/2. 

ELMG=2.*(AL1-AL2) 

X0-AL2+0.56-D  . .  ^  x> 

VSWR-ELMG/(3 . 3U36*(  ( (X1-X2)+(X11-X22)  )/2. ) ) 

RECIP-l./VSWR 
A=6.2832»X0/ELMG 
I>SIN(A)/COS(A) 

BlaCMPLX(0.  ,13) 

CEJZ-{*C  1/(6^332*0)  ),(RECIP-B1)/(1.-(RECIP*B1)) 

C-CABS(CEJZ) 

RBCPC-l./C 
P-HEAL(CEJZ) 

G-AIMAG(CEJZ) 

ZRA-ATAN2(G,F) 

Z*ZRA*l8Q./3.1^l6 
WRrrE(6,200)C,Z,RECPC 

i  200  !ST'C-'F9.'l,bX'Z(DB5liaiS)-'P12.3,^'ONIi  CM®  CF9.V) 
>h  END 


The  program  HIPELKU  uses  the  experimental  data  for  the  frequency 
•ange  12, 000-18, OCr  :5b,  with  line  15  of  the  program,  to  calculate  Urn 
jarametero  C  and  Z  of  the  complex  number  Ce^2.  »  oaaple  ;hicknes.,  Di 


la  entered  in  inches,  while  the  slotted  waveguide  readings  arc  enter- 1 
in  centimeters.  The  readings  XI  and  X2  are  twice  power  points  on  either 
side  of  the  null  AL1,  the  same  being  true  for  m,  X22,  and  A12.  In 
line  8  the  position  of  ;he  voltage  null  (X0)  with  respect  to  the  sanplc 
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face  IB  calculated.  The  quantity  (0.56-D)  centimeter  la  added  to  this 
line  because  the  shorting  plunger  was  fixed  at  a  position  0.56 
centimeter  Inside  the  south  of  the  waveguide  short,  and  the  sample,  of 
thickness  D.  was  placed  In  this  space,  n*  VSWR  was  calculated  by  the 

two-point  method  as  3hown  on  line  9« 

me  values  C  am  Z.  calculated  on  limes  16  and  21.  are  entered 

into  tables  which  allow  one  to  read  T  and  t.  the  parameters  of  the 
cooplex  number  Te^x.  T  and  t  are  then  entered  In  the  program  mm 
shown  below,  which  calculates  the  dielectric  constant  and  loss  tenant. 

Finalku 


6  t 

7  < 

8  i 

9  I 

10  i 

11  ! 
’2 

J-3 

1U 

19 

16 

17 

18 

19 

20 

21  100 
22  200 

23  300 


COMPLEX  TAURC,GAM2,EPSLN 

HEAD(5,100)H,D 

REA1)(  9 , 100 )  AL1 ,  AL2  ,T  ,TAU 

TAUH^rAU*  3  •  I1* 159/180 . 

TAURC-CMPLX(0. ,TAUR) 

DCM*2.5^#^ 

0AM2*(T*CEXP(TAURC) )/0CM 

A*0.622 

B-0.311 

ACM“2.91<*A  <x 

KL1P"1 ./(*< .*ACM*ACM) 

^(^-(S/6.2832)»2)/(FUPe(l./EL^)«2) 

041  '*RKAL(  EP3LN ) 

IvMPP-AIMAG  ( KPSIH ) 
lW)LM*01PP/i'MP 

ECF-tMF*H/Cls-( 

WDLC1WU<«U/(B-((b-H)W)) 

W!<m(6,200)i-i:p,'ir«LC 

WhriE(6,300)B*1P,TNDUi 

l  JSSS^Wnuc  mrtcoRRi-m.T.sx.-usi 

0  'UliaiCTWC  CrXt;rAHr(KtA3)-'F11.7.5X,’LOSS  TANO(wae)-' 

K11.7/) 

BW 


211 


3* 


the  Input  data  fur  '.his  prvy^T&fi  Include  sW*|de  bright  ant 
thickness,  II  and  D,  which  -ire  entered  In  in:!**,  slotted  line 
itKasurcrsenia  In  cenilserters,  are  values  for  T  and  i  obtained  fit**  starts 
.uii  the  proven  hlltxKU.  The  dielectric  constant  and  loss  tanrtnl  ev 
calculated  usli>*  equations  (20)  through  (22),  In  lines  13*16  or  HIlMd-W). 
these  values  are  Ms*n  corrected  by  using  fo tttumz  given  Uy  Wlwl  awl 
Itspuport.  L  33  J  which  correct  c*  tw*i  tan  A  for  the  air  gap  between  u#/ 
aarple  and  tbs  bixadwall  or  tlw  waveguide.  The  corrected  values  for 
dielectric  constant  and  loss  tangent  are  calculated  on  lime  1?  awl  lb. 
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